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XXX. BAND SECHSTES (SCHLUSS-) HEFT 1961 


Note sur les vitesses critiques d’un arbre muni de disques 
Par J. Kestens 


1. Introduction. Nous nous proposons d’examiner les vitesses critiques d’un arbre muni de 
disques. Nous supposerons que l’arbre est de révolution autour d’un axe rectiligne (dans sa position 
non déformée). Des disques ou volants sont calés sur cet arbre. Considérons un de ces disques, 
le i-me p. ex.: G; est son centre de gravité et J,;, Jo; = Js; sont les trois moments d’inertie prin- 
cipaux. Ce disque est calé sur l’arbre de telle facon que l’axe principal portant J); est paralléle 
a l’axe x de l’arbre non déformé, ce qui revient 4 supposer que l’angle de calage du disque sur 
Varbre est nul. L’intersection de l’axe x de l’arbre avec le plan principal perpendiculaire a l’axe 
portant J,; définit un point A; qui est en général différent de G;: les disques sont imparfaitement 
centrés. Soient ¢; = A; G; l’excentricité et O; l’azimut de A;G; mesuré A partir d’un repére com- 
mun a tous les disques. 

Le rotor est donc formé d’une série de n disques calés sur un arbre élastique de masse linéaire 
variable et de moment d’inertie géométrique des sections variable. Le probléme dynamique général 
d’un tel rotor en rotation est complexe. Nous simplifions le probléme en négligeant les vibrations 
de torsion d’une section de l’arbre par rapport 4 une autre section et en supposant donc que ce 
rotor est animé d’un mouvement de rotation en bloc. L’azimut du disque i en rotation est donné par 


V=~vtO0,, avee p=a, 


ou @ est la vitesse instantanée de rotation. 

Considérons un instant, avec Grammel [1]* un arbre sans masse et supposons ¢; = e?/k? <1 
avec m; k? = Ji; (m;: masse du disque i, k;: rayon de gyration). Lorsque ce rotor tourne a une 
vitesse w = const., l’application des théorémes de la résultante cinétique et du moment cinétique 
montrent qu’en premiére approximation (parce que ¢; < 1) le mouvement de G; est solution d’un 
systéme de deux équations différentielles avec seconds membres, soit donc la solution générale 
(somme de mouvements périodiques de pulsations q;) et une solution particuliére de pulsation w 
(avec w+ w;). Pour que w = 0 rigoureusement il faut appliquer 4 l’arbre un couple M qui est 
petit et en premiére approximation M = 0. Lorsque w = @; est une des pulsations solution des 
équations différentielles sans second membre, la solution particuliére croit linéairement avec le 
temps. Cette croissance est trés rapide aussi ces vitesses de rotations sont-elles trés dangereuses, 
c’est la raison pour laquelle on les appelle les vitesses critiques ([1]. t. 2., p. 155—159). es 

Grammel [1], [2] a établi ce qu’il appelle le Principe d’ Equivalence: Les vitesses critiques @; 
coincident avec les vitesses de rotation d’un systéme équivalent ow les excentricités e; sont nulles 
(systéme parfait) et qui est tel que cet arbre ait une fléche stationnaire dans son équilibre relatif, 
sous l’effet des réactions d’inertie d’entrainement (forces centrifuges et couples gyrostatiques) et 
des forces de rappel élastique. E vat 

Grammel a démontré cette proposition lorsqu’on ne tient pas compte de Peffet raidissant du 
couple gyrostatique du disque ({1], t.2, p. 191). Ziegler a étudié le méme probléme en tenant 
compte du couple gyrostatique et est arrivé aux mémes conclusions ([8], n 2, 3 et 4). Dans un 
ordre d’idées plus générales, Ziegler [7] a étudié de fagon critique les divers critéres de stabilite 
(méthode statique, du défaut initial, énergétique, dynamique) et il a montre que le principe d équi- 
valence doit étre utilisé avec prudence, car il n’est pas valable d’une fagon absolue. En ce qui 
concerne le probléme que nous traitons ici, le systéme parfait est stable quelle que soit la vitesse 
de rotation w mais il ne posséde de configuration non triviale que pour certaines vitesses de rotation 
@; Ce sont ces configurations non triviales (fléches stationnaires) qui, par lexistence d’imper- 
fections (balourds, e; + 0) eréent des forces d’inertie entrainant une instabilité ((7], P- 396). ue 
contre, le systéme réel (¢; + 0) est stable pour w == «; et le principe d’équivalence est bien d’appli- 

cation dans notre cas ({7], n° 31 et 32, p. 394397). 


~ * Voir les références page 384. - 
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Ces vitesses critiques sont également appelées vitesses critiques synchrones directes (,,Gleich- 
lauf‘’). Ce sont les seules vitesses critiques qui sont excitées par les forces d’inertie des balourds. 
D’autres forces excitatrices, telles que des vibrations transversales des appuis, peuvent exciter 
d’autres vitesses critiques, dites vitesses critiques rétrogrades (,,Gegenlauf™, [8], p- 383). Ces 
vitesses critiques sont beaucoup moins dangereuses car elles impliquent une importante dissipation 
d’énergie interne de l’arbre qui ne tourne plus d’un bloc ((1], t. 2, p. 164). 


Fig, 1. 


Rappelons enfin que nous supposons que l’arbre est de révolution; en fait, ce que nous exigeons 
est que l’inertie et le coefficient a(x, s) (voir ci-dessous, n° 2) ne dépendent pas de l’azimut. Si 
par exemple, a(x, s) ne présente pas cette symétrie (influence de clavettes de fixation des disques 
par exemple), Grammel a établi qu’il y a dédoublement des vitesses critiques w; et que le systéme 
est instable dans les intervalles des vitesses qui en résultent ({1], t. 2, p. 166—168). 

Le principe d’équivalence permet de simplifier la formulation mathématique du probleme qui 
nous occupe: établir l’existence des vitesses critiques synchrones directes (,,Gleichlauf*). 


2. Coefficients d’influence. Soit F x(x, s) la fléche de arbre en x due a une force F appliquée 
en s. L’application du théoréme de Betti montre immédiatement que a(x, s) = a(s, x). Outre le 
coefficient d’influence a, introduisons les coefficients f, y et 6 définis de la fagon suivante (figure 1): 


Bee da(x, s) io a(x, s+ As) a(x, s) 


= fléche en x due a un couple 


ds As As eens 
0a(x, Ss) ect = 
WS) es a lim — [a(x + Ax, s)— a(x, s)] = pente en x due a une 
be Ax 
force + l ens, 
67a (x, s) 


6(x, s) = 


Duds 7 Dente en x due a un couple + lens. 


Fig. 2. 


Ces coefficients d’influence ont été introduits par van den Dungen [3] (voir également [1] t. I 
p- 95). Ils jouissent de propriétés de symétrie que nous allons établir briévement. L’arbre, ou 
le rotor, est de longueur / et la fonction a(x, s) est définie dans un carré de cotés I. Si P et P, sont 
deux points du carré symétriquement placés par rapport a la diagonale x = s (figure 2), on a, 
comme il a été rappelé ci-dessus, 


W(x, s) = a(s,x) ou a(P)=a(P,). (2) 


Si on porte la valeur de la fonction «(P) perpendiculairement au plan x, s on obtient une surface 
représentative de la fonction «, qui est symétrique par rapport au plan x = s et on en déduit 


53 0 P) = = 0 Py) 


soit 


y(P) = B(P;) ou (x, s) = Als, x), (3) 


(x, s) = 6(s, x) . (4) 


et on a évidemment 


4 
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3. Equation de \’équilibre relatif du rotor. En vertu du principe d’équivalence, exprimons que 
le rotor en rotation est en équilibre relatif. A cet effet, considérons un disque. Soient (figure 3) 
Ox la ligne d’arbre (axe du rotor non déformé), Oy perpendiculaire 4 Ox et tel que l’arbre fléchi 
soit dans le plan (x, y), Oz complete ce triédre tournant autour de Ox a la vitesse angulaire w. 
Attachons en outre au disque considéré le systéme d’axes (G; €, 7,2), le vecteur unitaire i est 


tangent a la ligne élastique, 1, est dans le plan (x, y) et 1, est paralléle a l’axe Oz. 
Le théoréme de la résultante cinétique appliqué a ce disque 
donne Z 


0=mj,= Y F—mj,—mj, (5) 


ou j,, J, et j, sont les accélérations relatives, d’entratnement et de Y 
Coriolis et 3) F est la résultante des forces extérieures appliquées 
au disque. 

On a j, = — w’ y, j, = j, = 0 et (5) donne 


Se AS ree y . (6) 


, 


ou — )/ F est, en vertu du principe de l’égalité de l’action et de la réaction, l’action du disque sur 
larbre. 

Pour appliquer le théoréme du moment cinétique, commengons par calculer ce dernier. Le 
disque est soumis 4 une rotation instantanée absolue 


On =O, Ov— 0. 
puisque w, = 0 (équilibre relatif). 
Le moment cinétique Wg (= A p1e+ Bq1,+ Cr1, suivant les notations classiques) vaut 
Me = J; w, cos y 1,— Js @, sin y 5 : 
Par dérivation on a 


ANE ms dle O dl, 
Ha J Pe 608 PF — Jaw, sin p 
et 
dl = dl 7 
= sin ~ w, 1,, Fr = COS P We 1,, 
donc 2 
=e = (J, — J) w’ sin ¢ cos Thee 


La pente @ étant faible, sin px yw y' et cosgx 1. Posons C= J,—J,). Ona 


= Cw’ y'l, 


= couple appliqué au disque (= 0) + couple de réaction de V’arbre sur le disque (qui vaut — couple 
de réaction du disque sur l’arbre). 
L’action du disque sur l’arbre est donc un couple 


BC gsi 1, (7) 
Les équations de l’équilibre relatif de arbre s’obtiennent immédiatement avec (1)57(6)"e0 (7): 
Varbre est soumis aux forces (6) et aux couples (7), on a donc 


y(x) = w* » a(x, s;) m(s;) y(s;) — o 2 B(x, si) C(si) ¥'(s) » 
y' (x) = @? BD y(x, s;) m(s;) ¥(s;) — oo 2 d(x, s;) C(s;) y"(si) -* 


3 bak iti i i i jante ona C & J,/2 & Ja). 

1 C est généralement positif dans le cas des disques (pour les disques minces sans Jan x/ 
La théorie ae nous développons est également valable pour des disques tels que C < 0 (Ceci pourrait se pro- 
duire si le disque est muni d’une jante importante). ce ' ie 

2 Dans le as des vitesses critiques synchrones rétrogrades (,,Gegenlauf*’), on a @, = — Wg; soit a 2 Oa. 
Il est aisé de voir que dans ce cas on obtient les équations aux vitesses critiques en remplagant | (si) par 
— C'(s;) = — [Ju(si) + Jo(s;)] (LL p- 186). Enfin, le cas des oscillations de flexion du systéme or vibre a 
un plan, sans tourner, s’obtient en remplagant C(s;) par — J,(s;). @ est, dans ce cas, la pulsation de ces oscilla- 
tions ([6]). - 

({61) 
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Si nous nous limitons a l’approximation de la résistance des matériaux, on sait que les sections 
de l’arbre initialement planes restent planes aprés déformation (le gauchissement de ces sections 
est di aux contraintes tangentielles, nous négligeons cet effet). Dans ces conditions l’arbre lui- 


s 


méme peut étre considéré comme la juxtaposition de disques d’épaisseur ds dont la masse vaut 
u ds et dont la différence des moments d’inertie principaux vaut ¢ ds = j, ds —j, ds = TH R? ds> 0 


avec R, le rayon de l’arbre. L’équation générale de l’équilibre relatif s’écrit donc 
I 
ye) = oF Sale, &) miss) 1 of aCe, #) Bis) Xe) ds 
l 
— oF 5: Ble, 51) Cle) ¥'(si\— oo | Ble, 8) fs) ¥ (5) ds 


et une équation analogue pour y’(x). On peut écrire ces équations sous la forme suivante: 
y(x) = w? S a(x, s) m(s) y(s) — w* S B(x, s) e(s) ’(s) » 


¥'(s) = 08 S yl, 8) m(s) y(s) — 0 d(x, 8) e(8) ¥(s) ®) 


ou les sommes S sont étendues sur les masses concentrées et continues, ce sont des intégrales de 
Stieltjes, déja considérées par B. de Saint-Venant. 

Nous obtenons ainsi un systéme homogéne de deux équations intégrales pour les fonctions 
y, y, ou encore, une équation intégro-différentielle [une ou l’autre des équations (8); la seconde 
est en effet la dérivée par rapport a x de la premiére si on se souvient de la définition de +(x, s) et 
6(x, s), (n® 2)]. 

4, Réduction du systéme (8) a une équation intégrale 4 noyau symétrique. Multiplions la premiére 
équation (8) par )/m(x) et la seconde par i }/c(x), i= —1. Le systéme (8) s’écrit, en posant A = w?, 


V/ma(x) (x) = AS a(x, s) V(x) m(s) Y/m(s) y(s) + AS i B(x, 8) Ve(s) m(x) i Ves) A ‘s) 
i Ve(x) (x) =A S i-y(x, 8) m(s) e(x) /'m(s) y(s) —A S A(x, s) Ve(x) e(s) t Ve(s) ¥'(s) - 
Effectuons le changement de variables 
€ =x-+ pl=a2l, D= Vets (10) 
nH=stpls2l, p =0,), 


et posons 
ie: , Fmt) a) pour p= 0, soit pour 0X &€< 1, (11) 
i Ve(x) y’ (x) pour p=/, soit pur I< F< 21. 
Le systéme (9) s’écrit 
21 
Y() =A S, K(E,0) Y(o). (12) 
ott le noyau K est défini de la facon suivante | 
Ko) | O<éCI | l<é<21 
| 
0<o<l | x, s) Ynz) m(s) g(x 8) Prml sy ola) ae 
l<oK<2l i B(x, s) Ym(x) (8) | = —8(x, s) Ve(x) (s) 


En vertu de (3) on constate que le noyau K est symétrique: K(é,c) = K(o, &), il est complexe 
(mais pas hermitique). 

L’équation (12) est donc une équation intégrale homogéne mixte noyau symétrique complexe 
Nous montrerons que cette équation admet des solutions Y non triviales (dites fonctions propres) 
pour des valeurs discrétes de A (dites valeurs propres). Y(&) défini par (11) est réel pour 0 <E< | 
et (lorsque ¢(x) > 0) imaginaire pur pour 1< & < 21 si le systéme (8) ou (9) admet une solution 
réelle (y, y’) (ce que nous montrerons étre le cas). . 
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5. Orthoganalité des fonctions prepres. Supposons que l’équation (12) admette les valeurs pro- 
pres A, + A, auxquelles correspondent les fonctions propres Y, et Y,: 


Y,(¢) =4,S K(&,o) Yi(c) et Y,() =A, S K(é,c) Y,(0). 
Ona ‘ : 
2 Y€) YE) = 4, $ S KE, ¢) ¥,(0) 26) 
=4,SS K(é,c) Y,(c) Y4(6) 
=A, S S K(&, 0) Y,(c) Y,(é) (14) 
puisque K(é, oc) = K(o, é), done 
0 = (A, —A,) SS K(é, 6) Y,(c) Y,(6), (15) 
et puisque A, =: /,, on obtient 
| SS K(, 0) Yi(o) Y2(6) = $ Yi) ¥2(4) = 0, 
soit explicitement ; 


(16) 
S m(x) ya(2) ya() —$ ef) y4(x) yx) = 0 


qui exprime que les fonctions Y,, Y,, soit (y1, y{), (yo; Vs), sont orthogonales. 


6. Réalité des valeurs propres. Supposons que le systéme (9) admette une valeur propre com- 
plexe A a laquelle correspond une fonction propre complexe y= u--iv, y’ =u’ +iv’. Les 


coefficients du systéme (9) étant réels, ce systtme admet également la valeur propre conjuguée A 
et la fonction propre conjuguée y = u—iv. Posons 


Y,(6).—= | 'm(zx) y(x) fonction propre a laquelle correspond 
NEO rey ee ae la valeur propre A, 
i Vela) y’ (x) ae 
a \/m(x) y(x) fonction propre a laquelle correspond 
; j V/e(x) y' (x) la valeur propre A, } 
et appliquons (15) rs 
0 = (A—A) SS K(é, a) Y,(c) Y,(é). (18) 


On a 
SS K(E, 0) Yi(o) ¥o(E) = SS [a(x, s) m(x) m(s) y(s) ¥(x) + i y(x, 8) m(s) e(x) y(s) 7 ¥(®) 
+ 4 B(x, s) m(x) e(s) t y'(s) (x) — d(x, s) e(x) e(s) t y'(s) t y'(*)] 
= SS ja(x, s) m(x) m(s) y(s) y(x) + (x, s) (x) e(s) ¥’(s) ¥'(*) 
— B(x, ) m(x) e(s) [y’(s) yx) + ¥'(s) y(#)]} 
= $ S {a(x, s) m (x) m(s) [u(x) u(s) + v(x) v(s)] 
+ A, 8) o() o(8) [w/(x) w’(s) + ve) 0'()] 
— 2 B(x, s) m(x) c(s) u(x) w'(s) + v(2) v'(9)]} 
en vertu des relations de symétrie (2) (3) (4). 


Ce second membre est formé de sommes réelles. Il s’ensuit que (A—A) = 0 et A est réel, et 
considérant le systéme (8) ou (9), les fonctions y(x), y’(x) sont également réelles. 


7. a) Définition. Considérons une fonction arbitraire 
m(x) u(x) , 0O<é<l, 
oe 
i V/c(x) w’(x) , Ve Sat, 
et u, u’ réels, suffisamment réguliére pour que 
SS K€, «) U(g) U(e) 


existe. Nous dirons que cette fonction U est admissible. 


1 On remarquera que Y,(€) = — Y,(&), Y2(&) n’est donc pas la fonction conjuguée de Y,(é). 
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b) Lemme 1. Le nombre SS K(&, 0) U(&) U() est positif. 
En effet 


SS K(E,0) UlE) Ulo) = $$ fax, 5) mls) m(s) u(x) w(s) —Blx, 8) m(x) e(8) w(x) w’(s) 
| — yl, 8) e() m(s) u’(x) u(s) + d(x, 8) e(2) e(s) w’(x) u’(s)] - 
Des forces F(s) et des couples + M(s) donnent les fléches u(x) et les pentes WAx): 
u(x) = S a(x, s) F(s) + S B(x, s) M(s) , 
u’(x) = S (x, s) F(s) + S$ 6(x, s) M(s) . 
Le travail (virtuel) de ces forces et de ces couples est positif. Il vaut 
S F(x) u(x) + S M(x) u(x) = S S [a(x, s) F(s) F(x) + B(x, s) M(s) F(x) 
| + ple, s) F(s) M(x) + (x, s) M(s) M@)] > 0. 
Il suffit de prendre les ,,forces‘‘ m(s) u(s) et les ,,couples“ -+- M(s) = — e(s) u’(s) pour que le lemme 


soit démontré. 
Ceci est une démonstration physique. Nous montrerons ultérieurement (n? 12) que 


l 
SS K(é, co) U(é) U@) = f E Iv’? dx 
0 
ou v(x) est la fonction itérée de la fonction u(x). Il est manifeste que cette derniére intégrale est 
définie positive. 
8.a) Définition. Soit U,(é) une fonetion admissible, la fonction itérée V,(&) est définie par 


V(é) = S K(é, 0) U,(o) - (19) 


b) Existence des Valeurs Propres. Théoréme I. Calculons le nombre réel R = R(Up, Vo) 
(rapport de Rayleigh) défini par 
__ SS K€, oe) Ug(a) Uol(E) _ S U(é) Volé) 


RU Yo) = SVG) SS K(E0) Ualo) Pal) os. 
Suivant que R(U,) est > 0 ou < 0, ona 
0< 1, = RU) oa 0S 427 Res (21) 


En vertu du Lemme | (n° 7. b), R(U,) a le signe de 
S Vo(E) = $ m(x) vo(x) —$ e(x) v9%(x) . 


Nous commencons par démontrer le 


c) Lemme II. II existe des fonctions itérées V(&) telles que SV2(é) + 01. (22) 


1. Cas o& u(x) = 0 ou e(x) = 0. C’est en particulier le cas ot il y a une masse répartie. Toute 
fonction itérée v est une fonction de classe C’ [a, b] (c.a.d. continue ainsi que sa premiére dérivée, 
et dont la dérivée seconde est continue sur des segments, dans l’intervalle fermé [a, b] limitant 
larbre) et qui satisfait les conditions aux limites aux appuis (soit, par exemple si l’arbre est simple- 
ment appuyé, v(a) = v’(a) = 0 et v(b) = v’’(b) = 0) comme cela ressort de la définition (19) des 
fonctions itérées qui sont les fléches que prend Varbre sous l’effet des charges ou des couples arbi- 
traires, ces fléches sont de classe C’ puisqu’elles peuvent présenter des discontinuités de moment 
fléchissant (v’’) et d’effort tranchant (v’”’). 

Supposons (démonstration par l’absurde) que S V2(é) = 0 quel que soit V,. Si V, est une 
fonction itérée, V, + V, Vest également et 


S[Vi(é) + Vi(é)P = S Ve(é) + S VZ(E) + 2 S V,(é) V,(€) = 0, 


soit 
21 l 
S,ValE) ValE)= S,max) 04x) va(x) — $ of) v{(x) v(x) = 0 (23) 


? Le mode de raisonnement de cette démonstration est di A Kamke [4]. 
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quels que soient V, et V,. Or une intégration par parties donne 
l 


1 I 
J ve v})' dx = [v,(c v4]. — fev, oj dx = — fev v, dx 
0 6 6 


ou les termes intégrés disparaissent puisque v est continu, et que nous supposons la poutre appuyée 
a ses extrémités. Donc 


1 
Smv,v,—Sev v= fv, (ur, + (cv) ]dx + Y mv,v,— Cv, v4 (24) 
0 


quelle que soit la fonction v,, les sommes étant étendues sur les disques. 


En particulier, pour v, telle que 7 m V1 Vg — J} Cv, v, = 0 (de telles v, sont manifestement des 
fonctions itérées, il est en effet possible de solliciter l’arbre par des forces et des couples tels que 
les fléches et les pentes aux droits des disques soient nulles) on a 


[vl rs + (cvi)'] de = 0 (25) 


quelle que soit la fonction itérée v,, ce qui implique jv, + (c vj)’ = N(v,) = 0. 
En effet, supposons N(v,) == 0, il existe alors un intervalle [«, 6] (ne contenant pas de disque) 
tel que N(v,) > 0 et on peut choisir v, > 0 dans cet intervalle, et v, = 0 en dehors de cet intervalle, 
I 


et dans ce cas { v, N(v,) dx > 0, en contradiction avec (25). 
0 


Ainsi il faut N(v,) = 0 pour toute fonction itérée v,, or il est manifeste qu’il existe des fonctions 
itérées v, telles que N(v,) == 0. En effet, si c == 0, ou si w == 0 (dans le cas ou on supposerait c= 0) 
il existe un intervalle [a’, 6’] tel que N(v,) = 0, on peut prolonger v, en dehors de [a’, 8’] de fagon 
a ce que les conditions aux appuis soient satisfaites et v, est une fonction itérée admissible. Il y 
a donc contradiction, ce qui démontre le lemme. 

2. Si uw = c = 0, reprenant la démonstration par l’absurde, on a 

SVEE) = > mF Co = 0, (22’) 
les sommes étant étendues aux disques. Or on peut choisir une fonction itérée telle que v(x;) = 0, 
v (x;) = 0, v(x;) = v'(x;)) = 9, 7 =1,2...i—1,i+1,...,m ot les x; sont les abscisses des n 
disques, sauf le disque i. Cette fonction est manifestement une fonction itérée admissible pour 
laquelle (22’) n’est pas satisfait. Cette contradiction démontre le lemme. 


d) Démonstration du Théoréme I’, Ayant choisi la fonction admissible U,(&), on définit la 
fonction NV par 


Uo(é) — R(Up) S K(E, o) Up(o) = N¢E) (26) 
avec 
| Ve@inG),- 0=6<1, 
life’), 1 S21. 
Si N(é)=0, le théoréme est démontré. Supposons done N(é) = 0 et considérons le probléme 
auxiliaire défini par l’équation intégrale non-homogéne dépendant d’un paramétre réel A 


U(é) —1 S K(&,0) U(o) = N@). (27) 


N(E) est la fonction définie par (26) et U(&) = U(§; A) dépend du paramétre A. Cette fonction U, 
solution de (27) existe tant que A n’est pas une valeur propre de l’équation sans second membre 
(car dans ce cas, il n’y a pas de solution a l’équation (27) puisque le premier membre serait = 0 et 
le second membre =¢ 0). A la fonction U(&; A) associons la fonction itérée V(&; 1) = S K(é, oc) U(o;A). 
En particulier, pour A = R(U,) on a U = U, et V= Vo. 
Considérons la fonction auxiliaire H(A) définie par 
H(A) = $ NE) ViEsA) =$ UEsA) ViEsA) —1$ SK(Es0) Us) ViEs2). 
Pour A = 0, 
H(0) = S U€é; 0) Vg; 0) = $$ K(é, «) U(E, 0) Uo; 0) > 0 en vertu du lemme 1 (n° 7.b). 
& o 


N(é) 


1 Cette démonstration est une adaptation de la démonstration du méme théoréme établi par Kamke [4] 
dans le cas du probléme auto- adjoint défini par une équation différentielle avec ses conditions aux limites, 
mais sans condition de passage [n° 12, équation (46)]. BY 
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Pour A = R, 
H(R) = $ NG) VolE) = $ Ul) Vole) — B SS K(&, c) Up(o) VolS) = 0 
en vertu de la définition (20) de R. 
Calculons la dérivée de H(A). 
H(z) — Ha) = $ [V(E5 42) — VEs i)] NEE) 
= $ V(Es As) [UEs a) —Ar S K(E; 0) U(o341)| 
—$ Vigsda) [Us te) —Ay $ K(E, 0) UlO; 4] 
= $ [V(E5 42) UEs ay) — V(Es Ax) U(Es As) 
+A, ° S K(&, 0) U(o34,.) V(EsA,) —A : S K(&, o) U(o3A,) V(E3 As) 
= $S KE, 0) [UE Ay) U(o3 22) — UlEs 2) Ulos Ay)I 
+ As S V(E5 Ay) V(Es Aa) — Ay $ VUE 3 Aa) VES Ay) 
= (lo —h) § VE) VSI) 
car K est symétrique, d’ou 
H(i) =$ V4E32). 
1). Supposons que R(U,) > 0. 
On a H’(R) = S VE; R) = S V2> 0 et comme H(R) = 0, H < 0 a gauche de R, soit 
H(A) = SS K€, 0) U(E3A) Ulos2) —AS V7(Es2) <0 
a gauche de R, ce qui implique, puisque J > 0 
S V?(E;2) > 0 a gauche de R. 


Il faudrait donc que H(0) < 0 s’il n’y avait pas de valeur propre entre 0 et R, or H(0) > 0, 
il existe donc une valeur propre A, telle que 


0<A4, S R(U,) 


etd, > 0, car sid, = 0, il faudrait que le numérateur de (20) soit nul, ce qui n’est pas le cas (lemme 1), 
2). Supposons que R(U,) < 0. Un raisonnement analogue montre que dans ce cas 


OS Aa 2 RU). 


Le théoréme est démontré: il existe au moins une valeur propre réelle 2, (ou A_;) A laquelle est 
associée la fonction prepre 


Yani =| mere 0SE<1 


7 (+ 1 pour A et —1 pour j_ 
iVe(x) yui(x), I< E<21 1 P i) 


soit 


Ye 1(x) = Ay S a(x, s) m(s) ¥+1(s) —Ax1 S B(x, s) e(s) ¥4.3(s) , 
(28) 


yale) = Asca S px, 8) m(s) ¥4.1(8) Ar $ (x5 8) es) y4ar(s) 


e) Corollaire. S’il existe une fonction itérée V, issue d’une fonction admissible U, telle que 
S V2(é) > 0, il existe une valeur propre /, définie par 
A, = minimum R(U, V) (29) 
(U, V) 
quand on considére l’ensemble des paires de fonctions (U, V) telles que S V2(E) > 0. 
De facgon analogue, si S V?(é) < 0, il existe une valeur propre A_, définie par 
A_; = maximum R(U, V) (30) 
(U, V) 
quand on considére l’ensemble des paires de fonctions (U, V) telles que S V(é) < 0. 


En vertu de (8) il est en effet évident que (29) ou (30) prennent leur valeur extremum A, ; pour 
la paire de fonctions (U, V) = (y+.1,A41 y41)- 
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9. Lemme III. Nous supposons explicitement qu’il y a un éparti 
et éventuellement c(x) = 0. a ee Be 


Soient V,(é),..., V,(€), r fonctions itérées telles que 


SVG V6) V()=0,° pq, p.qsr. 


Il existe des fonctions admissibles V telles que § V2(é) + 0 et S V(&) Wee a0 p= 1, 6... 
Démonstration. Construisons la fonction itérée 


V(E) = Vilé) + 3 0p Vil). 


Ona 
S V(E) V,(E) = S Vo(é) VE) + op SVRE), = p= ,.--s 7. 
Exprimons que V est orthogonal 4 V,: S V(é) V(E) = 9, soit 


___ SHE) Vp) 


Pp S V3 (é) (31) 


Il reste 4 démontrer qu’il existe une fonction V,(&) telle que S V2(é) + 0. Ona 


S V(é) = S VE) + 2 Y op [— 0, S VAE] + XY 02 S V2E) 
=a S-Ve(c) 05.9 Vee) 


— svg — » SME) VMOP 
Be a SFOs ces 


Nous démontrons cette proposition par contradiction: Supposons que S V2(é) = 0 quelle que 
soit la fonction itérée V,(&). On a 


are) — yr IS Vols) Ve()P 


L’expression (32) étant vraie quelle que soit la fonction itérée Vy, on a en particulier, avec la 
fonction itérée V, + V’, V’ étant une fonction itérée arbitraire, 


S(Vo + VP =SV2EL2SV,V 4S V2 
r l : 
=. x S V3(E) [S (Vo a V’) VF 


l i 
aie VBE) [(SVoV,P +2SVV,SV'V, + (SV Vp) 


et comme (32) est également applicable a la fonction V’ au lieu de Vo, il reste 


S ve) Ve) = 5 EET Me, 


soit 


SV’ (Vet S0,V,)=0 ou SV V=0, (33) 


Op étant défini par (31). (33) est vrai quel que soit V’; la démonstration du lemme II [n° 8. c) 1.] 
peut étre appliquée a (33), elle implique que 


peut (ev) = pw (v% + 2) Op Up) ae [e (v% + D Op *p) |) = Nv) = 0 


quelle que soit la fonction itérée V,. Mais il est manifeste qu’il existe des fonctions V telles que 
N(v) = 0, car il existe un intervalle [x, $] ne contenant pas de masses concentrées dans lequel 
N(v) + 0, il suffit de prolonger la fonction V en dehors de cet intervalle en choisissant v= 0; 
soit V)=— )) 0, V,, les conditions aux appuis sont satisfaites, Vo et V sont des fonctions itérées. 
La contradiction démontre le lemme ITI. 


10. Definition. ([4] p. 259). 
Si A est une valeur propre et Y une fonction propre qui lui est associée, nous définissons la 
fonction Sgn par ; SU ae 
Sees SH X) Ye — OCR) 2 34, 
Send = Sy = Smia)y@—S aera — Nemes YE). ©) 
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Plusieurs fonctions propres peuvent correspondre a cette valeur propre et il convient de noter 
que Sgn A est indépendant de la fonction propre choisie qui est associée a A et 


Send = dais signed. (34’) 


En effet, nous savons que A et S Y? ont le méme signe (n° 8). 

Le cas des valeurs propres multiples peut se présenter lorsque l’arbre a des appuis d’encastrement 
intermédiaires de sorte que ni les moments fléchissants, ni les efforts tranchants ne se transmettent 
au travers de ces appuis. Un tel arbre est alors formé de parties élastiques indépendantes les unes 
des autres. Si on exclut ce cas (ou sil’on traite chaque partie de l’arbre indépendamment des autres) 
Biezeno [5] a montré que les valeurs propres sont simples (il néglige l’effet raidissant des couples 
gyrostatiques Cy’, mais le résultat est évidemment applicable au cas général traité ici). Nous 
supposerons par la suite que les valeurs propres sont simples. 


11. Théoréme II. Nous supposons que la masse répartie n’est pas nulle: (x) 4= 0, ou que c(x) = 0, 
(ou que simultanément yu == 0 et c == 0). 
1) Il existe une infinité de valeurs propres, elles peuvent étre ordonnées en une suite 


cin Alig = aa hy Ag ee (35) 
avec 


Sgn A, = signe p . (36) 


Cette suite est limitée pour les valeurs propres positives et A_,—> —0o pour p—> 00. 
2) A la suite (35) est associé un systéme complet de fonctions propres orthonormales 


eee? Y_2(€), Vege) Yi), Y,(é), Se (37) 
les conditions d’orthogonalité et de norme étant 
Dy ka = 05, Sen A, (38) 
3) Ona 
Ap = minimum R(U, V) (39) 
(U,V) 


pour l’ensemble des paires de fonctions admissibles U, V lorsqu’il existe une fonction itérée V, 
telle que 


S20) 95 1G bese AS De eevee Ls (40) 
et que les fonctions V satisfont 4 ces conditions; - 
A_p, = maximum R(U, V) (39’) 
(U,V) 


pour les fonctions itérées V telles que 
SMI -S VY RO: k= —1, —2,...,—(p—1), (40’) 
R(U, V) étant défini par (20). 


4) Il n’y a qu’un nombre fini de valeurs propres > 0.1 Ce n’est que lorsqu’il n’existe plus de 
fonction V, satisfaisaut les conditions (40) que la suite des valeurs propres positives s’interrompt. 

Démonstration? I. En vertu du lemme II [n°8. a)] il existe une fonction itérée V telle que 
S V2: 0. Supposons p. ex. que § V2(é) > 0; en vertu du théoréme I [n° 8.b)], il existe une valeur 
propre A, a laquelle est associée la fonction propre Y,(é). Nous supposons que /, est une valeur 
propre simple, et A, est la valeur propre positive minimum. 


II. Le lemme III a montré que si V,,..., V, étaient des fonctions itérées telles que 


S V2 0 (si V, est une fonction propre, cette condition 
est automatiquement satisfaite) 


SV,V;=0, p+q, p.qSr 


1 ss . Piet So wet. he 
is Ai den Dungen a attire notre attention sur cette propriété caractéristique du probléme qui nous occupe. 
ans le cas d’un arbre continu sans disque, elle se démontre aisément. La démonstration de Grammel ({1], 


Pp. 200—202) ne peut etre transposée au ¢€as ou T Tr de 1S 
le oto osséde Ss d es, aussi nous emontrer ons directe- 
Pp qu 2 0) d 


2 FuBnote 1 s. S. 374. 
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il existe des fonctions itérées V telles que 
SVv2s- 05 Se 0% a nS (41) 


Supposons le théoréme vrai pour /,,...,A,, il existe donc une fonction itérée V telle que (41) 


soit satisfait, V est orthogonal aux r premiéres fonctions propres. Supposons § V2> 0, nous 
montrons que 


Argan SR(U,V) et Ay = minimum R(U, V) 
(U,V) . 


et la fonction propre Y,,, est orthogonale aux fonctions propres’-Y,,),, 2304s 
R(U, V) est défini par (20), soit 
SS K(&, 0) U(é) U(c) / 
R(U, = eee: 
(Oh) 5 VE) (20’) 


La paire de fonctions (U, V) étant choisie, R est un nombre défini par (20’). Il est évident 
que R = /, et nous montrerons en fin de démonstration que R > /,. Nous considérons a nouveau 
la fonction N définie par 


U(e) — RU, V) S K€, «) U(c) = N¢é). | (42) 


Si N(é) = 0 le théoréme est démontré, 1,< R=4/,.,. Supposons done N(€é) = 0 et considérons 
le probléme auxiliaire 


U(é) —AS K(é,c) Uo) = N. (43) 


Ce probléme admet la solution U(é;A) tant que A n’est pas valeur propre. 
Remarquons d’abord que la fonction U(é) dont est issue la fonction itérée V est orthogonale a 


Dee > ¥- Of a 


A V(é) = $K(é,0) Uo) 
et 


S$ U(é) Y,(6) = Ap SS K(E. 0) Y,(c) U() =4,S V(o) Y,(0) =9, pa=h...or. 


Remarquons ensuite que la solution U(é;A) du probléme auxiliaire (43) est également ortho- 
gonale aux r fonctions propres Y,,..., Y,. En effet on a, avec (43) et na) =A, S K(é,c) Y,(0), 


done S U Y, =1, SS K(é,0) Y,(0) U¢), 
AS S K(é,0) U(o) Y,(é) = S$ UE) Y,(é) —S NE) Y,(6) | 
=/, SS K(,9) ¥,(0) U) —$ Y,(€) [UE — RS KE, ) U(o)] 


= hy SS K(E, 0) ¥,(0) H(€)—S Y,(€) UE) + RS SK(E, «) Ulo) Y,6) 
=A, $$ K(é,0) Y,(c) U(é) +0+0 
d’ot 2 
or A >A, donc U est orthogonal a Y,P=1,...51. = 
A la solution U(é;A) de (43) associons la fonction itéreé V(€; A) et notons que pour A = R(U, V), 
ona U= U, v= Vy. 
Considérons encore la fonction auxiliaire H définie par 
H(A) = S N(é) O(E;a) = S TEs) V(Es2) —ASS KEE, 0) Ulos4) (Es A). 
Pour A= Rh, (Rh) = 0. 
Pour A = /,, “ n 
H(A,) = S O(Es4,) V(E54,) —4, SS K (&, 0) U(osA,) V(Es Ar) - 
Calculons le nombre R pour la paire T, Vet pour J = /,, soit 
| S H(é;A,) Ves A,) = R(2,) SS K(E,0) U(os4,) ViEs4,) « 
Ona acted. ba 4 
| H(A.) = [RU,) —A] $$ K(E, 0) Ho3%,) Hes A) 
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Or R(A,) =A, donc, en vertu du lemme I, H(A,) = 0. On démontre, comme pour le théoréme I, 
que H’(A)=S V(E;A) et que H(A) < 0 et diminue de fagon monotone a gauche de R(U) tant que 
A=: valeur propre et pour 1 >/,. Or H(A,) = 0, il existe donc une valeur propre /,;1 telle que 

yey res #1 CU (44) 


Il reste 4 montrer qu’on a bien/, << R(U). 
Puisque U et U sont orthogonales 4 Y,,..., Y, et que U-+ ¢ U est admissible quel que soit 


le nombre réel ¢, on a R(U + ¢ U,V +¢V) =A, pour tout |e| suffisamment petit (dans ce cas en 
effet R reste > 0; autrement il pourrait changer de signe), soit explicitement 


SS K(é,0) [U(o) +e 0()] (US +2 TO]—4, SIV +e VOP = 0 
ou, en développant 
SS K(é,) U(o) Ué) + 2e S$ K(é, 0) U(o) H(é) + 2 SS K(E, 0) He) UE) 
—A,[S VE) + 2S VE) WE) + & S VG)] SO. 


Nous démontrons maintenant (44) par ’absurde. Supposons que /, = R(U), les termes soulignés 
disparaissent et on a 


2 e[SS K(é,0) U(o) (O —A,$ VO VO] = —#[SS KE, o) Ho) HE) —4, S VE]. 


Les termes entre crochets sont des nombres fixes, comme le| est arbitrairement petit, et que le 
signe de ¢ est arbitraire, on obtient la contradiction qui démontre (44). Le cas des valeurs propres 
négatives se traite d’une fagon analogue. 

En vertu du lemme III, il existe une infinité de valeurs propres lorsqu’il y a des masses réparties. 
Il reste 4 montrer qu’il y a un nombre fini de valeurs propres > 0. Nous reportons cette question 
plus loin (n° 14). 

12. Autre mise en équation, équation différentielle. Nous reportant aux relations (5) et (7), 
considérons une tranche d’épaisseur 4x de l’arbre (figure 4), qui se comporte comme un disque. 
Si T est Veffort tranchant, le théoréme de la résultante cinétique s’écrit 


0= T,— T_—yp Ax (—o?’y), soit a Sars 


dx 
Zz rf. 
Oje=ty 
om 18 Az 


Fig. 4. 


et le théoréme du moment cinétique s’écrit, M étant le moment fléchissant 


ote = ¢Axaty 1,= —(M, —M_—T Ax], 


soit 
= =—coy+T. 
On al ; 
(E a Sai Ss M”’ 
= wey —T 


= w(cy’)’ + wy 


(Ely) Alp yt (ey): (45) 
C’est l’équation différentielle cherchée. 


; Eee eee 3 ; : 
Si nous considérons une tranche d’épaisseur Ax, localisée en x = s ov se trouve calé un disque, 
ona 


soit, en posant J = w? 


(ELy\,0—(ELy”),_o =Amy(s), } (46) 


(ETy”).40—(ETy”),0 =A Cy'(s). 


; ade ae a ae f 2 eae ; 
a Rappelons encore ici que nous nous limitons a utiliser les équations simplifiées de la résistance des ma- 
tériaux. Nous négligeons donc la déformation de l’arbre résultant des contraintes tangentielles. 
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Ce . . . . 
: sont les conditions de passage au droit des disques. I] faut encore adjoindre les conditions 
aux limites. Si p. ex. l’arbre est simplement appuyé: 


¥(0) = ¥"(0) =x) = yO =O. (47) 
systeme (45), (46), (47) est équivalent au systéme (8), ou a l’équation intégrale (9). 
es théorémes I et II s’écrivent, si U est une fonction admissible et V la fonction itérée telle 


que S V2(é) > 0 
(<) 


SUV _ S(muv—cuw'v’) 
Si S(m v2 — ¢ v’?) i= Riu, v)) : (48) 


Ona 
V =S K¢(é,c) U(o), 


v(x) = S a(x, s) m(s) u(s) — S B(x, s) e(s) u’(s) , 
v'(x) = S y(x, s) m(s) u(s) — S 6(x, s) e(s), u’(s) 


ya , > . 
et v est la déformée de l’arbre soumis aux ,,forces“ m u(s) et aux ,,couples‘‘ — c u’(s). Calculons 


soit 


(49) 


1 
779, . , . , . . 
f E Iv’? dx en transformant cette intégrale par des intégrations par parties, en tenant compte 


des conditions aux limites et des conditions de passage: 
l l 
Loe Tv’ vw dx = J [EI vw’). — SY (ET o’y of + f(El0”) vdx. 
0 
Les termes intégrés sont les sommes 
s,—0 s,— 0 s.+0 
ets | es 


ou la somme Y’ est étendue aux disques. On a donc 
i 


t 
i EIv'2dk= [»(s) (EI ane Be [»’(s) (EI a) has ee i (EI v")’vdx. (50) 
0 


On dira que la fonction v est la fonction itérée d’une fonction admissible u du systéme (45), (46), 
(47) lorsque 
(E I vy" eS. Lb u - (c u’)’ 
[(E rey |, = m uls) » [Ere [iP = ews). (51) 


0 


(50) s’écrit 
JE Lv"? de = 5 m(s) u(s) 26s) — 5 ols) w(s) 0%) + flaw (ewy] ode 


= idem idem ++ Sti uv—cu'v') dx 
=S(muv—cu'r) 
et (48) s’écrit 
fELv'2 dx 


0 ae Dv) Bess , 
ee ed R(v) (48’) 


<a 
, (>) 
od les fonctionnelles ® et Y sont définies par (48’) et R(v) = R(u, v) > 0. 
(<) 

Les théorémes I et II ont été démontrés directement 4 partir de l’équation différentielle auto- 
adjointe, de ses conditions aux limites et conditions de passage par Kestens [6] (Le systéme (45), 
(46), (47) n’est qu’un cas particulier d’une théorie générale). Il faut remarquer que les démon- 
strations de ces théorémes données ici sont beaucoup plus simples et plus courtes, ce qui en fait 
Vintérét. Par contre les démonstrations de la réalité des valeurs propres et des propriétés d’ortho- 
gonalité des fonctions propres sont plus rapides par la méthode qui part de l’équation différentielle 
avec ses conditions aux limites et de passage [6]. L’interprétation physique de la fonction itérée v 
issue d’une fonction admissible u, suivant le systéme (51) est évidemment identique 4 celle donnée 


par (49). 
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13. Interprétation physique des valeurs propres négatives. Soit A_p =—A(A> 0) une valeur 
propre négative. Le systéme (8) s’écrit 


y=+ASa(—my)+ASfey’, 
y =ASy(—my)+ASdcy’. 


Il traduit l’équilibre d’une poutre soumise a des couples +- Acy' et a des réactions élastiques 
du sol A my sur lequel elle serait posée (figure 5). 


Zz 
Ac y' Ac y’ 
aa i 
y 
Am y Amy 
Fig. 5. 


14, Fin de la démonstration du théoréme II: I] y a un nombre fini de valeurs propres 
positives. 

Cette démonstration nécessite l’emploi de deux résultats complémentaires que nous rappelons, 
sans en reproduire les démonstrations. 

1) Théoréme de comparaison ([6]): Les conditions aux limites étant les mémes, supposons 
deux problémes ayant méme P(v) et tels que 


P*(v) < Pv), 


P(v) et Y(v) étant définis par la formule (48’). (Il s’agit done de deux arbres qui ne different que 
par la grandeur des masses ju(s) et m/(s;) et des différences de moments d’inertie c(s) et C(s;), mais 
qui ont méme E I(x) et qui sont géométriquement les mémes). On a 


i eee 


Ce résultat s’obtient immédiatement de (48’) pour la premiére valeur propre. Pour les autres, 
il faut faire appel au théoréme du minimum-maximum de Courant [6]. ; 

2) Un arbre sans volant et tel que (x) == 0, c(x) == 0 posséde un nombre fini de valeurs propres 
positives (vitesses critiques) ({1] t. 2, p. 200). 


Démonstration: Reprenons le systéme (45), (46). 


I. Considérons wm = y A (A> 0) comme une fonction du paramétre jw, soit plus explicitement, 
en posant u(x) = k u(x), fonction du paramétre k que nous faisons varier de 0 4 + co. Le nombre 
de vitesses critiques ~, variera d’une unité chaque fois qu’il existe une asymptote verticale aux 
»courbes propres“ ~,(k), c’est-a-dire s’il existe des k pour lesquels @ = 00 ({1] p. 200)!. Or il est 
aisé de voir que dans le cas général il n’existe pas d’asymptote. En effet, s’il y en avait, il fau- 
drait, en exprimant que 1 =©o, 


py tey) =0. ou. — (cy) =k wy. (52) 


2 
conditions aux limites: 


par exemple, conditions aux disques: as 

: 

¥(si) = y'(s;) = 0. (52") 

k apparait comme valeur propre d’un probléme auxiliaire. I] est manifeste qu'il n’existe pas de 

solution non triviale au systéme (52) qui satisfasse aux conditions supplémentaires (52). Il n’y 
a donc pas d’asymptote et le nombre des vitesses critiques est constant quel que soit k. 

Remarquons cependant que les conditions (52’) peuvent étre satisfaites individuellement par 

une solution non triviale et si par exemple tous les C(s;) sont nuls (les volants se réduisant a des 


masses localisées), il ne reste comme condition supplémentaire que y(s;) = 0, il peut y avoir des 
solutions, et le nombre des vitesses critiques variera dans cette éventualité. 


1 Ceci suppose les valeurs propres simples. 
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Ecartons pour l’instant ce cas: 

II. Supposons m(s;) C(s;) 4: 0. Puisque dans ce cas le nombre de vitesses critiques reste cone 
stant quel que soit k (c’est-a-dire le nombre de valeurs propres positives) appliquons le théoréme 
de comparaison rappelé ci-dessus d’une part au cas k = 0 (le systéme étudié ayant le méme nombre 


de valeurs propres 2; > 0 que celui pour lequel k = 0) et d’autre part a l’arbre ne portant que 
les volants, et supposé tel-que 4~—=c=0. Ona 


l 
Py) = — fey’de + ¥ mls) Xs) — ¥ Cls) ¥%(5) SY mls) y%s) — Y Cla) ys) = W" 
et donc 
Nees 

On en déduit que le probléme dont est issu Y a un nombre fini de valeurs propres positives 
comme Villustre le schema figure 6. Certaines valeurs propres 1* du probléme * de comparaison 
peuvent étre rejetées 4 + 0o pour A, comme certaines valeurs propres 1* < 0 peuvent devenir 
des valeurs propres 1 > 0 (pointillés de la figure 6). Puisqu’il n’y a pas de point d’accumulation 
des valeurs propres a distance finie (nous avons supposé les valeurs propres simples, mais le théoréme 
II est également vrai dans le cas des valeurs propres multiples [6]), il ne peut y avoir qu’un nombre 


fini de valeurs propres positives et le nombre de valeurs propres positives du probléme considéré 
(k == 0) est le méme. 


‘Cai 


7}, 


ACA) 
(Ay) (As) 


yeas Pe 
(Aj) Az) (As) 
Fig. 6 


—) 


Il reste 4 examiner les cas m(s;) C(s;) = 0 quel que soit 1. Ces cas III, IV et V ci-dessous échap- 
pent a la démonstration IT ci-dessus. 


| 


III. Supposons m(s;) = 0, C(s;) == 01. Comparons ce probléme a un probléme ot m(s;) + 0. 
Ona 


W = | (uy*—cy) de— ¥ Cbs) y's) 5 J (uyt—ey’) de + 3 ms) sD Cla) ¥%s) = ¥* 
soit 
Aj zat. 
Or en vertu de II ci-dessus, il y a un nombre fini de valeurs propres Ai > 0 et la proposition 
est donc également démontrée dans le cas ot m(s;) = 0. 
IV. Supposons m(s;) C(s;) = 0 avee m(s;) = 0, C(s;) + 0 
m(s;) += 0, C(s;) = 9 


Introduisons des masses m*(s;) - 0 dans le systéme de comparaison, on a 


} bk 


V = | (yey) de +E mls) y%6) — B Clo) y6) Sf (ry*— ey) de 
ae m/(s;) y?(s;) ++ J m*(s;) v*(s;) — 2 C(s;) y'7(si) = ide 
A; = AF 
et en vertu de II, le systéme de comparaison n’a qu’un nombre fini de valeurs propres positives ; 
= aay ee ee Se ne ne peut conclure comme ci-dessus en appli- 
a dL ne a eo Petals comparaison tel que Y* > Y et pour 


lequel il serait possible de tirer une conclusion est ici le cas ou c(x) = 0, et il y a une infinité de 
valeurs propres positives pour ce probléme. 


soit 


1 Ce cas est purement théorique et n'a pas de réalité physique. 
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Par contre, en appliquant le théoréme de comparaison au systéme pour lequel m(s;) = 0, on a 
1 l 
Wet — f(uy?—ey’) dx Sf (uy?—ey) de + 3 m(s) (si) = ¥ 
0 0 


soit 
Ap eeahy 

Il y a un nombre fini de valeurs propres 1#* > 0 mais le nombre de valeurs propres 2; > 0 
n’est pas le méme que le nombre de valeurs propres y heomiss (18 

Reprenant le raisonnement du début de la démonstration (en I.), il peut exister des valeurs 
propres k du probléme auxiliaire (52) (52’), mais comme les abscisses des asymptotes verticales 
aux courbes propres @,(k) n’ont pas de point d’accumulation a distance finie, les nombres des 
valeurs propres positives A; et A** ne peuvent différer que d’un nombre fini. 

La proposition est donc démontrée dans tous les cas. 


15. Conclusions. Les démonstrations que nous avons été amené a développer au sujet du 
probléme particulier des vitesses critiques, s’adaptent sans plus a une large classe de problémes 
aux valeurs propres [6]. C’est plutét comme illustration d’une théorie générale 4 un cas simple, 
qu’il faut considérer les développements précédents; ce n’est que le n° 14 qui est particulier au 
probléme traité. 

Au point de vue physique, l’introduction de disques, comme d’ailleurs de fagon générale l’intro- 
duction de masses concentrées dans un systéme continu quelconque, ne modifie pas l’aspect quali- 
tatif; aussi il n’est pas étonnant que nous retrouvions les propriétés des vitesses critiques d’arbres 
continus [1], [2]. 

Au point de vue mathématique par contre, une nette séparation est généralement faite entre 
les systémes continus et les systémes discrets; les premiers sont souvent considérés comme un cas 
limite des seconds. Du point de vue de l’esthétique des démonstrations, il est cependant satisfaisant 
de traiter simultanément les deux cas. 

A cette fin une premiére méthode consiste a traiter l’équation différentielle 4 laquelle sont 
jointes les conditions aux limites, et les conditions de discontinuité de certaines dérivées aux endroits 
de masses concentrées. C’est particuliérement dans cette voie qu’on constate que les démonstrations 
ne sont pas une simple généralisation du cas continu [6]. Une autre méthode consiste 4 traiter les 
problémes aux valeurs propres sous forme d’équation intégrale [3]. Une large classe de problémes 
peut ainsi étre traitée, en introduisant un noyau réel symétrique, grace 4 des changements de 
variables et de fonctions analogues 4 ceux du n° 4. 

La théorie des équations intégrales 4 noyaux symétriques est immédiatement applicable. La 


théorie est également développée pour le cas de noyaux complexes hermitiques: K(s, t) = K(t, s). 
Le probléme que nous avons traité est un peu différent puisque nous avons introduit un noyau 
complexe symétrique. I] est remarquable que les résultats sont également valables dans ce cas. 

La méthode des itérations et la détermination de bornes supérieures et inférieures peuvent étre 
aisément transposées. Pour ces questions, nous renvoyons le lecteur 4 [6]. 
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General Stress-Strain Relations in Non-linear Theory 
of Viscoelasticity for Large Deformations 


By T. Tokuoka 


1. Introduction. The most general phenomenological one-dimensional linear theory in visco- 
elasticity has been established mathematically by B. Gross}. Here a Maxwell- or Voigt model 
constituting series or parallel connections of a spring and a dashpot are used. Introducing the 
relaxation time spectrum and the retardation time spectrum on the basis of the superposition 
principle, the properties of one-dimensional linear viscoelastic substance are characterized. 

Recently, as the result of studying the mechanical behavior of highpolymer for large defor- 
mations, the Weissenberg-effect and Poynting-effect i.e. so-called normal stress effects have been 
observed; hence it is demanded that the theory of viscoelasticity be extended to three-dimensio- 
nality and non-linearity. Here it is required to establish the general three-dimensional non-linear 
theory of viscoelasticity. 

Several papers? have been presented to extend the viscoelastic theory to three-dimensiona- 
lity. M. A. Biot and H. Ziegler study the irreversible process thermodynamically; W. Olszak and 
P. Perzyna study the three-dimensional theory by means of operational method, but the above 
methods are based on linearity i.e. the deviations from the equilibrium state are sufficiently 
small. Furthermore several attempts’—!° to generalize the theory to non-linearity are presented, 
but these are formal generalizations of the fundamental relations of one-dimensional Maxwell- and 
Voigt models into three-dimensions, so it is not satisfactory from the physical point of view. 

In this paper, we generalize Novozhilov’s non-linear theory of clasticity for large deformations" 
into viscoelasticity, and we obtain the general stress-strain relations of the non-linear theory of 
viscoelasticity for large deformations. 


2. Equilibrium Equations. With respect to the orthogonal coordinate system X‘ (i = 1, 2, 3) 
fixed in space, a point x' of the unstrained substance for t = 0 transforms to a point x** for time t 
as a result of the deformation of the substance. 

Then we have 

Paik ahd (rane: Coe oa} () 
Awe (a, a, x; 2) 
and 
wat, 75.0°;-0). = 0. 
Generalized strain components are defined by 


1 


ej=— Bj—%)) (J =1,2,3) (2) 
where 

Bij = ait xX (3) 
and 

eh Ox* k 

st Oxi 


1 B. Gross, Mathematical Structure of the Theory of Viscoelasticity, Paris 1953. 
2M. A. Biot, J. Appl. Phys. 25 (1954) p. 1385. 
3 M. A. Biot, Phys. Rev. 27 (1955) p. 1463. 
4M. A. Biot, J. Appl. Phys. 27 (1956) p. 240. 
H. Ziegler, Ing.-Arch. 25 (1957) p. 58. 
W. Olszak and P. Perzyna, Ing.-Arch. 28 (1959) p. 246. 3 
T. Alfrey, Mechanical Behavior of High Polymers, Interscience Publ. New York 1948. 
W. Segawa, J. Phys. Soc. Japan 12 (1957) p. 996. 
9 M. Yamamoto, J. Phys. ge J aver Gaseous! | 
10 W. Segawa, J. Phys. Soc. Japan p- . a 
11 a scien ae pe ede of the Nonlinear Theory of Elasticity, Rochester, N. Y. 1953. 
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6 
7 
8 
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From its definition the generalized strain tensor is symmetric: 


Eij = Fji- (4) 


In these and following expressions we adopt the so-called Einstein summation convention to 
simplify mathematical expressions. 

Now we consider the equilibrium of an Insert elementary rectangular parallelepiped in the defor- 
med substance, which has the side dx*‘ parallel to the coordinate axis X‘. Moreover o*' is the stress 
on the plane perpendicular to the X;-axis, F* is the body force per unit volume and we neglect 
the inertia force. We can next obtain the generalized equilibrium equations in the same way as 
in the classical linear theory of elasticity: 


bo*t 
Ox*i 


ole FM O's (5) 


Equating to zero the resultant of the moments of all forces acting on the parallelepiped 


otii = oti (6) 


where o*'/ is the component of o*' projecting on the X/-axis. 

The equilibrium equations (5) and (6) do not differ in form from the conditions for equilibrium 
of a parallelepiped derived in the classical linear theory of elasticity. They differ in essence, how- 
ever, because in the linear theory of elasticity, the change in position of the points of the body 
due to its deformation are neglected. 

Equation (5) is expressed in Kuler’s representation, so we will transform the equilibrium equation 
(5) into Lagrange’s representation. 

Now from (1) we obtain 


dx*i — aX dxi 


(7) 


or ‘ ; oe A] 
— = *L 
dxJ = x), ; dx =z & ’ 


where 


. ax) . 
HS yi => Oxti ? D = det (x) 
is the volume dilatation factor and A/ is a cofactor with respect to the aS element of D. 
According to the mathematical relations 04//dx/ = 0, the equilibrium equation (5) is transformed 
into 


a( Ad o*') 
MEA 
a DF*=0; (8) 
here, the independent variables are x1, x?, x3 in Lagrange’s representation. 
Let us denote the fixed unit vector along the X‘-axis by e; and the unit vector in the direction 
of that axis after deformation by é;; so from (7) 


Peet 
4; Tan e; 
vu 
or at (no summation on i) (9) 
cos (7;,€;) = Va 
ti 


where (i;,.€;) denotes the angle between the vectors 4; and e;. 


Let us denote the unit vector perpendicular to two directions ¢, and i; by m, [here, (i, j, k) 
is a cyclic permutation of (1, 2, 3)]. From the relations 


and 


Vsii 8ij — (8: 5 


wiSjj 


(no summation on i and j) 
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it follows that 
cos (14. e;) = Ny, . e; 


1 Sie 
== e,: (i, X 4) 


sin (a;; ) 
*m kN 
Xi xij 
— i _ 7 e| : (2, x e€,,) 
V8es233 — (ij) 
Aj se 
a (no summation on i and j) 
Vsss837— (ay 
PS eae tionton:k 
Sed no summation on k) 
where 
S*k De oR OO ees a EE 
—— V gi: gj; sin (i;, t)) = Vg: 8j7 — (8:;)" (no summation on i, j and k) 


is the ratio of the area perpendicular to the X*-axis before deformation to its area after defor- 


mation, 4 
Then the stress o* on the plane which have a perpendicular n, is 


S, = o*! cos (n;, €;) 


Ik 
ES Aj *1 
V8: 875 — (83 
ae : 
= GEE Aj o* (no summation on k) . 


Therefore the equilibrium equation (8) is reduced to 


*j n. 
as si)+ DRY =0, 


axt \ Si 


We oP Ae . ° ° . n : 
Resolve the vector o/ in the direction 7, and denote its component by o/!: 


*7 

n, BN ail X,1 n, 

ai = oF! 4, = Ose, 
81l 


* J RP on 
gf 5/2. ued g-\ cD B® <.0:, 
dx? \ SI Ven 


Projecting this equation on the X'-axis, we have 


© (ati ctl) + D F*i= 0 


axJ 
where 
L nm 
Sa a a Mer : : ; 
tes. (no summation on i and j). 
Ss Ves 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


Now the magnitude 1*‘/ can be represented in terms of the stress o*'/. Multiplying the both sides 


of (12) by a vector e; and a cofactor A} and then summing with respect to i, we obtain 


n V A Line ah 
git — sui gi.e; (no summation on /) 
D 


and substituting (10) into the above equation, 


a j i 4h 
Oe ee SEES Gt @, 


St 
Mee pee 
= ae : Ei An A Cae (no summation on j and I). 
S 


27* 
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Therefore substituting this relation into (15), +*'J is reduced to the following expression: 
; J » 
t* bjies Ban. gtim : (16) 


so that 
whi rtit = Al tim, (17) 
These are the relations which relate the magnitudes t*'/ to the stresses o*'/. 
The magnitudes t*‘J defined by (16) are generalized stesses and are symmetric with respect 
to.suffixes i and j from (6). 
Writing (14) in more compact form, we define 
Tess (hia le, (Sah 22.3) 5 (18) 
then we have 
div T'+ DF*'=0. (19) 


This is the equilibrium equation, expressed in Lagrange’s representation, for the case of the most 
general deformation of the substance. 


3. Hamilton’s Principle. Now we consider an arbitrary region V (before deformation) in the 
substance and its surface S; the surface force is given on the portion S; of S. The substance which 
is in unstrained state at time t = 0, deforms viscoelastically to time t. A(t) denotes the cumulative 
elastic energy in the region V from the unstrained state (f = 0) to the strained state (time 1); 
R,(t) and R,(t) are the work done on the substance by the body and surface forces respectively; 
B(t) is the rate of dissipation of energy in the region V at the time ¢; @ is the elastic potential i. e. 
the elastic (recoverable) energy per unit volume before deformation; and Y is the dispation func- 
tion, i. e. the rate of dissipative energy per unit volume before deformation. 

Then 


ies Nn es pal US (20) 


From the law of energy conservation, we have 


A(t) + i B(t!) dt’ —[R,() + R,()] = 0, (21) 


where we suppose that the deformation is adiabatic. 

Now we assign to the displacements u‘ arbitrary virtual increments du‘ (du' = 0 at time 0, 2) 
and assume that the state of the body satisfies the following variational equation from Hamilton’s 
principle: 


F140) = j Be’) a!’ — [R,(e’) + Rwy} d= 0, 


t t’ 
uf {0 + J OB(t'') dt!’ — [6 R(t’) + SRy()]| de =); (22) 
We may regard this equation as generalized Hamilton’s principle. 
From (20) 
6A = fff dOdV, 6B = fff d¥dvV, (23) 
OR, = fff (F** 6u’) DdV, 
OR = [J dul) dO* = [ff ou’) Se dO | 4) 
f | f 


where f** is a component along the X’ -axis of the surface force acting on a unit area of the surface 
of the deformed body, d(2* is the differential of this area, and the coefficient S *n/S" is the ratio 
of the elements of area in the terminal and initial states. 

Furthermore we adopt the important assumption that the elastic potential depends only upon 
the strain and not upon rate of strain. However, the dissipation function is a function of both 
strain and rate of strain: 

@D => D(E; ;) 5 (25) 
= Y (8:5 ij) . (26) 


~~ 
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Then, these variations are 


oD 
6@ = = 0; j (27) 
and 
OE he ov 
Ov — 68: 06; ; | 88 j 0&5; (28) 


From (2), (3) and (1), the strain ¢; ; 18 represented by the displacements u’: 


I 0 ; : 
A 5) (wi; | ur; | ui ur) 


its variation is 
1 , p 

0&5 = > [(du),; + (du'),; + (du*),; uk; + uk (du*), |] 

a. 1 

ane 


[x**(du*) “4% *(6u") 5] : 

Substituting this into (27), we obtain the variation of elastic potential O: 
Lee 

6D = xi Ben j (du i 


- = P’. grad (6u') 


= div (P' du‘) — (div P*) du’ , (29) 
where 
, 0 
P' — (wi oe ej (i = is Ze 3) . (30) 
So the variation of cumulated elastic energy in the region V is 
6A = — fff (div P’) du' dV + ff (P'), dui 62, (31) 
v Ss 
f 


because the volume integral of the divergence of a vecter can be reduced by Gauss’ theorem to a 
surface integral; where (P'),, is the projection of the vector P' on the exterior normal to the bound- 
ing surface. 

Namely the variation of rate of strain is 


Bec = 8 (it, + thy + wh, why + uh, why) 
= > [xt H(0i"),; + (S04), 
+> LM (Out),5 + et MOu"),c)- 
Substituting this into (28), we obtain the variation of dissipation function Ue 
OY = ah! Fo (Oi) + (2 Ga + AM Ga) Oda 
= U'- grad (6u') + V‘- grad (du’) 


= — (div U‘) du‘ — (div V‘) dui 
+ div (Ui du’) + div (Vi du’) , (32) 


where 


(i = 1,2, 3). (33) 


Therefore the variation of the rate of dissipation of energy in the region V is 


5B = — fff (div U) die dV — fff (div V*) du! dV 


+ ff (U4, di dQ + ff (V'), dui dQ. 
re Sf 
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Since the variation of displacement du‘ vanishes at the time 0, 


foB(e') de” = i= J ff (div U) dui dV + [f (On du! an 


t= 0 


f 
+f {fs [div (U!— V')] dui dv — fi (ui—V'), dui iQ) at! 
=|— If! (div U') dui dV-+ ff (U7), dui dQ} 
a i Sf J (div W') du! dv — ! (W'),, dui dQ| a (34) 
where 
ld rat bet Aone fa e; . (=1,2,3). (35) 


Substituting (31), (34) and (24) into (22), we obtain 


Re i 
(Pt OS an ae 


fi-ss [div (P! + U') + DF*'] sui dV + ff 
0 Sa 


+ fae f [JJ (div W%) bu‘ dV — [f (W'), du’ an dt’ =0. (36) 
0 0 V f 


je a9] dt! 
t=t’ 


t=t/ 


The above equation should be valid for arbitrary value of the variations of displacement du‘; then 
we obtain 


div (Pi + P)) + DF*i=0, (37) 
(Pi + U), — 5" pri = (38) 

and 
wi=0. (39) 


4, General Stress-Strain Relations. By comparison with the equations (19) and (37), we can 
deduce the stress-strain relations: 


T=P UY. (40) 
In order to express these relations in terms of components, we substitute (18), (30) and (33) into (40), 
snr . | oD Ch 
Kip jK __ ..%*i 
Xk vie! XK ai 
Tees 
. aD ow 
ie 

B dej— 087 ey 


These are the required general stress-strain relations, which show clearly that the generalized 
stress tensor (t*‘/) is symmetric. This equation expresses the so-called generalized Lagrangian form 
of viscoelasticity. 

It can be seen from (41) that the stress t*‘J(t) is expressed as a summation of a term only depend- 
ing on the elastic potential © at time t, and a term only depending on the dissipative function V 
at time t. Now we denote 


Cm: od Us ow 
Pte eS (42) 
then 
: ring +7, (43) 


These relations express no more than a generalized Voigt model, and indicate that the stress-strain 
relations should reduce to a Voigt model in the limiting case of the assumptions (25) and (26). Equa- 
tion (43) is the fundamental assumption that M. Yamamoto! and W. Segawa? adopt formally 
as the stress equation for the three-dimensional Voigt model in their papers. 


1 Cf. footnote 9 on page 385, 
* Cf. footnote 10 on page 385, 
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From (38) and (40), the boundary conditions on the body surface S; are 


A S*n ee 
Se aera =0. (44) 
The viscoelastic condition (39) is reduced to 
CELA ae 
dt FeAl 08; aie 3 (45) 


The physical meaning of the above condition should be investigated thermodynamically. 


5. Isotropic Substance. We consider now the isotropic substance with respect to the rate of 
strain as well as the strain itself. We denote that 


ay = El» 

1 lm 
325 Imn 
UES 31 Onan Ep Emq Enr 
and 

by = & 
ieee 

by = ar Ove Elp Emq 2 (47) 

—— 1 6lmn @ - ‘ 

Siam Sif pqr “lp Emq Enr 


whare a, da, a3 are the strain invariants with respect to an orthogonal transformation of coordi- 
nates, b,, 6,, bs are the rate of strain invariants with respect to the same transformation, and eG 2 
is a generalized Kronecker delta. 


Assuming that 


D = D(a, ag, as) , ne; 
P= Y(b,, bo, bg; ay, Ao» a3) 3 
then 
oP ae aD OQ, ow ty ow Obm (49) 
G6 j 8m 8645” Bj ODm O8: j 
and 
aa, gs 0a, sim dag ted Secs a 
def % > nie Oy Ema? a Ofer Erngnee 
ab, ra) ab, ee aes dbs ie a gimn 3 nie 
06: j A] 06: j JY “mq? 08 j Q “iar “mq “nr 


Substituting (49) and (50) into (41), the stress-strain relations of the isotropic substance can be 
characterized by the six invariants 0®/0a,,, 0Y/0b,, (m = 1, 2,3), which are specified by the 
properties of a given substance. 

We can rewrite the Lagrangian form (41) in matrix representation by putting 


(c*i/) =T, (51) 


i 7 1 imn a 
(05) =I ? (O;¢ Ema) ts A, ? & Ovar Emq én iA As ? 


(52) 
i 5 1 imn ¢ H Ss 
(51 8,4) = By, G 5m” Bag Ene oy ae 
Then the stress-strain relation is 
oD, , @ aD ow, , ow aw = 
Sees as oa, 8 | a the ee tn, BP (53) 


‘Namely, the stress is a linear function of I, A,, A, B, and B,, and this coefficients are the strain 
invariants, or the rate of strain invariants, which are specified by the substance. 
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6. The Case of Small Strains and Small Rates of Strain. When the elongations, shears and angles 
of rotation are small, the second order or more of strain components and angles of rotation can be 
neglected. So from (2) the generalized strain component ¢;; is reduced to 


1 
ij oe 


(ui; + u' ;) = ij» (54) 
where e;, is the ordinary classical strain component. It is easily verified that the rates of the area 
S*i/Si (i = 1, 2, 3) and the volume dilatation factor D = V*/V differ from unity only by magni- 
tudes of the same order as the elongations and shears. Also the matrices (x*") and (Aj) differ 
from the unit matrix only by magnitudes of the same order as the elongations and shears. 


Then from (16) and (17) 
ati gti! wots! watt wot, (55) 


where o'/ is the ordinary classical stress component which acts on the plane before deformation. 
The specific body force per unit volume after deformation F’* is approximately equal to F which 
is the specific body force per unit volume before deformation. i. e. 


DF* ~ F* ~F. (56) 
Substituting (55) and (56) into (14), the equlibrium equation is reduced to 


dod : 

Ct 
ra ia) (57) 
This is identical with the equilibrium: equation of the classical theory of elasticity. 

Let us expand the functions ® and V in a power series in the strains e;; and rates of strain e; } 
respectively. No negative powers can appear in the series, for otherwise the elastic potential and 
the dissipation function would tend to infinity for infinitesimal displacements or infinitesimal rates 
of displacement of the points of the substance from their initial undeformed stationary position. 
Furthermore, if the substance is free of all stresses in the vanished strain or vanished rate of 
strain state, then we must begin the two series with the second power terms. 

Thus we can write 


d ieee i lag es 
eS (oe vkl, 
Urals Cj et gp ej Ct Cmn T ++ +> 


fais =) eee oo) 
oy ij Oat ay ya ae 
so that 
a@ bs ed 
26g) EE eth EE ey Onin += 5 | 
ov ijkl ij,klmng ye 
ae Pe ey POO” Cry mn Toes 


Neglecting the terms of the second or higher power in (59) and substituting (55) and (59) into (41), 
we have 


oti = nie, + iil a, (60). 


which are the linear stress-strain relations of a three-dimensional anisotropic Voigt substance. 
The viscous coefficient y‘j>*' is in general a function of the strains; but when the dissipation func- 
tion is independent of the strain, the coefficient y‘j-*! is a constant which is specified by a given 
substance. 

Furthermore, if the substance is isotropic with respect to strain and rate of strain, the elastic 
potential and the dissipation function are reduced to 


1 1 
D =F M145 + GoM, P= > pr bi + po be, (61) 


where 9, Po Y, and y, are all constants which are specified by the properties of the given sub- 
stance. 
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Then (53) is reduced to 


T= 9, 4,1+ 9, A, + y,6,1+ y, B,, (62) 
namely 
Tey: Saas 
oO oO oO a, 0 C20 + 33 —ey5 — 215 
Pe aaa ee Va 
entre PSs = Coane ay aienen ee S370. Cag C74 5 
Si A8a- 93 
eas Pe 0 “5 oa 31 — 32 C14 1 eae, 
by 0 Coo + 33 — eyo =P 
+ Yr: by + Ye ° —€51. yg C14 €ag_| (63) 
0 by me C3T mea Ci, + C25 
After simple calculations we obtain the following compact stress-strain relations: 
wy = / Lard 
ot = 2G6e¢,;+2ué;;, (64) 
where 


Ln, = dees 0 di, 


Cea. 

1 

Be hg Ou (65) 
26=—gq, 
21 =—Yo> 


and o'/’, e;; are the deviatoric stress component and the deviatoric strain component respectively. 

This is a Voigt model which consists of a purely elastic spring with shear modulus G in parallel 
with a dashpot containing a purely viscous liquid of viscosity w, and the spring obeys Hooke’s 
law for elasticity and the dashpot obeys Newton’s law for viscosity. 


7. Conelusion. General phenomenoligical stress-strain relations in non-linear theory of visco- 
elasticity for large deformations have been presented. 

In the first place, according to V. V. Novozhilov! we express the generalized equilibrium equation 
for large deformations in the Lagrange representation, and we apply the generalized Hamilton’s 
principle to the equation of energy conservation, which denotes that the sum of the elastic energy 
and the dissipative energy is equal to the work done by the body force and the surface on the 
substance; so that we obtain the required general stress-strain relations in comparison with the 
above two equations. 

On the condition that the elastic potential is a function only of the strain, and the dissipation 
function isa function of therate of strain and of strain; such asubstance is reduced to the Voigt material 
necessarily, and the stresses which act on the substance are given by the sum of elastic- and viscous 
stresses, and the stress-strain relations are reduced to the so-called Lagrangian form. 

If elongations, shears and angles of rotation are small and also the strains and rates of strain 
are sufficiently small, the stress-strain relations are expressed by a linear Voigt model constituting 
a Hookian spring in parallel with a Newtonian dashpot. 

Non-linearity in the theory is classified into two groups i.e. the geometrical non-linearity 
and the physical non-linearity. The former is introduced into the theory through the definition of the 
generalized strain and of the generalized stress and through the equilibrium equation for large 
deformation, and the latter through the general stress-strain relations. 

The main result of this paper is that the general stress-strain relations in viscoelasticity are 
deduced necessarily from the physically appropriate assumptions. 


(Eingegangen am 11. Oktober 1960) 
Anschrift des Verfassers: Dipl.-Ing. Tatsuo Tokuoka, Kyoto (Japan) Shijo-horikawa agaru, Nakagyo-ku. 


1 Of footnote 11 on page 385. 
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Grundsitzliche Lésung des Problems der freien Kippung 
rechteckiger Rahmen mit doppelt-symmetrischen Querschnitten 


Von L. Kraus 
Herrn Prof. Dr.-Ing. K. Karas zum 70. Geburtstag 


1. Einleitung. Das Problem der Rahmenkippung wurde von Chwallat, Schoblik? und Eflinger® 
untersucht. Chwalla und Schoblik behandeln einen rechteckformigen Dreigelenkrahmen mit 
Hilfe der elastostatischen Grundgleichungen. Zur Berechnung der Kippkrafte eines solchen Rah- 
mens mit Rechteckquerschnitt (Abb. 1) verwendet Chwalla Potenzreihen’; Schoblik bringt die 
Kipp-Gleichungen in die Form einer Whittakerschen Differentialgleichung?. EBlinger ermittelt 
nach der Methode der kleinsten Quadrate die kritischen Krafte der symmetrischen Kippung eines 
eingespannten Rahmens mit Rechteckquerschnitt (Abb. 2). 


Abb. 1. Das in den Arbeiten von Chwalla und Schoblik untersuchte Abb. 2. Das in der Arbeit von EBlinger untersuchte System mit 
System mit Rechteckquerschnitt (siehe FuBnote 1 und 2 von Seite 394). Rechteckquerschnitt (siehe FuGnote 3 von Seite 394). 


Das Kipp-Problem von Rahmen ist durch diese Abhandlungen jedoch noch nicht erschépfend 
geklart. In allen drei Untersuchungen‘ wird der EinfluB der Normalkraft in der Grundgleichung 
der Torsion 


M;r=—D Ci N/D) UP W Vass 


nicht beriicksichtigt (fiir die Bezeichnungen sehe man Ziff. 2). Ferner steht die Behandlung des 
I-Querschnittes aus. SchlieBlich bediirfen die Berechnung des Horizontalschubes bei statisch 
unbestimmten Systemen und die Ubergangsbedingungen an den Rahmenecken einiger erganzender 
Betrachtungen. 

Ks ist das Ziel der vorliegenden Arbeit, die noch offenen Fragen zu klaren. Zur Kippzahl- 
berechnung der méglichen Lastfalle haben sich die Integralgleichungsmethoden als giinstig er- 
wiesen, weil sie fiir die numerische Rechnung mit Hilfe von elektronischen Rechenautomaten sehr 
bequem sind. Sind die Integrationskonstanten bestimmt, so kann die numerische Lésung den 
Programmierern tiberlassen werden. ~ 

Im folgenden werden die Integralgleichungen der Kippung infolge eines quadratisch verlaufen- 
den Biegemomentes und einer konstanten Normalkraft fiir den geraden Trager mit doppelt-sym- 
metrischem Querschnitt (I- und Rechteckquerschnitt) hergeleitet. Dann wird das geometrische 
Verhalten und das Kraftespiel der freien symmetrischen und antimetrischen Kippung symmetrisch 
ausgebildeter Rechteckrahmen betrachtet. Daraus werden die Randbedingungen und Ubergangs- 
bedingungen hergeleitet. Weiterhin wird untersucht, wie der Horizontalschub ermittelt werden 
muffs, Ferner werden die zur Programmierung fertigen Formeln fiir zwei Lastfalle angegeben. 
SchlieBlich wird das Vorgehen bei der numerischen Berechnung erlautert. 


1 E. Chwalla, Stahlbau 11 (1938) S. 161. 
? E. Chwalla und F’. Schoblik, Proc. of the 5th Int. Congr. of Appl. Mech. 1938. S. 44, 
3M, Eflinger, Stahlbau 23 (1954) S. 53. 


4 Dsgl. in dem wahrend der Drucklegung dieses Aufsatzes erschienenen Beitrag zur Kippung von Rahmen- 
ecken K. Kléppel und W. Protte, Stahlbau 30 (1961) S. 169. 
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2. Vollstaindige Differentialgleichungen und zugehirige Integralgleichungen der Kippung. Die 
beiden fiir doppelt-symmetrische Querschnitte geltenden, gekoppelten gewéhnlichen Differential- 
gleichungen! lauten 

B 


= —_Nu— M,MG)0— K,s— K, (1) 
und 
W D N ip as M ial 
Tia ewes ¢ Do | i M(s) u’ = pd. (2) 


In (1) und (2) bedeutet W die auf den Schwerpunkt bezogene Wilbsteifigkeit, B die Biege- 
steifigkeit fiir die Minimumtragheitsachse, D die de St. Venantsche Drillsteifigkeit, 1 die Stablange 
und z, den polaren Tragheitshalbmesser. Es ist N die iiber die Tragerlange konstante Druckkraft, 
M, ein Vergleichsbiegemoment, M(s) das dimensionslose Biegemoment als Funktion von s,7 = x h/2 
der Abstand der Streckenlast p vom Schwerpunkt, K, und (K, + K,) die Biegemomente an den 
Tragerenden, deren Richtungspfeile mit der Achse des kleinsten Tragheitsmomentes zusammen- 
fallen. Ferner bedeutet 2 den dimensionslosen Eigenwert, a und @ dimensionslose KenngréfBen, 
A, B, C, bis C, Integrationskonstanten, # den Drillwinkel und u* die auf die Stablange bezogene 
dimensionslose seitliche Verschiebung u des Schwerpunktes. SchlieBlich ist z die Langenkoordinate 
und s = z/l. Striche sind Ableitungen nach s. 

Wir setzen uw” aus (1) in (2) ein und erhalten mit der Methode von Lagrange? die beiden Integral- 
gleichungen 


u*(s) = u(s)/d = | (s—1) r,() d&@— 4 [K, 8 + 3K, 8] + 4s+8 (3) 
und ‘ 


(s) = | {Sinh [a(s — 1)] — a(s — #)} r,(t) dt + —4' _% 


jEDQ—o @ k [Sinh (a s) — a s] 


+ ™ {2 [Cosh (as) — 1] — (as)} + {6 [Sinh (as) —as]— (a 8) 
dl ) Cosh (a s)—1] +.“ [Sinh(as)—as = 2 [Cosh (a s) —1]— as) \ 
bape my Ce Cont(as) 1] +2 [Sinh (as) as} + 5 2 [Cosh (as) —1] (0h) 


Al K, K, 
= or a (10 ks? + 5mst-4 3ns8) + 2 (6hst4+2ms+4nsh)] 
+ C, Cosh (as) + C, Sinh(as) + Cg3s + C,, (4) 
- wobei 
Vike i NG Sie Fy 
a= W Sei So. CS a = ED” 
r(s) = —4 |/ 3 [M(s) Hs) + u*(s) NMI , 
ra(s) = 2 {[M%(s) + % B p YB D 1/M3] 8(s) + M(s) u*(s) N1/My} 
d 
E M(s) =k+ms+ns3? 
ist. 


Fiir den Stab mit Rechteckquerschnitt ist in (2) die Grie W = 0, und wir erhalten mit Hilfe 
des Grenziiberganges a—> oo aus (4) die folgende Integralgleichung fiir den Drillwinkel: 


9) =>; [ (9 T(t) Sea, G (10 ks? + 5 ms* + 3 ns°) 
6 


+ 926k + 2mse+ nsh)|+ Cs + Cy. (5) 


1 Die Gl. (1) und (2) sind Sonderfalle der von Chwalla hergeleiteten allgemeinen Kippdifferentialgleichungen; 
vgl. E. Ghent Pee Cu boidiin Wortschaity, Wien 1950, S. 125. Im Gegensatz zu Chwalla zahlen wir 
den Drillwinkel # im positiven Uhrzeigersinne. 

2 Bzgl. deren Erlauterung siehe etwa L. Kraus, Ing.-Arch. 26 (1958) S. 1. 
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Fiir die seitliche Verschiebung u* gilt (3). Bei Beachtung von 

s 1 s ae 
[6.00] = 2609+ [PE] 
0 0 


ist leicht zu erkennen, da (3) und (4) baw. (5) die Differentialgleichungen (1) und (2) erfiillen. 


3. Rand- und Ubergangsbedingungen bei der Rahmenkippung. Wir betrachten einen rechteck- 
formigen symmetrischen Rahmen aus ideal-elastischem Werkstoff. Seine beiden Stiele mégen 
entweder vollstandig eingespannt oder in der Lastebene gelenkig gelagert und senkrecht dazu 
eingespannt sein. Die Querschnitte des Riegels und der Stiele seien I-Profile. Die duBere Be- 
lastung sei symmetrisch (Abb. 3). Sowohl im Riegel als auch in den Stielen mégen sich die Haupt- 


trigheitsmomente wesentlich in ihrer GréBe unterscheiden (E J, > B und E J,> B). Steigert 


ii 
cs 


Abb. 3. System, Abmessungen, Belastung und Lagerungsarten. Es sind EJ, und B bzw. EJ, und B die 
Biegesteifigkeiten bzg. der Maximum- und der Minimumtragheitsachse der Stiele baw. des Riegels (symmetrische I-Querschnitte). 


man die 4uBeren Lasten, so wird der Rahmen zuerst in der Lastebene verformt, um beim Erreichen 
der kritischen Laststufe seitlich auszukippen. Je nach den geometrischen Abmessungen werden 
sich antimetrische oder symmetrische Verformungsfiguren einstellen (Abb. 4). Die Randbedin- 
gungen an den Stielenden sind von der Verformungsart unabhangig. Die Lagerung der beiden 


Abb. 4. Antimetrische und symmetrische Rahmenverformung (schematisch), Koordinatensysteme und Kriftespiel 
(bis zur Riegelmitte) eines antimetrisch bzw. symmetrisch ausgekippten, in Riegelmitte mit 2 N belasteten, 
rechteckférmigen Zweigelenkrahmens (die Kraft 2 N greift um e P oberhalb der Schwerachse an). 


Rahmenstiele soll so beschaffen sein, daB® die dort bei der Kippung auftretenden Torsionsmomente 


My und die Biegemomente K, aufgenommen werden kénnen (Abb. 5a, b und c). Fur den voll- 
stindig eingespannten Stiel sowie fiir die Aushildung des Gelenkes gema Abb. 5a und Abb. 5b 
(die Torsionssteifigkeit der Stirnplatte wird als unendlich gro8 angenommen) gilt 


80) = 0. (6) 
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Dagegen haben wir bei der Ausbildung des Gelenkes nach Abb. 5c 


8'(0) = 0. (7) 
Weiterhin ist bei jeder der beschriebenen Lagerungsarten 
u(0) = 0, (8) 
=al/ 5 
und eee 0) 
(10) 
Abb. 5. Konstruktiv mégliche Ausbildungen des Gelenkes. Abb. 6. Konstruktive Ausbildung der Rahmenecke. 


Setzen wir voraus, da die Rahmenecken nach Abb. 6 konstruiert sind, so 148t sich in guter Nahe- 
rung Wélbverhinderung annehmen, d. h. 

(1) =0 (11) 
und 

8'(1) = 0 (12) 
setzen'. Da die Rahmenecken voraussetzungsgemaB nicht seitlich festgehalten werden, ergeben 
sich die Kontinuitatsbedingungen 


(1) + w(L/h = 0, (13) 
; Al) + u(y = 0 (14) 
u(1) + u(1) = 0. (15) 


Die restlichen Ubergangshedingungen an der Rahmenecke und die Randbedingungen fiir die 
Riegelmitte werden nun fiir antimetrische und symmetrische Kippverformungen getrennt her- 
geleitet. 

a) Antimetrische Verformung des Rahmens. Die Riegelmitte wird sich wohl infolge 
der Biegebelastung in der Lastebene senkrecht nach unten verschieben, jedoch weder verdrehen 
noch seitlich ausweichen. Weiterhin verschwindet dort die Kriimmung der Funktion u. SchlieB- 
lich wird der Querschnitt an dieser Stelle keine Wélbspannungen aufweisen.. Es ist also 


a0) = 0, (16) 
WO) = 05 (17) 
a0} =0", (18) 
w' (0) = 0. (19) 
Aus (1) folgt mit (16) und (18) stets 

Keo 0. (20) 

Aus Griinden des Gleichgewichts miissen die seitlichen Querkrafte gleich grof sein, das heift 
K,=K, All. (21) 


1 Chwalla (s. FuBnote 1 und 2 von Seite 394) schreibt #’’(1) = #’(1) = 0 vor; d. h. die Querschnittsebenen 
der Rahmenecke kénnen sich frei verwolben. 
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Die beiden restlichen Ubergangshedingungen an der Rahmenecke liefert die Grundgleichung des 
Torsionsmomentes+ : 
My = My, — M, M(1) w'(1)/l + (Ky + Ky) o'(L)/t = 0°") W/B — 02) (1 —@) Dit. (22) 
In (22) bedeutet My, das durch die 4uBeren Lasten und die Reaktionen bedingte Torsionsmoment 
und v die Verschiebung des Querschnittsschwerpunktes in der Rahmenebene senkrecht zur Stab- 
achse (Abb. 4). Schreiben wir (22) fiir das Riegel- und das obere Stielende an und gehen wir zur 
Bestimmung von My, vom Gelenkpunkt aus vor, so ergeben sich bei Beriicksichtigung von (11), 
(12), (20) und (21) die beiden noch fehlenden Bedingungen zu 
K,(hjl) + Ky— Nu(1) — My MQ) w'())ll + Ky (l= 8") WIP (23) 
und 
K, — Nu(1) — M, M(1) w'(1)/h + [Ky(h/l) + Ky] o/h = 0°") WR. (24) 
b) Symmetrische Verformung des Rahmens. In Riegelmitte gilt aus Symmetrie- 


grinden 


u'(0) = 0 (25) 
und wegen des Verschwindens der Flanschverwélbungen 
J (0) = 0. (26) 
Wirkt dort eine Last 2 P in der Entfernung ep oberhalb des Querschnittsschwerpunktes, so muB 
W''(0) W/B = H0) P ep (27) 
sein [vgl. (22)]. Bei der Beanspruchung des Riegels durch eine Streckenlast haben wir statt (27) 
WY (0) =O. (28) 
Aus Symmetriegriinden verschwinden die Querkrafte K,/l und K,/h, das heiBt 
: K,=0 (29) 
und 
K, =0. (30) 


Das Stieleinspannmoment K, ergibt sich aus der Gleichgewichtsbedingung (Verschwinden der 
Momente beziiglich der durch die beiden Gelenke gelegten Achse) zu 


Ky = P [ul0) + ep HO +1 p ful +9 OO] at. @l) 


Mit Hilfe von (30) schreiben wir (22) fiir das obere Stielende an 
K,— N [u() —u(})] — My MQ) wh + Koh = 8") (32) 
und haben damit alle Bedingungen zur Bestimmung der Integrationskonstanten gewonnen. 
c) Sonderfall: Rechteckquerschnitt. Es ist W = 0 baw. a =o, und es gelten (3) und (5). 
Es sind zwolf Konstanten aus Rand- und Ubergangsbedingungen zu ermitteln. Diese Bedingungen 
sind fiir den I-Querschnitt in den Abschnitten a) und b) angegeben und kénnen zum Teil unver- 
andert iibernommen werden. Es entfallen alle Aussagen iiber die Verwélbung und die sekundaren 


Normalspannungen [(6), (7), (11), (12), (17) und (26)]. Statt (22) gilt nun 
My = — 8'(1) (1—o) Dil, (22a) 
und deswegen sind (23), (24), (27) und (32) abzuandern in 
K,(hl) + K,—Nu(1) —M, MQ) w'()/t + K, ol + 6) (le) Dl =0, 23a) 
K, — Nu(1) — M, M(1) w'(1)/h + [K,(h/l) + K,] 0'(1)/h + 81) l—0) Dih=0, (24a) 
0(0) Pep + #0) (l—e) D/Jl=0 (27a) 
und 
K, — N [u(0) — u(1)] — My M(1) a'(1)/h-+ Ky o'()/b + HC) (1—A) Dll =o. (82a) 
SchlieBlich gilt statt (28) 
O10) —0.. (28a) 
* Siehe etwa die FuBnote 1 auf S. 21 der im Forschh. a. d. G. des Stahlbaues Nr. 6. Berlin 1943 verdffent- 


lichten Arbeit von E. Chwalla, sowie Gl. (B12) auf S. 32 des als Bull. Nr. 10 f i ildi 
; f : RING rom The Div. 
Statics and Struct. Eng. at The Royal Inst. of Techn., Stockholm erschienenen Hiehes Yoh 0. Paste 
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4, Ermittlung des Horizontalschubes. Fiir die Verschiebung v (Abb. 4) gilt? 
Dae Vazex = Nos a P(e) S (33) 


Das Gleichungspaar (1) und (2) bildet zusammen mit (33) die Grundgleichungen der raumlichen 
Stabilitat von Staben mit doppelt-symmetrischen Querschnitten. Da voraussetzungsgem4B B,, 1. 
> B ist, tritt ein Stabilitatsversagen des Rahmens in seiner Ebene nicht ein. Dennoch ist (33) 
bei statisch unbestimmten Systemen bedeutsam. Denn aus dieser Grundgleichung errechnen 
sich der statisch unbestimmte Horizontalschub H und die Stieleinspannmomente. Da in (33) der 


2= 
DES 


#2 ——+j 


Abb. 7. System und Verlauf des Horizontalschubes; - - - - - qualitativer Verlauf nach der Theorie erster 
Ordnung [H,) = M/h(1 + J,h/3 J, 1); qualitativer Verlauf nach der Theorie zweiter Ordnung. 


Kinflu8 der Normalkrafte auf die Biegung beriicksichtigt wird, ergeben sich die statisch Unbe- 
stimmten stets aus einer transzendenten Gleichung fiir den Horizontalschub H. Tragen wir H 
ber der auferen Belastung auf, so ist zu erkennen, daB die gegenseitige Abhangigkeit nicht mehr 
linear ist. So wird etwa der Horizontalschub H eines Zweigelenkrahmens, der an den Rahmen- 
ecken mit den gegengleichen Momenten M belastet ist, fiir gréBere M kleiner baw. gréBer als der 
nach der Theorie erster Ordnung berechnete Wert, wenn H im Riegel als Zug- bzw. als Druck- 
normalkraft wirkt (Abb. 7). Bei Verwendung des Horizontalschubes der Theorie erster Ordnung 


(jae z/( 4 zal (34) 


erhalten—wir also falsche Kippzahlen. Um dem zu entgehen, empfiehlt es sich, den Horizontal- 
schub durch eine Potenzreihe 
H =a, M + a, M? + a, M? + a, M* (35) 


anzunahern und mit (35) in die Stabilitatsberechnung einzugehen. Die Koeffizienten a,, a,, a3 
und a, in (35) sind von den Tragheitsmomenten und den Langen der Stiele und des Riegels ab- 
hangig und werden aus der aus (33) resultierenden transzendenten Gleichung gewonnen. 


5. Lastfall: Rahmeneckmomente. Wir untersuchen einen Rahmen mit Rechteckquerschnitt, 
der an den Ecken durch zwei gegengleiche 4uBere Momente M helastet, in der Lastebene gelenkig 


Abb. 8. Systemskizze des Rahmens mit Rechteckquerschnitt und Eckmomenten M. 


gelagert und senkrecht dazu an den Stielenden eingespannt ist (Abb. 8). Die transzendente Glei- 
chung zur Bestimmung des Horizontalschubes H lautet 


y= (a) ects + oP (36a) 
beziehungsweise 


ver (ss) é Ctgh (e) + 4 > (36b) 


1 Siehe die in FuGBnote 1 von S. 395 genannte Arbeit von E. Chwallia. 
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wenn H im Riegel als Druck- bzw. als Zugnormalkraft wirkt. In (36a) und (36b) bedeutet y = M/A, 
w= JolJ,, x=1fh und e2 = HP/E J,. Fir x =1 und verschiedene -Werte sind die (y, ¢)-Kurven 
nach (36a) und (36b) in Abb. 9a und b eingezeichnet. Fiir ¢ = 2/2 ist y = 1/x (Abb. 9a) und 
unabhingig von yw. An diesem Punkt ist die Normalkraft im Riegel gleich der Eulerlast Py 
— 72 E I,/(21)2, die Verformungen werden unendlich groB, das Tragvermégen des Rahmens in 


seiner Ebene ist erschépft. Fir ¢ = 0 gehen (36a) und (36b) in den Wert yo = [(u/3 x?) + 1/x] 


MHL M/Ht 


| System: 
System M M 
M M 7 7 as 
a r ; : uy t 
z aS) 
i L 5 J, 
6 + ——Z1 ——+ 
+2] ——+ iy, 
Bezeichnungen: 
Bezeichnungen: y=M/HL 
Nill peee 
45 pi s side, 
forest 2 x=l/h=1 3 
pad] J 


Abb, 9a und b, (y, ¢)-Kurven nach (36a) und (36b) fiir x = 1. 


der Theorie erster Ordnung iiber. Es ist also in Abb. 9a nur der Bereich 0 < ¢ < z/2 interessant. 

Aus Abb. 9a und 9b erhalten wir die (H, M)-Kurven bei konstanten J,, J, h und / durch einfache. 
Umrechnung (Abb. 7). Dafiir benétigen wir in Abb. 9a und 9b den Bereich kleiner Werte ¢. Denn 

die den gréBeren ¢ entsprechenden Wertepaare H und M erzeugen im Rahmen Spannungen, die 

oberhalb der FlieBgrenze liegen. 

Da im Riegel ein konstantes Biegemoment M,= M—H h wirkt, lassen sich fiir ihn die Lé- 
sungen von (1) und (2) geschlossen angeben. Wir beschrainken uns dabei auf diejenigen Falle von 
M, bei denen stets 0< <1 baw. > 0 gilt, wenn im Riegel eine Druck- bzw. eine Zugnormal- 
kraft wirkt. Mit N = H (positiv, wenn Druckkraft), 


M,=M—Hh, M()=1, W=0, y»=0, A=MUVBD 
und 5 
= [(H P/B) + #/(1 —e)] 
lauten die Lésungen von (1) und (2) fiir den Riegel 
B(s) == a {A, cos ks-+ A, sink s + [1— (1—e) (k/a)?] (As Aj 


und 


| 
AVBD 


(K,s+K,) (37) 


u*(s) = A,cosks + A,sinks + A,s + A,. (38) 


XXX. Band 1961 L. Kraus: Problem der freien Kippung rechteckiger Rahmen 401 
Mit 7 = 0, M,= Hh, A= H h2/B D, m=—lundk=n=N 0 = 0 ergeben sich aus (3) 
und (5) die Integralgleichungen und deren erste Sea fiir } und u* des Stieles zu 
Hs) = 2 j (¢ — s) M%(t) S(t) de + Ee (ist-+ 2K, 8) +654 C,, (39) 
is) =—# [ Me() 8) de + 2% (K 24K, ")4€ 40 
, /BD £ 3 2 ) Le ( ) 
u*(s) =A)/D/B [ (e—s) Mi) He) & —*. (K, 8 +3K,8)+ 454+ B 
s) M(t) H(t) de 5g (Kit +3Ki8%)+ 45+ 8B (41) 
0 
und 
u*(s)! = — AV DiB | M(t) H(t) eae ck s)+ ve (42) 
a) Antimetrische Kippform. Aus (16), (18), (19), (37) und (38) erhalten wir 
A, = Ay=K,=0. (43) 
Ferner ergeben (8), (9), (10), (21) und (39) bis (42) 
A=B=C,=0 (44) 
und 
K, = K, hil. (21a) 


Mit Hilfe von (37) bis (42) folgt aus (14), (15), (13), (23a) und (24a) bei Beriicksichtigung von 
(43), (44) und (21a) das lineare Gleichungssystem zur Bestimmung der restlichen fiinf Konstanten 
1 


Ah? : = ‘/ 
ae Vas at+kcoskA,+A,+C,=4 Be Ore 
— oe Ky — Ky + sin ke Ay + Ay = A) D/B (Hl) j (1—1) MW 0 at, 
weet eet ad hee oe | Be 2 BID ieee 
eee gt eee, [oe a(z)|4= 
=1/D/B J M(t) He) dt, | (45) 


AEP 


fon 1/50 PP Kt Ky EP o_o (th4.=0, 
Voie eee 


ne 


j M2(1) 3(t) de + i t M(t) (t) a. 


Weil die Normalkraft fehlt, ist die Integralgleichung (39) von (41) unabhangig und reicht mit den 
aus (45) zu ermittelnden Konstanten K,, K, und C, allein zur Berechnung der Eigenwerte 4 und A 
aus. Um (39) in der fiir die numerische Penecn bequemen Form 


~” = [ (t — s) M(t) H(t) dt + X,s*+ X, 8? + Xs (46) 
0 : 
schreiben zu kénnen, miissen die GréBen 
Bee es py a Sehoh binant, £47) 
12A1/ BD 61/BD i 


28 
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anstatt K,, K, und C, als Unbekannte in (45) eingefiihrt werden. Bei Erweiterung der beiden 


letzten Zeilen mit h2/ 124 1/BD D ergibt sich somit aus (45) in Matrizenschreibweise 


Aq= (45a) 
[22] [ke cos k | [1] 
[0] [sin k] 1] 
ae | 9 Rane] (EP a —-a0(3)]] 
0 io |—om ag epee I 
ayo/Bim2y| (0) [0] 
[2] [72] 


[2 4h") VD/B| [— 3 any V0/B] 
E 12 (ih)> \ 67 yi8 | [— 6a p/B| 
con 3. DB je P| [h/2.1 


h+@niyaoses)  [w2nifeoises,| 


2 { (1—t) Mt) Ho) at 


oe 
a 1 DIB (h/t) [ 1—1) MQ) O@ at 
q= A, 9 t= iyore [ Mi 3 d 
x, ped (t) H(t) dt 
x, : 


1 1 
2 DIB (h/12 1) | [i HO ae + f MH DE i 
6 6 
und v'(1)/h = — v'(1)/l = H h?/3 E J, ist. 
b) Symmetrische A ABA dts Aus (25), (28a), (37) und (38) b 

Ree (25), (28a), (37) und (38) bzw. (8), (9), (10), (37) bis (42) 
bzw. A,= A,=0 (48) 
A=B=C,=0. (49) 

Mit Hilfe von (29) bis (31) und (37) bis (42) liefern die Uber i 
| gangsbedingungen (32a), (13), (15 
und (14) das folgende Gleichungssystem zur Bestimmung der restlichen fe aS pa 


Po4 kK; EEL! 


(cos k—1) A, = (D#/h) , if M(t) 9(t) dt + M(t) H(t) i, 


~ 


ae spo lnie 1B cosh A, + BPE, —a—e(*)| A=i|/2 | Mo ona, 


A 
(50) 


cos k A, +A, = 4 (h/l) 2 f (1—t) M(t) He) dt, 


1 
C,—ksink A, =# [ (11) M(t) H(t) dt. 


0 
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Nun erweitern wir die erste Zeile von (50) mit h/ D und fihren statt C, und K, die beiden dimensions- 
losen Unbekannten 


C, 


me 


_ Kl 

AVBD 
ein. Entsprechend (46) laBt sich nun die zur Berechnung von A und / allein notige Integralgleichung 
fiir } wie folet schreiben: 


be, und —sX, 


(51) 


OWS ={ (t — s) MX%t) Di) dt + X, 8, (52) 
und das Gleichungssystem (50) lautet in Matrizenform 
Dabei ist ee Sey 
— h\ AVBD Sri 
[22] (7) 2 | fe A \z (cos k — | [0] | 
ie 1 ea ee eo GI, 
[0] [0] [cos k] [1] 
[?] [0] [— k sin k] [0] 
22 | il M(t) A(t) dt + f M(t) H(t) a 
X, ' Ss 
- A \2 f Mo 5 a 
oe und — ae 2 
A rT sy) — +t) Ma) 8d 
| A(h/2) \2 i (1 —t) M(t) H(t) dt 
A, 1 


22 f (1 — t) M(t) A(t) dt 
f 
6. Lastfall: Einzellast in Riegelmitte. Wir untersuchen denselben Rahmen wie in Abschn. 5. 
Statt der Rahmeneckmomente M sei er jedoch mit einer an der Riegelschwerachse angreifenden, 
lotrechten Last 2 P beansprucht (Abb. 10).. Den Horizontalschub ermitteln wir aus der trans- 
zendenten Gleichung 


Gi) — wine sia Vrlu=9) _ (iy) + y [ectg(e) —1] =0 (53a) 
cos (Vy/u x é) 


Abb, 10. Systemskizze des Rahmens mit Rechteckquerschnitt und Einzellast 2 P. 


bzw. 
(53b) 


(/y) + Vylw¢ Sinh Voluxe) _ (1/y) +y[e Ctgh(e) —l1]=0, 

Cosh (Vy/u x €) 
wenn P im Rahmen Druck- bzw. Zugnormalkrafte hervorruft. In (53 a) und (53 b) ist y = H/P, 
p= Jo/Jy, * = Yh und 2 = PRE J. 
28* 
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Abb. lla gibt fiir x = 1 und verschiedene y den Verlauf der (y, ¢)-Linien nach (53a) wieder 
Auch hier diirfen die Drucknormalkrafte in den Stielen und im Riegel die Werte der Eulerkrafte 
nicht erreichen [(P h?/E J,)< 2? und (HP/E Je) < (x/2)?]. Das heiBt, daB fiir die y-Werte 
stets y << 2x/7 unde < gilt. Die Kurven wu = konst. beginnen auf der Ordinate mit dem Wert 
1,5 x/(3 + u/x) und haben dort eine waagerechte Tangente. Diese Kurven verlaufen 


HP Yorena = 22 (Knicken des Riegels) 
= ; : HIP 
06 eal a System: Es 
a] 2p | 
~S 
1 
ép. 
05 
per fn, 
Bezeichnungen: ea 
y=H/P > Bezeichnungen: 
oy j = e PRE, AS y=H/P 
x=l/h=1 g 
b=b/d; § 
< 
eed 
93 i = 18 
fe: 
el > 
ie 
y | 
02 | 
| 
| 
Of ee I | 
| 
| 
|) [Pie 
ANVE4 
0 05 40 15 20 30 0 05 40 45 40 30 
a b 


Abb. lla und b. (y, ¢)-Kurven nach (53a) und (53b) fir x = 1. 


a) fiir kleine Werte ~ monoton steigend, um die Gerade y = 2 x/z bei ¢ = /(x/2 x)? 4 zu er- 
reichen; sie verlaufen 

b) fir groBe Werte ~ monoton fallend und schneiden die Abszisse im Punkt A(z; 0). 

Die Kurven c) mit mittleren y-Werten steigen zuerst maBig, erreichen ein Maximum und 
miinden dann in den Punkt A. Die unter a) und c) genannten Kurven zeigen in ihrem Verlauf, 
daB der sich?aus der Theorie zweiter Ordnung ergebende Horizontalschub Hy gréBer als der nach 
der Theorie erster Ordnung berechnete Wert Hy ist. Beziiglich der Horizontalschiibe wirken also 
die zugehérigen Rahmen! wie Bogen. Genau das Entgegengesetzte (namlich Hy > Hy) ist bei 
den Rahmen der Fall, zu denen die unter b) beschriebenen Kurven gehéren. Hier ist also ein Ver- 
gleich mit den entsprechenden Horizontalschiiben des Kabels von Hangebriicken angebracht. 
Fir kleine Werte x streben auch die Kurven mit kleinen zu dem Punkt A. Das bedeutet, daB 
selbst bei kleinen Riegeltragheitsmomenten solche Rahmen in ihrer Ebene durch Knicken der 
langen Stiele versagen. 

Abb. 11b zeigt den Verlauf der (y, ¢)-Linien fiir x = 1 und verschiedene fs nach (53b). Die 
Kurven mit kleinen y fallen monoton. Es ist Hj; << Hy, und es gilt die Hangebriicken-Analogie?. 
Jedoch fiir gréBere ~ steigen die Kurven (Hj > Hy; die Rahmen verhalten sich wie Bégen). 


1 Ebenso die Rahmen der Abb. 9a. 
* Desgl. auch bei den Rahmen der Abb. 9b. 


/ 
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Zur Umrechnung in (H, P)-Kurven benétigen wir auch in Abb. lla und b nur den Bereich 
kleiner Werte ¢. Fiir gréRere ¢ ist namlich das Tragvermégen des Rahmens praktisch durch Er- 
reichen der FlieBgrenze erschépft. 


Mit den Gréfen 4 = PP//BD, 1= Hh/VBD, 0 =i? H/D, 6=# P|D, k= 1— Hh/P1, 


m=m=—], n=k=n=0, N= HA N= PM, = Pl, M,= Hh, mit den Konstanten 
A= B=C,=0 und K, = K, h/l [siehe (8), (9), (10) und (21)] und schlieBlich mit den Abkiir- 


zungen 
r,(s) = —A//D/B [M(s) 9(s) + u*(s) H/P], 
r2(s) = 2 [M*(s) 0(s) + M(s) u*(s) H/P], 
(3) = —A /D/B[M(s) (6) + a*(s) P/H], 
72(s) = 4° [M?(s) B(s) + M(s) u*(s) P/H] 


(54) 


schreiben sich die Integralgleichungen (3) und (5) fiir u*, 9, w*, und @ und ihre ersten Ableitungen 
wie folgt: 


Ss 


u*(s) ea r() dr— © (K, 8 +3K, 9%) + As+B, (55) 
0 
u*(s)! = fn dt —* (Ks +2K,8) +4, (56) 
; 0 
s : K, 
AO erga ra(t) dt— aes ie 2h + ms!) +*2(3ks?+ ms%)+ Cs+ Cy, (57) 
1 ; At ie K, 
90) =—1=5 | 04 — Te RG HL 6 (Bke4+2ms)+A2Qks+ms)+C, (58) 
u*(s) = | = H7,() dt —~_[K, (hil OK, #4), (59) 
0 
a*(s)! = i 7.0) &— "TR b/s +2 Kal, 
; 2B7 (60) 
56) = _[ ¢—9 79 + —* |B (t) 0+ Bala Gs, (61) 
i—¢ /BD(1—9) 
0 
Foy — [70 geet Pa(4) #4 a oa} 4 Cx (62) 
ee /BD(1—o) 


Auch hier gelten wieder einschrankend 
1—oe>0 und 1—e>0, 
wenn P in den Rahmenstielen als Drucknormalkraft wirkt. 
a) Antimetrische Kippform. Aus (16), (18) und (19) folgt mit (55) und (57) wiederum 
C,=B=K,=0. oe 63) 
Die Gleichungen (55) bis (62) liefern in Verbindung mit den Bedingungen (14), (15), (13), (23a) 
und (24a) das zur Bestimmung von K,, K,, C,, C, und A erforderliche Gleichungssystem 
: Ng (64) 
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Dabei ist 


KHZ) bya] © [ve b-GYecbaggat| Forme 
A 
6 


x [way] fy] [()-2V5e—»—(i)35—+] [Fe /83 | 
V2] =o oo [es yeh+A| fl 
[oe hte fe 1) Fal) dt 
X | (1—1) r4(t) a+ (7) | (1—1) 7,(t) de 
q= z oe | fr <> sy gaeee 


0 


11 fae u— [em dt 


X= GQ(l—o/@, = X,=AYB/D/i, — Xs= C,(1—o)/?*, (65) 
X,=K,layBD, X, = K,h/AVBD 
und « = v'(1)/h = — v'(1)/l = y [1 — e ctg(e I fiir gréBere P nach der Theorie zweiter Ordnung 


[siehe (53a)] bzw. « = H h?/3 E a wenn der Horizontalschub H bei kleinen P und sehr groBen 
Biegesteifigkeiten E J, und E J, nach der Theorie erster Ordnung ermittelt wird. Bei Beriick- 
sichtigung von (65) haben die Integralgleichungen (55), (57), (59) und (61) die folgenden, fiir die 
numerische Rechnung bequemeren Formen: 


ur(s) (BPD/a = | (s—9 © yap de —X, 2+ X,0, (66 
0 


He) Aol = | 6-9) SP de —TEOkS—s) + Xys, (67) 


u*(s) /B/D/a = [ @—95? B/D at — my VB DIBD Gs eke (68) 


A(s) (1 — 9/22 = | (¢— 9 HO d+ = VB FFE ae ria (69) 
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b) Symmetrische Kippform. Aus (25) und (56) sowie aus (27a) mit ep = 0 und (58) ergibt 
sich 


A=C,=0. (70) 
Weiterhin ist nach (29), (30) und (31) 

K,=K,=0 (71a) 
und 

K,= PIB. (71b) 


Die GréBen C,, B, K, und C, errechnen sich aus dem durch (13), (15), (32a) und (14) mit Hilfe 
von (70), (71) und (55) bis (62) gebildeten Gleichungssystem 


Sq=t, (72) 
mit 
oS alee =) a | [0] x 
oF] Yo : 


hes = Aa Nia oe : q= , 
© VCE) b-(y2 2] [ren] : 
o-(aye] ©  Ettevz]] | 


Et — 6) [rout | (¢—1) ra(t) a 


0 


fo—one at (F )fa-s ia 
i) oe { F(t) a—il/2 | { 7,(0) a(t) | (1— 1) r,(2) al} 


1 


B/D 
ae 


I 


1 
jee e5 (1 — 9) | r4(t) dt + | (¢— 1) 7,(0) i 
aa _ plB/D . 
X,=C, z° Kei Bae, PA 
_ _Kl ee eee 
ee ED: eae je 


und denselben « wie in Abschn. 6a. Mit Hilfe von (73) lauten die Integralgleichungen (55), (57), 
(59) und (61) wie folgt: 


v2(9 EB = | (6 — 90 yD i Xe (74) 
Hs) gt= [vo Qa % (BkF—)+X, (75) 
ant) 12D _ [o —)2O/B B/D dt — At ie Tae X,8, (76) 
fo Zt fo 7 dat a \F agg me eyes 
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7, Numerisches Verfahren. In den rechten Seiten t der Gleichungssysteme 2 q = t und © q = t 


(vgl. Abschn. 5 und 6) kommen die unbekannten Funktionen @ bzw. u*, 0, u* und @ unter be- 
stimmten Integralen vor. Die Lésung muB mit Hilfe der Kehrmatrizen geschehen: 


= iy und G¢ =O". (78) 
Die Simpsonformel 
. t+1 
[ toa= = Sak (79) 
i=1 


0 


stellt die Grundlage der numerischen Rechnung dar, die wegen der komplizierten Bauart der 
Integralgleichungen nur auf elektronischen Rechenautomaten erfolgen kann. In (79) ist t die 
Zahl der Intervalle, g; der Simpsonkoeffizient und J; der Integrand an der Stiitzstelle t;. Durch 
Anwendung von (79) auf die Integralgleichungen der Abschn. 5 und 6 ergeben sich homogene 
Gleichungssysteme, die als Matrizeneigenwertprobleme, wie folgt, lauten: 


[%(6) —C] xy = 0. (80) 
Die Elemente des Vektors ¢ bestehen aus den 9(s;) bzw. aus den u*(s,), O(s;), u*(s;) und P(s;) an 


den Stiitzstellen s; (vgl. Abschn. 5 und 6). Weiterhin ist 6 der Matrizeneigenwert, © die Einheits- 
matrix und & die aus dem Kern bzw. den,,Teilkernen“ K;,, gebildete Koeffizientenmatrix (Abb. 12). 


Kou, Kap Kaz Kee 


Abb. 12, Matrix @ und ,,Teilkerne“ K,- 


In den K;; setzen sich die Koeffizienten links von der Hauptdiagonale aus zwei Anteilen zusammen: 
aus den von den X,, und von den Integralen mit den Grenzen 0 und s herrihrenden Gliedern [vel. 
etwa (74)]. Dagegen bestehen die Koeffizienten rechts von der Hauptdiagonale nur aus den durch 
die X, bedingten Gliedern. Die 6 erhalt man durch Probieren. Fiir angenommene 6, wird mit 
Hilfe des Gaufschen Eliminationsverfahrens der Wert AN der Nennerdeterminante von (80) 
berechnet und iiber 6, aufgetragen. Der zu AN = 0 gehérende 6-Wert ist der gesuchte Matrizen- 
eigenwert. 


Nun zeigen wir die Zusammenhange zwischen 6 und A, A, 0, 0, k2, k der Abschn. 5 und 6. 
Im Abschn. 5 gilt 


4 
Has Sa. Mr. 
r=1 


Wir wahlen als Kigenwert 


Mi//BD=6. 
a 4 ee, 
Mit a, = a, h(/B Dil)'—* und R(6) = 2 ao" wird H = (//B D/h1) R(6). Damit ergibt sich 


A=5—R(5), 1 = (h/l) VBD/BDR(3), 9 = (ih) /BID R()| 
| (81) 


und 


i = jh) YDB R(@) + © ROK) 
Im Abschn. 6 haben wir 
Hi = x a, BF. 


r=] 
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: =a = ——- == 4 
Es sei 06 =A = PP//B D. Mitye = a. (/B D/I)’ r2//B D und R(6) = Ya, 6" erhalten wir H 


r=1 


= (//B D/h) R(d). Damit ist 


1=R(6), e=(,/hVBD/DR’), 9 (1) B D/D* >| a 


und 
k = 1—(I/h) /B D/BD R(0)/5. 


Bei sehr schmalen Rechteckquerschnitten sind die Kipplasten sehr klein, so daB man den 


‘Horizontalschub in guter Naherung nach der Theorie erster Ordnung ermitteln und 9 und 0 ver- 


nachlassigen kann. Es entsteht dann die lineare Beziehung R(5) = a* 6, und wir erhalten aus (81) 


be eee): ey PDE Det 3, =| oe 
k2 = 6[(I/h) /D/B a* + (1—a*)? 6], 
bzw. aus (82) 
A=a*d, o@=0=0 wd k=1—(h) VBD/BDa*. (82a) 


AbschlieBend ein Wort iiber die Wahl der Koeffizienten g; [vgl. (79)]. Bei Rechteckquerschnitten 
verwendet man fiir die ,,Teilkerne“ K,,, Kgs, K.- und Kg j (Abb. 12) wegen des Knickes bei 
s=t (Hauptdiagonale!) am zweckmafigsten die Trapezformel. Bei I-Querschnitten gilt das 
nur fiir K,,, und K;-. Fir die restlichen ,,Teilkerne“‘ K;, empfehlen sich auch die andern Simpson- 
formeln. Lediglich die Rechteckformel [in (73) sind alle g; = 1] sei ausgeschlossen. Denn sie 
nahert die Integrale oft schlecht an, und dadurch werden die Eigenwerte 6 zu ungenau. 

Der Verfasser dankt den Herren Prof. Dr.-Ing. K. Marguerre und Dipl.-Ing. E. Ludwig fir 
wertvolle Hinweise und freundliche Ratschlage. 


(Eingegangen am 17. Oktober 1960.) 
Anschrift des Verfassers: Dipl.-Ing. L. Kraus, Wiesbaden, Idsteinerstrafe 11. 
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Zwei Extremalprinzipien der irreversiblen Thermodynamik 
Von H. Ziegler 


Der Erinnerung an Hans Kauderer gewidmet 


1. Einleitung. In der Thermodynamik irreversibler Prozesse spielen, soweit es sich um lineare 
Vorgange handelt, die Reziprozitatsrelationen von L. Onsager! die zentrale Rolle. Diese Relationen 
sind schon von ihrem Begriinder, aber auch seither? und insbesondere im Zusammenhang mit 
ihrer rheologischen Bedeutung® als Ausdruck gewisser Extremalprinzipien gedeutet worden. Im 
Anschlu8 hieran sollen im folgenden zwei weitere Extremalprinzipien formuliert werden, die mit 
den Onsagerschen Beziehungen aquivalent sind und den Vorzug haben, besonders einfach sowie 
physikalisch einleuchtend zu sein. 

Man kann diese Prinzipien im Gegensatz zu den Reziprozitatsrelationen auch auf nichtlineare 
Prozesse anwenden. Sie fiihren in diesem Falle auf die vom Verfasser* vorgeschlagene Verall- 
gemeinerung der Onsagerschen Theorie. 


2. Die Hauptsitze. Es sollen hier einfachheitshalber sowie im Hinblick auf die Anwendungen 
in der Mechanik der Kontinua nur Systeme betrachtet werden, deren Zustand durch die mechani- 
schen Lagekoordinaten x, (k= 1,2,...,n) sowie die Temperatur ? > 0 beschrieben wird. 
Die an einem solchen System geleistete Elementararbeit ist, wenn hier und im folgenden die 
Summationsregel verwendet wird, durch 


dA = Xi, dx, (2.1) 


gegeben, wobei die X,, die zu den x, gehérenden verallgemeinerten Krafte bezeichnen. 

Der Inhalt des ersten Hauptsatzes der Thermodynamik besteht in der Aussage, da das 
System eine durch die Zustandsvariablen x, und # bestimmte Energie U(x,, #) besitze, 
deren Inkrement sich gemaB 


dU =dA+dQ (2.2) 
aus der mechanischen Arbeit (2.1) und der Warmezufuhr dQ zusammensetzt. 


Nach dem zweiten Hauptsatz gibt es eine weitere Zustandsfunktion, die Entropie S(x,, 9), 
deren Inkrement der Bedingung 


i dS = dQ (2.3) 


gentigt. Gilt in (2.3) das Gleichheitszeichen, dann wird die zugehérige Zustandsanderung als 
reversibel bezeichnet, andernfalls als irreversibel. Man kann (2.3) auch in die Aussage kleiden, 
daB sich der elementare Entropiezuwachs gemaB 


dS = dS" + dsi (2.4) 


aus der Entropiezufubr 


é il 
ds" = dQ (2.5) 


von aufen und der Entropieproduktion dS‘ im Inneren des Systems zusammensetze, wobei 
je nach der Art des Prozesses 
dS' = 0 (2.6) 


ist. 
3. Irreversible Prozesse, Die Elementararbeit 14Bt sich nach (2.2), (2.5) und (2.4) in den Formen 
dA = dU — dQ = dU — $ dS" = dU— 6 dS + 8 dSi (3.1) 


1. Onsager, Phys. Rev. 37 (II), (1931) S. 405 und 3 5 3 ‘ . B.G. imi. 
S gee eer nee 7 (94s) = Le ), (1931) 8 (II), (1931) S. 2665. Vgl. auch H. B. G. Casimir, 

2S. R. de Groot, Thermodynamics of Irreversible Processes, Amsterdam 1952. ; 

° M.A. Biot, J. appl. Phys. 25 (1954), S. 1385; Phys. Rev. 97 (1955), S. 1463; Variational and Lagrangian 
Methods in Viscoelasticity, [UTAM-Colloquium Madrid, S. 251, Berlin 1956. 

* H. Ziegler, Ing.-Arch. 25 (1957), S. 58; Z. angew. Math. Phys. 9b (1958), S. 748. 
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oder ausfiihrlich mit 


i: as eU aS . 
X;, Ax, oe 0) zs) dx, + (53 — 38) dd + 9 dS? (3.2) 
anschreiben, Bei einer Zustandsanderung, welche die Lagekoordinaten x, unverandert laBt und 
also in einer reinen Erwarmung oder Abkiihlung besteht, mu8 die Summe der beiden letzten Terme 


verschwinden. Das ist aber sowohl fiir Zu- wie fiir Abnahmen von ? wegen (2.6) nur méglich, 
wenn 


6U os 
wo aa ° (3.3) 
ist. Hieraus folgt, daB das letzte Glied in (3.2) die Gestalt einer Elementararbeit 
& dS' = Xj, dx, = dAi (3.4) 
besitzt. Man nennt sie die Dissipationsarbeit und hat nach (2.6) je nach Art des Prozesses 
dA 0". (3.5) 
Mit (3.4) und (3.3) gewinnt man aus (3.2) 
Fite Rep at OU OS .b)6 
Xj, = X,— Xi = an, O aS oe ee & S) (3.6) 
oder 
deo 
Pn re mit F=U—#S (3.7) 
und daher 
5 ae oF 


Jede Kraft X;, stellt sich damit als Summe von zwei Beitragen dar. Der eine, Xj, leitet sich ge- 
maf} (3.7) von der freien Energie F her, die mit U und S eine Zustandsfunktion ist. Er hangt 
also nur von den Zustandsvariablen x,, ? ab, und da er nach (3.4) keinen Beitrag zur Entropie- 
produktion dS' leistet, reprasentiert er gemaB (2.6) den reversiblen Teil des Prozesses. Der andere, 
Xj, bestimmt die Entropieproduktion und damit den irreversiblen Teil der Zustandsanderung. 
Er wechselt nach (3.4) und (3.5) mit dx, sein Vorzeichen, hangt also jedenfalls von der verall- 
gemeinerten Geschwindigkeit x, ab. 


4. Grenzfalle. Es ist niitzlich, wenn auch in diesem Zusammenhang nicht unerlaBlich, noch 
einige Begriffe zu definieren. Je nach der GréBe des Anteils dA‘ in dA kann man von Prozessen 
sprechen, die mehr oder weniger stark irreversibel sind. In diesem Sinne kann man ein System, 
bei dem fir beliebige Wahl der dx, die Dissipationsarbeit dA‘ = 0 ist, als reversibel und ein 
solches, bei dem stets dA‘ = dA, mithin dA’ = 0 ist, als vollkommen irreversibel bezeichnen. 
Beispiele fiir die beiden Falle sind das ideale Gas bzw. die inkompressible zahe Flissigkeit. 

In der Thermodynamik sind bis vor kurzem ausschlieBlich reversible Systeme behandelt worden. 
Erst die neuere Literatur hat sich irreversiblen Prozessen zugewandt, und zwar beschrankt sie 
sich meist stillschweigend auf vollkommen irreversible Systeme. Hier verschwinden die Xj, iden- 
tisch; nach (3.7) ist also die freie Energie nur eine Funktion der Temperatur, d. h. es gilt 


U—¢S = F(#). (4.1) 
Leitet man beide Seiten nach # ab, so erhalt man mit Riicksicht auf (3.3) 
oF 
S tie 00 . (4.2) 


Somit sind auch die Entropie und die Energie Funktionen der Temperatur allein und daher jede 
durch die andere gegeben. In der Tat folgt aus (4.1) und (4.2) 


U=0S— | Sap (4.3) 


a 


und umgekehrt 
dU 
sa {=. (4.4) 


1 Der obere Zeiger r baw. i bezeichnet den reversiblen bzw. irreversiblen Anteil der betreffenden GroBe. 
Man kann diese Unterscheidung entweder in den Lagekoordinaten oder in den verallgemeinerten Krafien 
treffen. In der zuletzt zitierten Arbeit (FuGnote 4, S. 410) hat der Verfasser den ersten Weg gewahlt, hier den 
bereits von B. T. Chu (Deformation and Thermal Behavior of Linear Visco-Elastic Materials, Report Brown 
University, Providence, R.I., 1957) beschrittenen zweiten Weg, der in diesem Zusammenhang einfacher 
erscheint. 


412 H. Ziegler: Zwei Extremalprinzipien der irreversiblen Thermodynamik Ingenieur-Archiv 


5. Die Onsagerschen Relationen. Die vorstehenden Uberlegungen beruhen auf den beiden 
Hauptsatzen sowie auf der Annahme, daB der Zustand des Systems durch die mechanischen Ko- 
ordinaten x, sowie die Temperatur # beschrieben werde. Dariiber hinaus sind keine Voraussetzun- 
gen eingefiithrt worden. 

Sobald die Zustandsfunktionen U und S bekannt sind, ist nach (3.7) der Zusammenhang zwi- 
schen den x, und den X{ gegeben. Hinsichtlich der Xj, (bzw. der X;, im vollkommen irreversiblen 
System) folgt indessen aus den beiden Hauptsatzen bloB, daB sie jedenfalls von den x, abhangen, 
und daB die Dissipationsleistung 


= Xix, =0 (5.1) 


ist, wobei das Gleichheitszeichen im allgemeinen nur im Fall der Ruhe gilt. : 
Die Onsagersche Theorie beschrankt sich auf Systeme, in denen die X; nur von den x; ab- 
hangen, und zwar in der linearen Form 


Xi =O, ; x; « (5.2) 
Ihr zufolge ist die Matrix (a,;) symmetrisch, d.h. es gelten die Onsagerschen Relationen 
Aj, = 4%; - (5.3) 
Fiihrt man (5.2) in (5.1) ein, so ergibt sich 
Li = a,; x, x; = D(a) 20. (5.4) 


Hiermit ist die Dissipationsleistung vermittels der Dissipationsfunktion D in den x, aus- 
gedriickt, und (5.2) kann mit Riicksicht auf (5.3) in der Form 
ee a), 

Maa (5.5) 
angeschrieben werden. Da umgekehrt aus dem Vergleich von (5.2) und (5.5) auch (5.3) folgt, sind 
die Beziehungen (5.5) den Onsagerschen Relationen Aquivalent. Man kann sie im euklidischen 
Raum der x, dahin deuten, daf der Kraftvektor! Xj der halbe Gradient der Funktion D(x,) und 
daher in jedem Punkt x, normal zu der ihn enthaltenden Hyperflache D = const sei. Mit der 
quadratischen Form (5.4) ist die Matrix (a,;) im allgemeinen positiv definit. Die Niveauflachen 
von D(x,) kénnen daher als Hyperellipsoide mit dem Ursprung als Mittelpunkt bezeichnet werden. 
Sie sind nach (5.4) ahnlich, und zwar enthalt jede von ihnen diejenigen mit kleineren Werten 
von D. Somit hat der Kraftvektor X} die Richtung der auBeren Normale », im Punkt x,. 

Abb. 1 veranschaulicht diese Zusammenhange in zwei Dimensionen. 


D=konst £, 


Abb. 1. Das Prinzip vom kleinsten Kraftbetrag im zweidimensionalen Geschwindigkeitsraum. 


6. Das Prinzip vom kleinsten Kraftbetrag. Die Beziehungen (5.5) und damit die Onsagerschen 
Relationen kénnen als Ausdruck zweier verschiedener Extremalprinzipien gedeutet werden. Jedes 
derselben stellt in physikalisch durchaus einleuchtender Weise den Begriff der Dissipationsfunktion 
in den Vordergrund. 

Das erste Extremalprinzip lautet wie folgt: 

Nimmt bei gegebener Dissipationsleistung Li und vorgeschriebener Richtung 
des Kraftvektors Xj dessen Betrag unter der Nebenbedingung D(x,)= Li ein 


* Darunter ist hier selbstverstandlich nicht ein Vektor im physikalischen Raum, sondern der Vektor mit 
den Komponenten Xj, im euklidischen Raum der x, zu verstehen. 
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Minimum an, dann bestimmt sich hieraus der Zusammenhang zwischen den 
Vektoren Xj und x. 


Dieses Prinzip wird durch Abb. 1 nahegelegt. Die Angabe von 
= -D(x;,) = XP x0, (6.1) 


schreibt namlich als geometrischen Ort fiir den Punkt %, ein bestimmtes Hyperellipsoid vor und 
fixiert gleichzeitig das Skalarprodukt der Vektoren X} und x,- Der kleinste Betrag des in seiner 
Richtung festgelegten Vektors X} ergibt sich fiir den Vektor xz, dessen Projektion auf die Rich- 
tung von Xj maximal ist, d.h. fiir denjenigen Punkt x, auf dem Hyperellipsoid, dessen auBere 
Normale », die Richtung Xj, hat. 

Fiir die strenge Begriindung sei der Kraftvektor X} durch seinen Betrag B = 0 und den festen 
Kinheitsvektor y, dargestellt. Dann gilt 


Xi a, = By, x, = Li =0, (6.2) 


wobei L‘ der vorgeschriebene Wert der Dissipationsleistung ist, und das Minimum von B wird fiir 
diejenigen Vektoren x, erhalten, welche das Skalarprodukt », x, zum Maximum machen. Dabei 
sind zum Vergleich nur die Vektoren x, zuzulassen, welche der N ebenbedingung 


D(x;,) == (6.3) 
gentigen; ferner wird man sich wegen (6.2) auf denjenigen Teil der Hyperflache (6.3) beschrinken, 


auf dem », x, 20 ist. Auf seinem Rand ist »,*, = 0 und somit das Skalarprodukt minimal. 
Alle weiteren Extrema ergeben sich aus den Bedingungen 


0 : : 
gtd eC a8 (6.4) 
baw. 
aD 
heat pes 0, (6.5) 


in denen A den Lagrangeschen Multiplikator bezeichnet. 
Macht man jetzt von der speziellen Form (5.4) der Dissipationsfunktion Gebrauch, so folgt 
zunachst die Eulersche Beziehung 


iy =2D. (6.6) 


oD 
Oxy, 
Sodann gilt nach (6.2), (6.5) und (6.6) > 
Li = Bry iy = BAS? % = 2 BAD. (6.7) 
Damit ergibt sich 4 mit Riicksicht auf (6.3) aus 
Bi= >, (6.8) 
und man hat zufolge (6.5) 
oe By => 5 (6.9) 


in Ubereinstimmung mit (5.5). ; e,.. ee 
_ Aus der Beziehung (6.9), die man auch in der Form (5.2) schreiben kann, folgt die Hindeutigkeit . 
der Lésung des Extremalproblems (6.4). Da ferner nach (6.5), (6.6) und (6.8) fiir den Lésungs- 


ektor 
‘ D 


5 aD . 
Me ee ig ther oe B= © (6.10) 
gilt, ist das Extremum ein Maximum. Damit ist die Aquivalenz der Onsagerschen Relationen 
mit dem zu Beginn dieses Abschnittes formulierten Minimalprinzip nachgewiesen und iberdies 
gezeigt, daB das Prinzip im vorliegenden Fall genau eine Lésung besitzt. 


7. Das Prinzip von der gréften Dissipationsleistung. Mit dem eben besprochenen ist das folgende 
Extremalprinzip aufs engste verwandt: 

Nimmt bei gegebenem Kraftvektor Xj die Dissipationsleistung L' unter der 
Nebenbedingung Li = D(%,) ein Maximum an, dann bestimmt sich hieraus der 
Zusammenhang zwischen den Vektoren Xj und %. 
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Dieses Prinzip wird durch Abb. 2 nahegelegt. Mit Xj ist némlich die Stellung der zum Ver- 
gleich zugelassenen Hyperebenen 
Bey i (7.1) 
gegeben. Die zugehérigen Hyperellipsoide 

D(x,) = Li (7.2) 


sind ahnlich, und da der Abstand Li der Hyperebenen vom Ursprung mit dem Quadrat der Linear- 
abmessungen der zugehérigen Hyperellipsoide wachst, 1aBt sich die Nebenbedingung nur fiir 
kleinere positive Werte von Li erfiillen. Der Hichstwert ist durch diejenige Hyperebene bestimmt, 
welche das zugehérige Hyperellipsoid gerade noch berihrt, und im Berithrungspunkt hat die 
-auBere Normale », des Hyperellipsoides wieder die Richtung von Xj. 


Abb. 2. Das Prinzip von der gréBten Dissipationsleistung im zweidimensionalen Geschwindigkeitsraum. 


Zur strengen Begriindung sei zunachst festgestellt, daB alle der Nebenbedingung (7.2) geniigen- 
den Punkte auf der Hyperflache 


F(x) = D(x) — Xi xy = 0 (7.3) 
liegen, von der wieder nur das Teilstiick mit Xj x, =O betrachtet werden mufB. Auf seinem 


Rand ist Xj x, —= 0 und daher die Dissipationsleistung minimal. Alle weiteren Extrema ergeben 
sich aus den Bedingungen 


(7) a8 5 


wobei 2 wieder der Lagrangesche Multiplikator ist. Fiihrt man (7.4) mit Hilfe von (7.3) aus, so 
kommt 


(+a Xi—iz =o. (7.5) 


Wenn man jetzt von der speziellen Form (5.4) der Dissipationsfunktion Gebrauch macht, 
so folgt zunachst wieder die Eulersche Beziehung (6.6) sowie aus (7.1), (7.5) und (6.6) 


aD. 2A 


Mags ra Sk = 
SP eee ea a m7?» (7.6) 
mithin nach (7.2) 
1 Slee lee (7.7) 
Damit geht aber (7.5) wieder in 
‘Mee 


d. h. in die Beziehung (5.5) iiber. 
Aus der Relation (7.8), die man auch in der Form (5.2) schreiben kann, folgt die Eindeutigkeit 
der Losung des Extremalproblems (7.4). Da ferner nach (7.1) und (7.2) fiir den Lésungsvektor 


gilt, ist das Extremum ein Maximum. Damit ist die Aquivalenz der Onsagerschen Relationen mit 
dem zu Beginn dieses Abschnittes formulierten Maximalprinzip nachgewiesen und iiberdies ge- 
zeigt, daf das Prinzip im vorliegenden Fall genau eine Lésung besitzt. 
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Pa ae < ee eee ist nur solange sinnvoll, als der Zusammenhang 
gxeiten x;, die lineare Gestalt (5.2) hat. Im Gegen- 
satz hierzu behalten die eben besprochenen Extremalprinzipien auch fiir allgemeinere Zusammen- 
hange zwischen den Kraften und den Geschwindigkeiten ihren Sinn. Es liegt daher nahe, diese 
Prinzipien als Verallgemeinerungen der Onsagerschen Theorie zu betrachten und nach dea oft, 
schrankungen zu fragen, welche sie den mdglichen Zusammenhangen zwischen den Xi und den x 
auferlegen. ; : : 
; Zu diesem Zweck seien Systeme betrachtet, in denen die Xi wiederum blo® von den x, ab- 
hangen. Ferner sei angenommen, daf die Dissipationsfunktion D(x,), die jetzt nicht mehr cue 
dratisch zu sein braucht, aber nach wie vor keine negativen Werte annehmen darf, glatt sei. Da 
man die praktisch zu erwartenden Funktionen D(x;,) wohl mit geniigender Genauigkeit durch 
glatte Funktionen annahern kann, bedeutet diese letzte Annahme keine wesentliche Einschrankung. 
Geht man jetzt vom Prinzip in Abschnitt 6 aus, so kann man die Beziehungen (6.2) bis (6.5) 
tbernehmen. Statt (6.7) hat man unter Beriicksichtigung von (6.3) 


(ee OD. 
ioe = Bh hy — De (8.1) 
mithin 
oD .\71 
Bi= eo i, D (8.2) 
und daher, da nach (6.5) 
: aD 
AL = By, = Bi iy (8.3) 
ist, 
| OD Nee 


Das ist die vom Verfasser in einer friitheren Arbeit! vorgeschlagene und dort durch eine semi- 
statistische Betrachtung begriindete Beziehung, welche auch im allgemeineren Fall die Ortho- 
gonalitat des Vektors Xj, auf der Flache D = const statuiert. 

In ahnlicher Weise kann man vom Prinzip in Abschnitt 7 ausgehen und die Beziehungen (7.1) 
bis (7.5) tibernehmen. Statt (7.6) kommt jetzt mit Riicksicht auf (7.2) 


oD . 


: es A 
Pe mpage (8.5) 
mithin 
A eD.\71 
Eee Or es i) 2 | e9) 


und daher nach (7.5) wiederum (8.4). Die beiden Extremalprinzipien sind also auch in diesem 
allgemeineren Fall aquivalent. Im ibrigen priift man leicht nach, daB (8.4) mit der Dissipations- 


funktion (5.4), d. h. mit (6.6) auf (5.5) zuriickfihrt. 


9. Bemerkungen. Es ist klar, daB die Verallgemeinerung der Prinzipien vom kleinsten Kraft- 
betrag und von der gréBten Dissipationsleistung auf nichtlineare Zusammenhange zwischen den 
Xi und den x,, wie sie im letzten Abschnitt vorgeschlagen worden ist, rein hypothetischen Charak- 
ter hat. Immerhin sind aber die Beziehungen (8.4), wie bereits bemerkt wurde, urspriinglich auf 
ganz anderem Wege gewonnen worden und ruhen damit auf zwei verschiedenen Stiitzen. Ferner 
darf darauf hingewiesen werden, daB insbesondere das Prinzip von Abschnitt 7 vom physikalischen 
Standpunkt aus recht suggestiv erscheint: wenn nach dem zweiten Hauptsatz jedes abgeschlossene 
System dem Zustand maximaler Entropie zustrebt, so erscheint die Annahme nicht unverniinftig, 
daB ein System unter gegebenen Kraften Xj seinen Zustand derart Andere, daB in seinem Inneren 
in jedem Augenblick ein Maximum an Entropie produziert wird. Gerade dies ist aber nach (3.4) 
der Inhalt des Prinzips von der gré8ten Dissipationsleistung. 

SchlieBlich darf wohl iiberhaupt nicht erwartet werden, da® eine Aussage iiber den Ablauf 
irreversibler Prozesse ohne Verwendung gewisser plausibler Annahmen zu formulieren sei, nach- 
dem schon die statistische Begriindung der Onsagerschen Theorie, wie H. B. G. Casimir? nach- 
gewlesen hat, von solchen Zusatzvoraussetzungen Gebrauch machen mu. 


1 H. Ziegler, Z. angew. Math. Phys. 9b (1958) S. 748. 
2 H. B. G. Casimir, a. a. O. 
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Man mag freilich einwenden, daB die Onsagerschen Relationen in besonderen Fallen nicht in 
der einfachen Gestalt (5.3) formuliert werden, sondern gelegentlich in der Form 


. a; ;,(b;) =a, (— b;) » (9.1) 
oder 
(9.2) 


LETS ecier  a hc 
Die Gestalt (9.1) bezieht sich auf Systeme mit gyroskopischen Kraften!, die bekanntlich keine 
Arbeit leisten und daher nicht in (2.1) eingehen. Sie treten zum Beispiel in einem Magnetfeld der 
Feldstarke 6; in Form der Lorentzkrafte oder als Corioliskrafte in einem mit der Winkelgeschwindig- 
keit 6; rotierenden Koordinatensystem auf?. Die Gestalt (9.2) dagegen ergibt sich dann, wenn 
gewisse verallgemeinerte Krafte und Geschwindigkeiten miteinander vertauscht werden*. Wenn 
man aber, wie das hier getan wurde und auf Grund der Analogie mit der analytischen Mechanik 
verniinftig erscheint, die Theorie auf dem fundamentalen Ausdruck (2.1) fiir die Elementararbeit 
aufbaut, der fiir jeden Satz von Lagekoordinaten die zugehérigen verallgemeinerten Krafte defi- 
niert, dann ist eine solche Vertauschung und damit die Form (9.2) der Onsagerschen Beziehungen 
im Vornherein ausgeschlossen. Was ferner die Relationen (9.1) anbelangt, so besteht noch immer 
die Méglichkeit, daB 
ay, (— bj) = ay, j(b;) (9.3) 


_ ist, so daB die Onsagerschen Beziehungen in ihrer normalen Form (5.3) giiltig bleiben. Es ist aber 
auch denkbar, daB die Giltigkeit der Extremalprinzipien von Abschnitt 6 und 7 auf nichtgyro- 
skopische Systeme beschrankt ist. 


(Eingegangen am 19. Oktober 1960.) 
Anschrift des Verfassers: Professor Dr. Hans Ziegler, Riischlikon (Schweiz) Weiherweg 6 


' Vel. H. Ziegler, On the Concept of Elastic Stability, in H. L arma i 
scoot hoe etal wy Po p seh on ability, in . Dryden und Th. v. Kérmén, Advances in 

2 Vel. S. R. de Groot, a. a. O., S. 8. 

3 Vel. S. R. de Groot, a. a. O., S. 216. 
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Der Spannungszustand in einem geschichteten Korper 
bei axialsymmetrischer Belastung 


Von Hans Bufler 


1. Einleitung. In einer kiirzlich erschienenen Arbeit des Verfassers! wurde der ebene Spannungs- 
bzw. Formanderungszustand in einer geschichteten Scheibe bei gegebener Randbelastung oder 
Randverschiebung untersucht. Hier soll das entsprechende axialsymmetrische Problem behandelt 
werden, wobei als Schichten lauter unendlich ausgedehnte, homogene und isotrop-elastische Platten 
mit beliebigen Dicken und beliebigen Elastizitatskonstanten zugrundegelegt werden. Wie frither 
unterscheidet man auch hier den Grenzfall vollkommener Haftung und den Grenzfall verschwinden- 
der Reibung zwischen den Schichten. 

Das ebene Problem konnte mit Hilfe der Methode der Fourier-Transformation in tibersicht- 
licher Weise behandelt werden, wobei sich sogar die Einfiihrung einer Spannungs- oder Verschie- 
bungsfunktion eriibrigte. Die axialsymmetrische Aufgabe wird im folgenden mittels der Hankel- 
Transformation nullter und erster Ordnung durchgefiihrt. Fir die bei der Formulierung der 
Ubergangsbedingungen bendtigten Transformierten der Randspannungen und Randverschie- 
bungen ergeben sich hierbei Ausdriicke, die jenen beim ebenen Problem véllig analog aind. Da- 
durch ist es méglich, einen groBen Teil der bei der geschichteten Scheibe gewonnenen Ergebnisse 
direkt zu titbernehmen. 


2. Mathematische Vorbemerkungen. Die in Zylinderkoordinaten dargestellte Bipotentialglei- 
chung (fiir die Lovesche Verschiebungsfunktion) wird bei Axialsymmetrie mittels der Hankel- 
Transformation nullter Ordnung in eine gewéhnliche Differentialgleichung vierter Ordnung mit 
konstanten Koeffizienten tbergefiihrt, deren Lésung an die an beiden Randern einer unendlich 
ausgedehnten Platte vorgeschriebenen Randbedingungen angepaBt werden kann. Darin liegt 
gerade der Vorteil dieser Methode. 

Im folgenden werden die anschlieBend benétigten Formeln aus der Theorie der Hankel-Trans- 
formation zusammengestellt. 


Wenn das Integral If (x) dx absolut konvergiert, und wenn die Funktion f(x) in der Nachbar- 
0 


schaft des Punktes x von beschrankter Schwankung ist, dann gilt? fiir n => — 1/2 
Flr) = f Gr? Lr) da f (2.x)! f(x) T(2-x) de. 
0 0 
An einer Sprungstelle erhalt man statt f(r) den Wert [f(r + 0) + f@—0)]. In vorstehender 


Formel bedeutet J,, die Besselsche Funktion erster Art und n-ter Ordnung. Ersetzt man f(r) 
durch r1/2 f(r) und f(x) durch «1? f(x), dann folgt 


eee fi Land fx Fail (aya | (1*) 
Dis Hankel-Transformation n-ter Ordnung der Funktion f(x) [oder f(r)] ist definiert durch 
(2) = f x f(s) I,(2.) de = F rf) Ty(er) dr. (2*) 


In Verbindung mit (1*) ergibt sich fiir die inverse Transformierte 


f(r) = faf@) I(r) da. (3*) 


Z : Z. angew. Math. Mech. 41 (1961) S. 158 
2 i Peal ee Tecnbiidasforns. S. 52. New York-Toronto-London 1951. 
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In den in Ziff. 3 vorkommenden Funktionen (Spannungs- und Verschiebungskomponenten) 
wird die Hankel-Transformation nullter und erster Ordnung benétigt. Schreibt man dafiir ein- 


facher f= f und ie : so erhalt man gema4B (2*) und (3*) 
=I rfl) I (Ar) dr bzw. fa) = = af rf(r) L,(Ar) (4*) 


und 


— fafa) (anda baw. f(r) =fa fla) Tar) da. (5*) 


Fir die Transformierten der Funktionen ‘3 (r) = df/dr und 
o? ih 2) 


AON |gat ea a St9) 


gelten folgende Beziehungen!: 


fa) =—af@ (6*) 

(Voraussetzung: r f(r) = 0 fiir r = 0 und r =o), 
Af(i.y) = lest i) f fy) (7*) 
(weitere Voraussetzung: r f(r) = 0 fir r = 0 und r = 00). Zweimalige Anwendung von (7*) liefert 
AAG.) = (sa *) £09). (8*) 


Ersetzt man in (6*) f durch Af, so folgt in Verbindung mit (7*) 
5xa Fe . 
Aft, 9) =—1(52—#) f0.9) (9") 


3. Spannungs- und Verschiebungszustand einer beidseitig axialsymmetrisch belasteten Platte. 
Im folgenden soll der Spannungs- und Verschiebungszustand einer unendlich ausgedehnten Platte 
(Schicht) der Dicke h bei gegebenen (normal und tangential gerichteten, axialsymmetrisch ver- 
teilten) Randlasten bestimmt werden (Abb. 1). 

Radialverschiebung u, Axialverschiebung v, Schubspannung T=t1,y, Axialspannung o,, 
Radialspannung o, und Ringspannung o, lassen sich wie folgt durch die Lovesche § ecekicbang. 
funktion ®, welche der Bipotentialgleichung 


AMP. =(Sat > a + ape) PO9) =0 (1) 
geniigen muB, ausdriicken?: 

aun =—2 mee, (2) 

Tp my) =— SEY + 20-1) A9G,y), (3) 

ars 9) = 5 |— ape + (19) ADE, y)]; (4) 

(6) = |— O69 0) APE, 9) (5) 

o,(t.y) = 5 |G + y VAP = 5 [pe mt oe (6) 
64(09) = 5 |= 1 ale 2”, » AB(r9)). (2) 


Hierin bedeutet E den Elastizitatsmodul und y die Querdehnzahl der Platte. 


Ae S. FuBnote 2 von Seite 417, dort S. 62. 
2S. z. B. K. Marguerre, Z. angew. Math. Mech. 35 (1955) Seite 242 [dort Gl. (29) und (30)]. 
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Wendet man auf (1) die Hankel-Transformati 
es a : at Ilter O 1 ; 
sewshnhiche Tate tcatidlolsiohuihg ae Blane nullter Ordnung an1, so folgt wegen (8*) die 


(D? — #2)? OA, y) = 0, (1) 
wobei D den Operator D = 0/dy bedeutet. Ih 6 : i 
als eee eka | eee Tt AR, ae Or 2 
D(A, y) = A*(A) Cof [Al y + B*(A) Sin Jal y + C*(A) |A| y Cof [Al y + D*(A) Ja] y Sin Jal y. 1”) 
Fihrt man nun mit (2) und (4) die Hankel-Transformati i 
Hub - mation erster Ordnung, mit (3) und (5) die- 
jenige nullter Ordnung aus, so lassen sich die rechten Seiten dieser Gtk ayn Ne a *) 


und (4*), (7*) durch die auf Grund von (1”’) als bekannt anzusehende Funktion (A, y) ausdriicken 
(die Ableitung nach y wird von der Transformation nicht beriihrt) : 


5 Hy) = 2D HEy), | 2" 
Ee - 

Tp 29) = [0 — 2») D2 (1—») 2] BG, 9), (3 

T(t, y) = Af DE + (L—») 2] BO, y), (“’ 

Gy(2, y) = [1 —») D9— 2 —») #2 D] Bl, y) 6" 


Das gleiche Vorgehen fiihrt indessen bei den Gleichungen (6) und (7) nicht zum Ziel. Vielmehr 
hat man dort auf den Ausdruck 0®/dr = © gemaB (6*) die Transformation erster Ordnung, auf 
die Ausdriicke ® und A® gemaB (4*) und (7*) die Transformation nullter Ordnung anzuwenden. 
Man erhalt dann (da die betreffenden Terme sofort wieder riicktransformiert werden miissen) 


a(r.9) =— 7 | BD Gy) her) d+ | AbD + Lv) #D] 86,9) Ler) dh, (6') 
0 0 
ols y) == | 12 D Oi, y) Iy(ar) dk-Ly | 4 [D? — 22 D] BA, y) Iq(4 7) da. (7/ 
0 0 


a ae et ee 


ween Ti (r) 
Ts 


Abb. 1. Axialsymmetrisch belastete Platte. 


Fihrt man statt der zunachst noch willkiirlichen ,,Konstanten“ A*, B*, C*, D* in (1'’) die 
,.Konstanten“ A, B, C, D gemaB 


A* [B+ (lL—2») C], 
Be| 1 JIA + (1 — 2) DI, 
(Os STORET a lea ape 
D* € 


1 Die Anwendung der Hankel-Transformation zur Losung des axialsymmetrischen Problems einer unend- 
lich ausgedehnten Platte geht auf I. N. Snedden zuriick (s. FuBnote 2 der Seite 417). 
2 Da bei der Riicktransformation gema8 (5*) nur iiber 0 < A S on integriert wird, kann man (im Gegen- 
satz zum ebenen Problem) in (1'’) (und damit auch sonst) auf die Absolutstriche bei / verzichten. 
3 Diese sind lediglich in bezug auf y Konstanten. 
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ein, was ohne weiteres zulassig ist, so gehen damit die Ausdriicke (2’) bis (7’) ttber in 


u(A, y) =— a [Ate (1—v) D] oj Ay-+[B+2 (1—v) C]Sinsy+CiryCojay+Diy Sindy}, (2”) 


v(A, y) = ae {{B—(1—21) C]€oj 24y+-[4—(1—2) D] Sinayt+DiyCojay+CiaySinday}, (3”) 


A 
T(A, y) =— [(B + C)Goj Ay + (A+ D) Sindy + Diyoj ay + Cay Smay], (4’’) 
6,(a, y) = AGoj Ay + BSindy + CayOojay+ Diy Sindy, (5) 


o,(r,y) = — | {[A-+2(1—») D]Coj Ay-+[B-+2 (1—») C]Sindy+ Cr yCojay+Diy Sina y} I(r) da 
0 
—[? [(A-+2 D)Coj ay + (B+ 20) Ginay + Cay Coj ry+Day Sindy] Iya 1) da, (6") 


iN) 


o,(t,¥) = —— | {L442 (1—») D]oj2y+[B +2 (1—») C]Sinay + Cay Gof ay + Day Sindy} 


0 
loc) 


x L,(Ar) di— 29 [A @Goiay + CSiniy) Lar) da. (7"") 


0 


Die Formeln (3'’) und (5’’) stimmen vydllig mit jenen des ebenen Problems (ebener Formande- 
rungszustand) iiberein, wenn man dort! die Fourier-Transformierten durch die Hankel-Trans- 
formierten nullter Ordnung ersetzt und die hier iiberfliissigen Absolutstriche bei A weglaBt?. Ebenso 
erhalt man Ubereinstimmung der Ausdriicke (2') und (4) mit den entsprechenden des ebenen 


Falles, wenn man dort® die mit i (i = /— 1) dividierten Fourier-Transformierten durch die Hankel- 
Transformierten erster Ordnung ersetzt und die Absolutstriche bei A streicht?. 

Dieser Befund ist nicht tiberraschend; wie K. Marguerre gezeigt hat+, laBt sich namlich auch 
beim ebenen Problem eine Verschiebungsfunktion @ einfiihren, die der Bipotentialgleichung 
AA® = 0 (hier natiirlich mit A= 0?/dx" + 0?/dy*) geniigt und aus der sich die Verschiebungen 
u(x, ¥), v(x, y) und die Spannungen 1(x, y), o,(x, y) nach den gleichen Vorschriften (2), (3), (4), (5) 
(mit r = x) errechnen lassen. Da auSerdem bei der im ebenen Fall anzuwendenden Fourier-Trans- 
formation die Funktion Af(x, y) sich analog (7*) transformiert, gelangt man dort zu den gleichen 
Ausdriicken (3') und (5’). Bei der Fourier-Transformation von (2) und (4) erhalt man wegen der 


dabei anstelle von (6*) geltenden Beziehung fi (4) =—iaf(’) (= y— 1) Formeln, die sich von 
(2') und (4’) lediglich durch die auf den rechten Gleichungsseiten hinzutretende imaginare Einheit i 
unterscheiden. Ferner fiihrt auch die Fourier-Transformation der Bipotentialgleichung 44D(x, y) 
= 0 auf die Differentialgleichung (1') und die Lésung (1’’), in der allerdings die Absolutstriche 
bei 2 beibehalten werden miissen® und folglich auch in (2’’) bis (6'’) auftreten. : 


Die ,,Konstanten“ A, B, C, D werden durch die Randbedingungen (Abb. 1) 
t(r, 0) =—q'(r), u(r, h) = — q(r), 
~ or, 0) =—p'(r), — o,(r, h) = — p() 


festgelegt. Geht man zu den transformierten Spannungen T und o, tiber, so lauten diese 


7(4,0)=—q'(), TA, h) = —G(), 
Fy(4,0) = —p'(4), Gy (2, h) = — pla) 


1S. meine in FuBnote 1 von Seite 417 zitierte Arbeit, dort Gl. (18) un 4 i iicksichti 

UES a eet ss an (18) und (14) bei Beriicksichtigung von (13). 

* in meiner zitierten Arbeit Gl. (17) und (15) bei Beriicksichtigung von (13). 

4 K. Marguerre, Ing.-Arch. 4 (1933) S. 332. 

® Dieses ist beim ebenen Problem erforderlich, da bei der (komplexen) Fourierschen Riicktransformation 
- tiber — co [A < + o integriert wird. 


¢ 


ae 
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und liefern mit (4’’) und (5) und den Abkirzungen ©@= Gin ah, C= Coj ah 


A=— p'(A), 
B=q7()—C, 
Ca TW Ght SOPM+GhC+ ©) p™+e7 QE): 
St — (hp 
Da SP A—Ah SpA + Gh—S O)GF'(4) —AhE— S) GA) 
G?— (hy? 


Damit erhalt man fiir die Transformierten der Spannungs- und Verschiebungskomponenten 


gemaB (2’’) bis a 


ty) = aay {[ <(1 — 2») ©? + (Ah)*> p’ + 2 (1—)Ah®@ p 
—2(1—v) @4h—SO)q’ + 2 (1 —») ALC — GS) q] Coj ay 
+ [(1—2) (Ah + ©©) p’ —(1— 27) AAC+ S) p 
+ <(A h)? — 2 (1 —v) Gq’ — (l1—2v) Ah GS Gq] Sindy 
+ [Ah + SG) p'— (Ah + ©) p— Sq’ —AhGS GJ Ay Gof ray 
+ [—G p’+AhS p—(Ah—GG) q’+ AhC—G) q]aySinay}, (2’”) 


Ts 9) = aap (2 C0) G+ SG) p’ 21») GE +S) p 
— (1 — 2 v) G2 + (Ah) 7! —2 (1—v) ANS GT] Cop ry 
+ [<(4 h)? — 2 (lL—v) S* p' + '(l— 27) Ah p 
—(1l— 27) (Ah— SC) q’ + (lL— 27) AAC — GS) q] Ginsay 
+ [G2 p’ —Ah@ p+ Ah— GSC) 7’ — (AhC—G) qT] AyGojay 
+ [—Gh+ SG) p’'+QAhC+S) p+ Sq’ +ahSqlay Sindy}, (3'") 
T(A y¥) = — enum {[@h)? p’—AhS p+ Ah—S C) q’—(AhO—S) q] Sind y + [S°—(h)"] q' Co} ry 
+[—G@h+ SG) p’'+ @hC +S) p+ Sq +¢AhSqjay Sindy 
+ [G2 p’ —AhG p+ Ah— SC) 7 —AhC—G) q]AyOojay}, (4) 
6,(4,¥) = — ea {[S2—(Ah)?] p’ Coj Ay—[(Ah+S C) p’—(AhC+G) p—CAh) q’-AhS q] Sindy 
+ [Ah+ SG) p’—Vh€ + ©) p—SY —AhS q]sayCojry 
+ [-S*p’ +AhS p—Gh— SC) 7 + @hC—GS) q]ry Sinay}, 6") 


ate —+ | wa 9] Ian 


a [ amy {[<SG2+ (Ah) p’—2AhS p+2(Ah—SG) q’—2 AAC—G) q] Coj ay 

é + [—(h+S©) p’'+ (hE +S) p+<2S— (Ah) q’+4AhS q] Sindy 
+ [—(@h+ SG) p’'+ @hE+ S) p+Srg’ +¢Ah Sq] ryGojay 

+ [S* p'—AhS p+ (Ah—SG) q’—(AhC—S) 4G] Ay Sindy} [y(ar) da, (6) 


og(tsy) == [ a ae y) I,(a1) da 
0 
29 {ea gays (C8 Ab S P+ Gh-SOT —GMC—S) Gq] Cojry 
3 + [-(Ah+©@6) p’+(Ah€+©S) p+-S2q’+AhS q] Ginday}Ip(ar) da. (7) 
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Die transformierten Randspannungen sind gemaB (4*) zu berechnen: 


ee 1 I,(A r) dr , 
PO! J leo i 


aE Bet L(ar) dr. 
q@ J We) : 


Aus den transformierten Verschiebungen und Spannungen (2'”’) bis (5'”’) erhalt man die wirklichen 
mittels Riicktransformation entsprechend (5*): 


u(r, y) Lf ea, I,(A 7) a, | 
o,(r, y) : oy(A, y) 


wry z | Fi GE ei ‘| L,(r) di. 
u(r, y) ; t(A, y) 


(8) 


(9) 


Abb. 2. Der axialsymmetrisch belastete Halbraum. 


Im besonderen gelten fiir die Spannungen im Halbraum die aus (4’’’) bis (7'") mit y = h—y7 
und h-—> oo (Abb. 2) folgenden Formeln: 


u(r, 0) = as A[—An p(t) + (L—Am) q()] e747 LAr) da, 
o,(7,4) = — if A[(L + 4m) pQ) +4 q(A)] e7" IA 7) a2, 


a(n) = { {(L—2»—An) PQ) + 2 Lv) An] FW} "Hr de 
i (10) 
— [ Al —An) p@) + @—Am) FO] e-*” LAr) da, 
agro) =— f {1 —2» 49) p@) + [2 =») — Am] GO} em" Lr) 


co 


—=2y { A [p(A) + (A)] e-47 I,(ar) da. 


0 


Liegt im Intervall 0 <r Sa eine gleichmafig verteilte Druckverteilung p, vor, so bekommt man 
mit dem Integral 


| rI(ar) dr= + I,(A a) 
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fiir die transformierte Belastung 


a 


Pol) = Po | r I,(A r) dr = am py Jy(Aa) 


Mt ha 
6 


Geht hlieBli i ; : . “5 
win a * a ae a—>0 zu einer Einzelkraft (oF nae (a2 7 pp) iiber, dann ergibt sich 
7a, a any 


wegen lim 
2 Po? 


a0 


Po(4) = 5 (11) 


und damit, wenn man auferdem von einer Randschub bsieht [q(A) = 
eee eae ubspannung absieht [q(A) = 0], das wohl- 


3 2 = 
U(r, 1) = 5— Py “Ee (R= yr? +77), 
3 3 
o,(r, 4) = — 5— Pye 
12) 
1 2 1—2y t ( 
6,(r, 4) = — == Ps 3 a 5 : ( a, ) el 
l1—2» 1 n n | 
Oo(T. 4) = Da Pal a(1— 8) we J 
Hierbei wurden die Integrale 
Ae—*n Tar) dA = 2 2 e—An TA) da = 2 
0 —a R3 ? e of r) re R ? 
é i 
[rem L@nda=, [Benen a3, 
6 6 


| e144 L..4r) dk= ica 2) 
6 


verwendet. 


_ 4. Der geschichtete Kérper im Falle vollkommener Haftung. Die im folgenden anzustellenden 
Uberlegungen sind jenen beim ebenen Problem vdllig analog®. Im Falle vollkommener Haftung 
aller Schichten miissen die Berithrbedingung 

v,(T, 0) = Vp 41(T, hy +1) 
und die Haftbedingung 


Uy(7, 0) = Up a(ts Ay-1) 


erfillt werden (Abb. 3 und 4), wobei der Index k zur Numerierung der einzelnen Schichten dient. 
Multipliziert man die erste der vorstehenden Gleichungen mit r I,(4 r), die zweite mit r I,(4 r) und 
integriert zwischen r = 0 und r = 00, so erhalt man wegen (4*) 


0, (A, 0) = v_41(A, hy +1) (Berithrbedingung) , (13) 
U,(A, 0) = Un. 1(Ay x41) (Haftbedingung) . (14) 


1 §. Biezeno-Grammel, Technische Dynamik 1. Bd. 2. Auflage, S. 118, Berlin-Gottingen-Heidelberg 1953; 
2S, FuBnote 1 von Seite 417. 
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Nach (2’) und (3/") ist aber mit den jeweiligen Randspannungen Pj, x2 Pk—1> Yk—1 bzw. Pris 
Gk+1> Pho Ix (Abb. 4)" 


4,0, 0) =? —"®) (4 Be + 54 Pri — ab Ge + BE Ge) » (15) 
pe ac emmy 7 
2, 1— y? = = = = 16 
v,(A, 0) — (%, Pr — Pe Pr—1 — Vk dx — Ox Ik—1) » (16) 
= 2(1— ) — = % = ers 17 
Up 1(A, hy 1) = oe (On41 PAY yk Pr— Pi+1 qk+1 + Oe qk) 2 ( ) 
2 (1 — v3) = e Ee al E3 18 
Dp4 (A, Api) = a (Oe Pid — Obs 1 Pe On+1 Uk+1 — Yet qk) (18) 
Por %p 
[Poh 
| P24 
Soe eens 
: ; | 
| 
\ Re 
| y, oe ke : Past, Pk, Pr-1| 
| P91 The Gk+1 Tk. Then: 


Py-2) W-2 
Py-1 IN-1 


Py Wy 
Abb. 4. Zur Herleitung der Beriihr- und Haftbedingung. 


Abb. 3. Der aus N Schichten bestehende Kérper. 


mit den Abkirzungen 
oo, S Aly, + Sin A hy, Co} 2h, , 
art = A hy, — Gina hy, Coj A hy, Py 
= A hyo) Ay, + Sindh, , 
Bt =1h,Coj Ah,— Sindh, , (ag) 
6, = Ah, Sina h, , 
ee = Sir? A h,, = (A h,)? 3 
_ = 2%) Sin? A hy + (A hy)? 
Ve ye aay) ae 


Damit lauten jetzt die Beriihr- und Haftbedingung 
Bist Pei — (Oe 41 + Oe Pk, k+1) Ph + Be Pk, k+1 Pk—1 = 
= Sk+1 Gert — (Yer1— Ve Ph, +1) Ue + Oe Pkse+1G—1 = 0 (Berithrbedingung), (20) 


Ok+1 Pe+1— (Yk+1— Ve Pkyk+1) Pk— Oe PE, k-+1 Pk—1 J . 
— Bier Gera + (1+ Ff Ge, 441) Ie — BE Pk, k+1 Gk—-1 = 0 (Haftbedingung) , (21) 


wobei 
uit ie : 
Maye, 


STS pele 99). 
Lf 5 E, ey ( ) 


Vk, k-+1 = 


bedeutet. 


* Diese Gleichungen gehen nach dem in Ziff. 3 Gesagten aus jenen des ebenen Problems (ebener Form- 


anderungszustand) dadurch hervor, da8 man dort die Fourier-Transformierten p und ov durch die Hankel- 


Transformierten nullter Ordnung, dagegen die Fourier-Transformierten g und w durch die mit i= /—1 
multiplizierten Hankel-Transformierten erster Ordnung ersetzt. 
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Die gekoppelten Differenzengleichungen (20), (21) kénnen bei insgesamt N Schichten je (N — 1) 
mal angeschrieben werden (namlich fir k= 1,2,..., N— 1) und liefern somit die zur Bestimmung 
der Unbekannten p,, py,..., Py—13 > Go+++>4N—1 erforderlichen 2 (N— 1) Gleichungen. In 
diesem linearen Gleichungssystem sind die Po Wo und Py, Iv bei gegebenen Randlasten als bekannt 
anzusehen. Sind statt der Randspannungen die Randverschiebungen vorgegeben (worauf hier 
nicht weiter eingegangen sei), so ist ebenso zu verfahren wie im ebenen Fall?. Nach Ermittlung 
der p;, und q, kénnen der Spannungs- und Formanderungszustand in jeder Schicht gemaf (2’”’) 
bis (7'’") bestimmt werden. AnschlieBend werden einige Spezialfalle betrachtet. 

hg (ee 


3” 
1 ¥22 
Ox 


| a bE RE ne ed 
tie 


® 


Abb. 5. Der Halbraum mit anschlieBender Schicht bei normaler Belastung. 


a) Der homogene Halbraum. Der Halbraum bestehe aus lauter gleichen, aneinander- 
geklebten Schichten. Die Lésung der beiden Differenzengleichungen (20), (21), die in diesem 
Sonderfall konstante Koeffizienten erhalten (a, = a, B, =P, 6, = 0, ye = y, 0 = a*, PE = B*, 
Pr,k+1 = 1), erfolgt genauso wie im ebenen Fall. Das Ergebnis lautet 


Pel?) = [1 + Ahk) pot) + Ahk gale ***, | 


7 = = Ahk () 
Gu(4) = [—A bk pold) + (LAL) Go(t)] o—***.f 

Setzt man hierin p, = — 6,, q, = —T und hk = y, so bekommt man — wie es auch sein muB — 
nach dem Riicktransformieren Ubereinstimmung mit den ersten beiden Formeln (10). 


b) Der Halbraum mit aufgeklebter Schicht. Zugrundegelegt sei hier der in Abb. 5 
dargestellte Verbundkérper bei normaler, axialsymmetrischer Randbelastung. Mit h,—>0o erhilt 
man aus (19) fiir den Halbraum 


euths 
Fes os = gas? 
ethe 
Bo = Be SORE Spas dane 
2855 ab 
PIG 8 (1 — »2) = ees 
und aus (22) 
_ 1+} E,e’*" 
ely Bb We cs 


Hiermit lauten die Berihr- und Haftbedingung (20), (21) fir k = 1 und q = 0, wenn man den > 
Grenziibergang h,—> 0o durchfiihrt, 


1—vj Ey a) - 1—2% 1— Ee v\= _1—%1 FA - 
Poehler et ee 
2 (1 — ») Le E, &y 


= 1—yv E = 1—vr? E, a¥\ = 1—v? E, 6, — 


F 1S. FuBnote 1 von Seite 417. 
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Die Auflésung dieser beiden Gleichungen nach p, und q, fuhrt zu 


— 2 1—2y»y, 1—v E 
Dy F 2 et ah at) | = He hg —1 , 
= Pra WES Sry Mc ae Ee Se D(A) » (24 
Pi Tea E 1— 2 eae | 1—2y, 1—22 E, , Lay 2 Pol) 

(T= 1 2 (1 —y7,)- 1l— 2? Ey 2 


pial S2ap, Eee 3 1—v} E, if 
A latieo 1—»? E, * V1 ar ob a Ea 1 


rf 1—2y, 1—»3 E, 


THN : eee ee ee s 
n() Reds See yeas. ee ee 
l= By 4) te FB, 2(1—») 1—»? E, 


a Pol4)+ (29) 
»| 


Bei gleichem Material von Halbraum und Schicht folgt hieraus 


pilA) = (1 -+ Ah) e~*" pA), 
H(A) = —Ahe—** pA) 


in Ubereinstimmung mit (23), falls man auch dort q, = 0 und k = 1 setzt. 
Hat man es dagegen mit einer Schicht auf starrer Unterlage zu tun, so folgt mit Lk, > 00 
aus (24) und (25) 
= _ bar+on - 


Pil) ay a* aS v2 PolA) ° 
= _ bin tore = 
q(A) = ~aat $y Pol) 


und mit den entsprechenden Ausdriicken (19) 
2 2(1—»,) CofAh, +Ah, Sindh, = — 
px(A) = 2 (1 —) (3 — 4 »,) CoPAh, + Ahye2+ G— 2) Pol) » 


= (1— 2») Gindh, —Ah, Cojrh eas 
qu(A) = 2 (1 —%) (—45,) Cof2 Ah, Ta hy? " a aes m4) Pol) - 


Besteht die Belastung aus einer Einzelkraft P, bei r = 0, dann bekommt man mit (11) nach Riick- 
transformation 


Pe eenene P, [2 2(1—») CofAh, + Ah, Sindh, Ian) da, - (26) 
0 


7 (3 — 4») CoPAh, + Ah)? + — 2%)? 


q(r) = Lar) da. (27) 


1— p [4 G—2) Sindh, —Ahy CofA hy 
nm : (3 — 4 ry) Cof?A hy + (Ay)? + (1 — 2 4)? 


0 


5. Der geschichtete Kérper im Falle verschwindender Reibung. Hier mu (natiirlich unter der 
Voraussetzung, da®B sich die Schichten nicht voneinander abheben) die Berihrbedingung (13) 
erfillt werden, wahrend die Haftbedingung (14) entfallt. Wegen der Reibungsfreiheit ist 


Ch i * =I —0. 
Beriicksichtigt man dieses in (20), so lautet die Beriihrbedingung hier 
Bist Peri — (He+1 + Oe Pe, b+1) Pk + Bi Ph, k+1 Pr—1 = 0. (28) 


Die Differenzengleichung (28), angeschrieben fiir k = 1,2... N—1, liefert N —1 lineare Glei- 
chungen fiir die N— 1 Unbekannten p,, po,..., Pn—1 (Pp und py sind die bekannten Randlasten), 
nach deren Bestimmung der Spannungs- und Verformungszustand [Riicktransformation der 
Gl. (2'’) bis (7'")] der einzelnen Schichten bekannt ist. 


a) Der homogen-geschichtete Kérper. Darunter soll ein Kérper verstanden werden, 
der aus lauter gleichen Schichten aufgebaut ist. Mit a, =a, B, =p, & = & Pr,k+1 = 1 erhalt 
die Differenzengleichung (28) konstante Koeffizienten: 


= Serge 
Pett =e Pu+ Pri = 9. (29) 
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Bei insgesamt N Schichten lautet die an die Randbedingungen (p, und py seien vorgegeben 
gepaBte Lésung! bei Beriicksichtigung von (19) 


) an- 


Ee 1 di a = ee 
PA) = qh gan pw — YP Po) vt + CY" Po— Pr) Y—*] (30) 
wobei 


Ah + SinsAh CojAh + VSinwtAh —(Ahy Sindrh 3] 
? AhGojfAh+ Ginth (3p) 


(Ah)? 
23 
Fir den homogen-geschichteten Halbraum folgt aus (29) mit N->0o 


Pe(A) = Pol) p—*(A) 


und firAh<lywl+ 


ist. 


und nach Ricktransformation 
PiAr) =i) A polA) p—*(A) [y(Ar) da. (32) 


Als Beispiel soll im folgenden der Spannungszustand in dem durch eine Einzelkraft Py am 
Rand beanspruchten homogen-geschichteten Halbraum untersucht werden. An der unteren Be- 
grenzung der kten Schicht tritt nach (32) und (11) der Druck 


co 


Pale) = arpa | Mt wt Let 1) da (33) 


0 


r 


auf. Hierbei bedeuten A* und r* die dimensionslosen GréBen A* = Ah und r* = i 


sowie gemaB (31) 


_A* + Gin dA* Cof A* + VSineA* —1* Sin d* (34) 
? 1* Coja* + Sind* 
Die Spannungen o, und o, am unteren Rand (y = 0) der kten Schicht betragen nach (6’”) und (7'”) 
mit q = q’ = 0 sowie p = p,_; und p’ = p, 


o,,(r*, 0) = =a | eae | [(1 2) Gin2a*+ 4*2] p,+ 2 (1—v) a* Sin a* py_a} Ly(4*r*) da* 
0 
— a | eiaege a prallSint at 42%) p22 Sin a pe | TT) a, 85) 
0 


r 1 : = ; ~ 
ogtlt™s 0) = segs | Siege ava {1G —2») Sind a+ 29%] By-+ 2 (Lv) a" Si 44 pia} HO*r*) die 
0 


co 


pase rf Soe: (Gin? a* p, —A* Sin A* py_a) A* Ig(A* r*) da*. (36) 


Rel SintA* 
0 


Fur r* = 0 folgt hieraus 


0,4(0, 0) = oy4(0, 0) = TB | Sa pe jae (2(1+9) a* Sin a pe_s[(1 4-2) Giwea*+4"*] py} ae ae) 
0 


. 4 = P. 
Im betrachteten Beispiel ergibt sich mit p, = oe yt 


| Py fat (Lv) A* y Gina*— [1 + 29) Sina" +A") ag 
6,1(0; 0) = Gpi(0, 0) = Gap | Sinea* — 1"? ph di*. (38) 
| 0 


ai Siehe meine in FuBnote 1 von Seite 417 zitierte Arbeit, dort Gl. (57). 
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Das Integral (33) (mit r* = 0) und das Integral (38), das auBerdem von der ee y 
abhangt, wurden fiir k = 1, 2,3... 10 numerisch ausgewertet und die Ergebnisse in Abb. 6 nie i 
gelegt. Die Spannung o,,(0, 0) wurde dabei jeweils fiir die zwei Extremwerte der Querdehnza 
y = 0 und y = 0,5 ermittelt. In Abb. 6 sind auSerdem die entsprechenden, im nicht-geschichteten 


06 ar 


| B 


Abb. 6. Verlauf der Spannungen P;(0) und 04,0» 0) in dem durch eine normal gerichtete Einzelkraft 


beanspruchten. homogen-geschichteten Halbraum lings der Wirkungslinie der Kraft fiir y= 0 und » = 0,5 
(Giltig fir k = 1, 2,3,.. 3 cecce Verlauf der entsprechenden Spannungen im nichtgeschichteten Halbraum. 


05 70 45 20 


Abb. 7. Verlauf des Drucks p,(r*) in dem durch eine normal gerichtete Einzelkraft beanspruchten homogen- 
geschichteten Halbraum lings des unteren Randes der ersten Schicht; - - - - - Verlauf des entsprechenden Drucks im 
nicht-geschichteten Halbraum, 


Halbraum langs der Wirkungslinie der Kraft P, auftretenden Spannungen p(0) und (o,),-9 = (5,)r=0 
eingetragen. Diese errechnen sich gemaéB (12) zu 


3 1 1—2 1 
P00) =— 00,9) = 3, Pores 9x0.) = 04(0,) = Poa =kh). 


Der Vergleich zeigt, daf der Druck p,(0) im reibungslos geschichteten Kérper wesentlich langsamer 
abklingt als im nichtgeschichteten Kérper; derselbe Befund ergab sich auch beim ebenen Problem. 
Ferner wurde die Druckverteilung unter der ersten Schicht (k = 1) p,(r*) nach (33) ermittelt 
und in Abb.7 kurvenmaBig dargestellt. Zum Vergleich ist dort wiederum die entsprechende 
Druckverteilung im nicht-geschichteten Kérper, fiir welche nach (12) 
3 P. 1 
Pilt*) = — (6y)n=n = on 3 (i + r*2)52 


gilt, eingetragen. 
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b) Der Halbraum mit anschlieBender Schicht. Gema8 (1 i 
ee chic emaf (19) und (22) folgt mit h,—> oo 


ae ert hs | 
= =f 
Z 4 hy 0o 
Ah, 
ele 
By 2 9 Reg 
l—v? EB, 24h 
Loo panes) 
¥ 7 Y Ei fe, 4 €, ee 


womit die Berithrbedinguag (28) (fir k = 1) tibergeht in 


fe 
0 G2 03 OF O5 06 07080970 2 S46 6 7EID 20. 30 40 06080 


Abb. 8. Spannungen p,(0) und o,,(0,0) in der durch eine Einzelkraft normal belasteten, auf einem Halbraum 
reibungsfrei liegenden Schicht in Abhangigkeit vom Parameter f fiir », = 0 und », = 0,5. 


Einsetzen von a,, /,, €, [nach (19)] fihrt zur Formel 
Ah, GojaAh, + Sindh = 
: Or (39) 
2 [Sin A hy — A) + ah, + GinAh, Cof Ah, 


p,(A 
P(A) = ae ae 


1— vy 
Besteht die Belastung aus einer Hinzelkraft P), so erhalt man aus (11) py) = P)/2 2 und mit 
den dimensionslosen Grié®en A* = A h,, r* = r/h, sowie mit der Abkirzung 
tae i — 7 E, 
eae 


aus (39) nach Riicktransformation 


ee A* CofA* + Gin A* © La%r*) dat 
Pe 2a he | f (Sin? A* —A*) + A* + Sin A* CojA* Es Ee i 
0 


Fiir die tangentiale Randspannung am unteren Rand der Schicht bei r = 0 gilt (37) (mit k = 1). 
Setzt man dort p, und p, gemaf obigem ein, dann folgt 


67,(0, 0) = o,(0, 9) 


P at [(L +2») Ginta* LA] (A* CofA" + Sind") og | 
2 i Spree (1b 9)" Sit ee ta | a a) 
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Die Integrale (40) und (41) wurden numerisch ausgewertet. Abb. 8 zeigt den nach (40) errech- 
neten Maximaldruck unter der Schicht p,(0) und die Spannung o,,(0, 0) = 05,(0, 0) nach (41) in 
Abhangigkeit vom Parameter f. Letztere wurde fiir die beiden Extremwerte der Querdehnzahl 
der Schicht (7, = 0 und 0,5) ermittelt. In Abb. 9 ist schlieBlich die aus (40) fiir f=1(d. h. Schicht 
und Halbraum bestehen aus gleichem Material) erhaltene Druckverteilung p,(r*) dargestellt 
und auBerdem wieder die entsprechende fiir den Halbraum ohne Schnitt giiltige Druckverteilung 


eingetragen!. 


GS 


G3 


QZ 


0 


Abb. 9. Verlauf des Drucks p,(r*) langs der unteren Berandung der durch eine normal gerichtete 
Einzelkraft belasteten, auf einem Halbraum reibungsfreiliegenden Schicht fiir f=1 ; - - - - - Verlauf des entsprechenden 
Drucks im Halbraum ohne Schnitt. 


6. Zusammenfassung. Es wird der Spannungszustand eines geschichteten Kérpers bei axial- 
symmetrischer Belastung fiir die Grenzfalle vollkommener Haftung und verschwindender Rei- 
bung mit Hilfe der Hankel-Transformation allgemein bestimmt. Als Schichten dienen dabei 
unendlich ausgedehnte, homogene und isotrop-elastische Platten mit beliebigen Dicken und be- 
liebigen Elastizitatskonstanten. Fiir den homogen-geschichteten Halbraum und den Halbraum 
mit anschlieBender Schicht werden geschlossene Formeln hergeleitet und numerische Ergebnisse 
mitgeteilt. 

Beim Vorhandensein lauter voneinander verschiedener Schichten wird die Rechnung natur- 
gem4B sehr umfangreich; in diesem Fall 1a8t sich das Problem mittels Ubertragungsmatrizen in 
ubersichtlicher Weise behandeln, wie an anderer Stelle noch gezeigt werden soll. 


(Eingegangen am 21. Oktober 1960.) 


Anschrift des Verfassers: Privatdozent Dr. H. Bufler, Miinchen 2, Arcisstr. 21, Lehrstuhl fiir Mechanik. 


ext Vel. hierzu K. Marguerre, Ing.-Arch. 4 (1933) S. 332. 
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Nichtlineare Tragflachentheorie 
insbesondere fiir Tragfliigel mit kleinem Seitenverhaltnis* 


Von K. Gersten 


1. Einleitung. Mit der Entwicklung der Flugzeuge zu immer héheren F luggeschwindigkeiten 
war zwangslaufig auch eine Entwicklung zu neuen, den jeweiligen Erfordernissen angepaBten 
Tragfligelformen verbunden. Verfolgt man die verschiedenen Tragfligelformen von dem ein- 
fachen ungepfeilten Fliigel mit gro8em Seitenverhaltnis und ziemlich dickem Profil der ersten Flug- 
zeuge bis hin zu den schlanken, extrem diinnen Deltafliigeln moderner Uberschallflugzeuge, so 
14Bt sich die Entwicklung im wesentlichen durch drei Merkmale charakterisieren?: 


a) Das relative Dickenverhaltnis der Profile des Fliigels wird kleiner. 
b) Die Pfeilung der F ligel nimmt zu. 
c) Das Seitenverhaltnis nimmt ab. 


Daher ist verstandlich, daB die moderne Tragfliigelaerodynamik sich in verstarktem MaBe mit 
diinnen, stark gepfeilten Tragfliigeln von kleinem Seitenverhaltnis beschaftigt. 

Die Strémung um einen schlanken Deltafliigel unterscheidet sich ganz grundlegend von der 
Umstrémung eines klassischen ungepfeilten Tragfliigels mit groBem Seitenverhaltnis. Der wesent- 
liche Unterschied besteht darin, daB im Gegensatz zum Fliigel mit groBem Seitenverhaltnis beim 
Fligel mit kleinem Seitenverhdltnis der Auftrieb und das Kippmoment nicht mehr linear mit dem 
Anstellwinkel anwachsen. Die Abhangigkeit des Auftriebs und des Kippmomentes vom Anstell- 
winkel 1a8t sich aufteilen in einen linearen Anteil und in einen nichtlinearen Anteil. Letzterer 
gewinnt immer mehr an Bedeutung, je kleiner das Seitenverhaltnis wird, und ist bei Seitenver- 
haltnissen von etwa A = b?/F = 1 bei maBigen Anstellwinkeln (a = 15° bis 20°) schon von gleicher 
GréBenordnung wie der lineare Anteil und daher nicht mehr vernachlassigbar. 

Die vorliegende Arbeit beschaftigt sich mit der theoretischen Berechnung der Auftriebsver- 
teilung und der aerodynamischen Kraftbeiwerte von Tragfliigeln mit besonderer Betonung der 
nichtlinearen Effekte bei Fliigeln mit kleinem Seitenverhaltnis. 


2. Die Strémung um Tragfliigel mit kleinen Seitenverhaltnissen. a) Wirbelschichten am 
Tragfligel. Die Ursache fiir den nichtlinearen Zusammenhang zwischen den Kraftbeiwerten 
und dem Anstellwinkel liegt in der besonderen Art der Strémung um Fliigel mit kleinen Seiten- 
verhaltnissen. Da die Seitenkanten bzw. die gepfeilten Vorderkanten groB sind im Vergleich zur 
Spannweite, wird das gesamte Strémungsbild von der Umstrémung der Seitenkanten bzw. Vorder- 
kanten bestimmt. Wie die Skizzen in Abb. la zeigen, fiihrt die Umstrémung der Seiten- bzw. 
Vorderkanten im allgemeinen zur Ablésung und dadurch zur Bildung von Wirbelschichten, die 
sich sehr bald zu stark konzentrierten Wirbelzépfen aufrollen. Wahrend bei sehr kleinen Anstell- 
winkeln (lineare Theorie) eine Wirbelschicht nach Abb. 1b nur von der Fliigelhinterkante abgeht, 
kommen bei gréBeren Anstellwinkeln die an den Seiten- bzw. Vorderkanten entstehenden Wirbel- 
schichten dazu, die das Bild grundlegend andern. Insbesondere liegen die an den Seiten entstehen- 
den Wirbelschichten und die dazugehérigen Wirbelzépfe nicht mehr in derselben Ebene mit der 
Wirbelschicht, die von der Hinterkante abgeht. Vielmehr andert sich der Abstand zwischen den 
Wirbelzépfen und der Fliigelebene noch mit dem Anstellwinkel, wodurch die nichtlineare Ab- 
hangigkeit der StrémungsgréBen vom Anstellwinkel hervorgerufen wird. Das Auftreten einer 


* Aus dem Institut fiir Strémungsmechanik der Technischen Hochschule Braunschweig (Leiter: Prof. 
Dr. H. Schlichting). Gekiirzte Fassung der von der Fakultat fiir Maschinenwesen der Technischen Hochschule 
Braunschweig genehmigten Habilitationsschrift, 1960 (Referenten: Prof. Dr. H. Schlichting, Prof. Dr. H. Blenk, 
Prof. Dr.-Ing. H. Schaefer, Prof. Dr.-Ing. E. Truckenbrodt). — Ausfiihrliche Fassung: Nichtlineare Trag- 
flachentheorie, insbesondere fiir Tragfliigel mit kleinem Seitenverhaltnis. Bericht 59/30 d. Inst. f. Stromungs- 
mechanik der TH Braunschweig (1959); auszugsweise vorgetragen auf der GAMM-Tagung, Saarbriicken, 
8.—12. 4. 1958, auf dem 10. Internationalen Mechanik-Kongref, Stresa, 31. 8.—7. 9. 1960, und auf der 
Tagung des ,,Fluid Dynamics Panel“ der AGARD in Briissel, 10. 4.—14. 4. 1961. Teilergebnisse findet man 
in K. Gersten, Z. Flugwiss. 5 (1957), S. 276 und K. Gersten, Jahrbuch 1958 der Wiss. Ges. f. Luftfahrt, S. 25. 

1 H, Schlichting u. E. Truckenbrodt, Aerodynamik des Flugzeuges, Bd. I (1959), Bd. II (1960), Berlin/ 
Gottingen/Heidelberg. 
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Nichtlinearitat folgt demnach aus der Tatsache, daB die vom Tragfligel ab- 
gehenden freien Wirbel nicht mehr samtlich in einer Ebene liegen. 


b) Einflu8 der Vorderkantenform. Wie aus Messungen! 2 3 bekannt ist, hangt 
die GréBe des nichtlinearen Anteiles sehr stark von der Form der umstrémten Vorderkanten 


Abb. 1. Entstehung der Wirbelflichen bei Rechteck- und Deltafliigeln. — a) grofer Anstellwinkel (nichtlineare Theorie); 
b) kleiner Anstellwinkel (lineare Theorie). 


bzw. Seitenkanten ab. Da die nichtlinearen Effekte hauptsachlich durch die Ablésung der Stré- 
mung und die Wirbelbildung an den Fliigelkanten hervorgerufen werden, ist verstandlich, daf 
die nichtlinearen Anteile bei zugescharften Kanten gréBer sind als bei abgerundeten Kanten. 
Abb. 2 zeigt fiir zwei Rechteckfliigel mit verschiedenen Formen der Vorder- und Seitenkanten 
die Abhangigkeit des Auftriebsbeiwertes vom Anstellwinkel nach Messungen von H. Winter?. 
Zum Vergleich sind die theoretischen Kurven nach der spater noch zu behandelnden linearen 
Theorie eingezeichnet. Die beiden MeBkurven fiir runde Vorderkante, aber verschiedene Formen 
der Seitenkante, liegen jeweils dicht zusammen. Durch Zuspitzung der Vorderkante dagegen 
werden die nichtlinearen Anteile ungefahr verdoppelt. Das riihrt daher, daB infolge der Zuscharfung 
die Strémung an der Vorderkante zwar ablést, jedoch etwas stromabwarts wieder anliegt, so daB 
die Umstrémung in Tiefenrichtung im grofen und ganzen mit der reibungslosen Strémung um 
ein Profil iitbereinstimmt. Das Totwasser in dem Ablésungsgebiet nahe der Vorderkante wird nun 
bei Tragfliigeln mit kleinem Seitenverhaltnis in das Unterdruckgebiet der Randwirbel hinein- 
gesaugt, was zu einer starken VergréBerung der seitlichen Ablésungszonen und damit der Wirbel- 
schichtbildung fiihrt. Daraus ergibt sich dann die Zunahme der nichtlinearen Effekte. Die Form 
der Seitenkantenausbildung ist deshalb nicht von so groBer Bedeutung, da es bei der seitlichen 
Umstrémung sowieso zu einer Ablésung kommt, unabhangig davon, ob die Kante scharf oder 
rund ausgebildet ist. 

Einen entsprechenden Vergleich fiir die Beiwerte von Auftrieb und Kippmoment bei einem 
Deltafliigel (Seitenverhaltnis A = b?/F =: 0,78, Zuspitzung 4 = 1,/1; = 0,125) und einem Pfeil- 
fligel jeweils mit zugescharfter und runder Vorderkante zeigt Abb. 3. Auch hier ergibt sich, daB 
bei Zuscharfung der Vorderkanten die nichtlinearen Anteile etwa doppelt so groB sind wie bei 
abgerundeter Vorderkante, ein Befund, der auch von anderen Autoren!:? bestatigt worden ist. 

Da aus den vorhandenen MeBergebnissen hervorgeht, daB durch Abrundung der Vorderkanten 
die nichtlinearen Anteile der Kraftbeiwerte gerade auf die Halfte absinken, sollen die folgenden 
Untersuchungen der Kinfachheit halber auf solche Tragfliigel beschrankt werden, bei denen die 


Vorderkanten zugescharft sind, da die Umrechnung auf abgerundete Vorderkanten keine Schwierig- 
keiten bereitet. 


1G. E. Bartlett u. R. J. Vidal, Journ. Aero. Sci. 22 (1955), S. 517. 

* A. H. Flax u. H. R. Lawrence, The aerodynamics of low-aspect-ratio wings and wing-body combinations; 
Anglo-American Aeronautical Conference, Brighton 1951, S. 363. : 

3 H, Winter, Forschung Ing.-Wes. 6 (1935), S. 40 u. S. 67. 
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c) Induzierter Widerstand. Der Unter- 
schied zwischen runder und zugescharfter Vor- 
derkante ist aerodynamisch dadurch gekenn- 
zeichnet, daB bei runder Vorderkante durch die 
Umstrémung Saugkrafte entstehen, wahrend bei 
zugescharfter Vorderkante infolge Ablésung der 
Strémung diese Saugkrafte fortfallen. Bei dem 
im folgenden immer zugrunde gelegten Fall zu- 
gescharfter Vorderkante wird daher bei reibungs- 
loser Strémung die resultierende Luftkraft stets 
senkrecht auf der Fliigelebene stehen. Aus dem 
Beiwert dieser Normalkraft cy erhalt man dann 
den Auftriebsbeiwert zu cy = cy cosa und den 
Widerstandsbeiwert zu cy; = cy sina. Fiir kleine 
Anstellwinkel ergibt sich daraus naherungsweise 


(1) 
(2) 


C4 
Cri 


Da man also aus dem Auftrieb sofort den Wider- 
standsbeiwert erhalt, sind im folgenden alle Be- 
trachtungen auf Auftrieb und Kippmoment be- 
schrankt. 


d) Einflu8 der Kompressibilitat. Die 
nichtlinearen Effekte treten besonders stark auf 
bei Tragfliigeln mit kleinem Seitenverhaltnis. Nun 
ist aber bekannt, daB die Strémung um schlanke 
Kérper, also auch um Tragfliigel mit kleinem 
Seitenverhaltnis, unabhangig von der Machzahl ist. 


ie 70° 15° 20° 02 30° 35° 40° 
Abb. 3. Einflu8 der Vorderkantenform auf Auftriebs- 
A=1,0,9 = 45°. — @ Vorderkante rund; ——O 


und Momentenkurven. — a) Deltafliigel, 4 = 0,78, 
Vorderkante scharf; 


ie 
08 


7 
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b Tye 4° Sif Wa 16 20° 24° 
Abb, 2, EinfluB der Form der Vorder- und Seitenkanten von Recht- 
eckfliigeln auf die Auftriebskurve. — a) Seitenverhaltnis A = 0,5; 
b) Seitenverhaltnis A = 1,0. — Messungen nach H. Winter (s. Fub- 
note 3 von S, 432) — e Vorderkante und Seitenkante 
rund; O—— Vorderkante und Seitenkante scharf; A 
Vorderkante rund, Seitenkante scharf; —  — Limeare Theorie 


nach E. Truckenbrodt. 


16° 80° = 24° 28° 38° 


= 0,125; b) Pfeilfligel, A=1,0, 
nach E. Truckenbrodt. 
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Dieses Ergebnis aus der von R. T.. Jones! begriindeten linearen Theorie schlanker Fligel (slender-body- 
theory) wird auch von den Experimenten gut bestatigt. Bei Uberschallanstrémung gilt der Fligel dann 
als schlank, wenn er sich etwa in der Achse des Mach-Kegels befindet und seine Vorderkanten 
weit genug vom Mach-Kegel entfernt sind. Da mit zunehmender Machzahl der Kegel spitzer 
wird, mu8 auch das Seitenverhaltnis des Fliigels bei gleichen Schlankheitsforderungen kleiner 
werden. Daf trotz der Uberschallanstrémung die Strémung um den Fliigel gleich der inkompres- 
siblen Strémung ist, beruht darauf, daB infolge der Schlankheit des Fliigels die wirksame Ge- 
schwindigkeitskomponente, z. B. senkrecht zur Vorderkante eines Deltafliigels, nur sehr klein ist 
und daher zu gleichen StrémungsverhaJtnissen wie im inkompressiblen Fall fiihrt. 

Da also die kompressible Strémung bei sehr schlanken Tragfliigeln im wesentlichen mit der 
inkompressiblen Strémung tibereinstimmt und bei Tragfliigeln mit groBem Seitenverhaltnis wenig- 
stens im Unterschallbereich auf den inkompressiblen Fall zuriickgefiihrt werden kann, sollen die 
folgenden Untersuchungen auf inkompressible Strémung beschrankt werden. 


3. Bisherige theoretische Arbeiten. Die meisten theoretischen Verfahren zur Berechnung von 
Auftriebsverteilungen an Tragfliigeln sind lineare Verfahren, d.h. ihre Anwendung ist auf kleine 
Anstellwinkel beschrankt. Man unterscheidet bei ihnen zwischen den einfacheren Traglinien- 
verfahren und den genaueren, dafiir komplizierteren Tragflachenverfahren. Da die Trag- 
linienverfahren ausschlieBlich auf Fliigel mit groBen Seitenverhdltnissen zugeschnitten sind, 
werden im folgenden nur Tragflachenverfahren benutzt. Die wichtigsten sind die beiden fiir 
beliebige Fliigelgrundrisse geltenden Verfahren von H. Multhopp? und E. Truckenbrodt?. | 
Das Verfahren von E. Truckenbrodt hat den Vorteil, daB es im Grenzfall sehr kleinen Seitenver- 


Betz 
(1935) Legendre (1952) 
Cross -Flow-Theorie ae (1900) 
angler (7956, 


Abb. 4. Wirbelmodelle der verschiedenen bestehenden nichtlinearen Tragfligeltheorien fiir Fliigel mit kleinen Seitenverhiltnissen. 


haltnisses die gleichen Ergebnisse liefert wie die Theorie schlanker Kérper von R. T. Jones}. 
Dieses Verfahren von E. Truckenbrodt wurde daher als Grundlage fiir das weiter unten entwickelte 
Berechnungsverfahren gewahlt. 

Versuche, auch fiir die nichtlinearen Anteile Berechnungsverfahren zu entwickeln, sind von 
verschiedenen Autoren insbesondere in den letzten Jahren unternommen worden. Es lassen sich 
diese Verfahren in drei groBe Gruppen einteilen (vgl. Abb. 4). 

. i nie NACA-Report 835 (1946), S. 59. 

- Multhopp, Method for calculating the lift distribution of wings (Subsonic liftine- 
ane RM 2884 (1955); vgl. auch H. Multhopp; Luftfahrtforschung Gee Stepaaemencnc 
E. Truckenbrodt, Jahrbuch 1953 der Wiss. Ges. f. Luftfahrt, S. 40. 


— 
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a) Gross Flow Theory. Die ersten rechnerischen Ansatze stammen von A. Betz! und 
sind heute unter dem Namen cross-flow-theory bekannt?. Bei dieser rein empirischen Methode 
wird davon ausgegangen, da® die seitliche Querumstrémung des Fliigels an Bedeutung gewinnt 
falls die Spannweite klein ist gegentiber der Fliigeltiefe. Der nichtlineare Anteil fiir den N nae 
kraftbeiwert ergibt sich dann im wesentlichen aus dem Widerstandsbeiwert dieser nahezu zwei- 
dimensionalen Querstrémung (Abb. 4a). Meistens wird dabei den Rechnungen der Widerstands- 
beiwert der sehr langen quergestellten Platte Cy = 2,0 zugrunde gelegt. 


b) Rechteckfliigel von kleinem Seitenverhaltnis. Die zweite Theorie stammt von 
W. Bollay®. Diese ist beschrankt auf Rechteckfliigel mit sehr kleinem Seitenverhaltnis. Dabei 
wird von einem Wirbelmodell ausgegangen, bei welchem die freien Wirbel nicht mehr in einer 
Ebene liegen, sondern an den Seitenkanten jeweils unter dem Winkel a/2 zur Fligelebene nach 
hinten abgehen (Abb. 4b). Da der Neigungswinkel gerade «/2 sein mu, laBt sich so einsehen: 
Da die Strémungsrichtung am Fliigel parallel zur Fliigelebene liegt (kinematische Strémungs- 
bedingung), mu die Strémungsrichtung beim Durchgang durch die seitlichen Wirbelschichten 
um den Winkel « gedreht werden. Bekanntlich liegt der Wirbelvektor in der mittleren Strémungs- 
richtung, folglich miissen die Wirbelvektoren unter dem Winkel «/2 zur Fligelebene geneigt sein. 
W. Bollay macht nun auBerdem die Annahme, da8 diese Richtung auch weiter stromabwarts ein- 
gehalten wird, so daB das seitliche Wirbelband aus geraden Wirbelfaden besteht. Da die Wirbel- 
starke lings Spannweite konstant angesetzt ist, bleibt jedoch die Anwendung auf sehr kleine 
Seitenverhaltnisse beschrankt. 


c) Schlanke Deltafliigel. In neuerer Zeit wurde von mehreren Autoren?: > 6 7,8 910 die 
Strémung um sehr schlanke Deltafliigel untersucht, bei denen sich infolge der Vorderkanten- 
ablésung zwei Wirbelschichten an den Vorderkanten entwickeln, die sich iiber der F liigel- 
flache in zwei mehr oder weniger konzentrierte Wirbel aufrollen (Abb. 4c). Von den genannten 
Autoren wird diese Strémung durch ein geeignetes Wirbelmodell ersetzt, was dann auch die nicht- 
linearen Glieder fiir die Kraftbeiwerte liefert. Der Grund fiir die Nichtlinearitat liegt darin, daB 
wie im Wirbelmodell von W. Bollay die Wirbel nicht mehr in einer einzigen Ebene liegen. 

Alle genannten Verfahren sind auf kleine Seitenverhaltnisse beschrankt (slender body theory). 
Sie gestatten daher keine Aussagen iiber den Einflu8 des Seitenverhaltnisses A auf die GréBe der 
nichtlinearen Glieder. Dariber hinaus sind sie nur auf bestimmte Fligelformen (Rechteckfligel 
oder Deltafliigel) oder spezielle Strémungen (z. B. konische Strémungen) anwendbar. 


4, Wirbelmodel und Grundgleichungen des Berechnungsverfahrens. Zur Aufgabe ist gestellt, 
fiir einen Tragfliigel von beliebiger Grundrifform und mit zugescharften Vorderkanten die Auf- 
triebsverteilung und die aerodynamischen Beiwerte fiir inkompressible Strémung zu berechnen 
mit besonderer Beriicksichtigung der nichtlinearen Abhangigkeit der Beiwerte vom Anstellwinkel. 
Der Einfachheit halber habe der Fliigel symmetrisches Profil. 

Obwohl bei der Umstrémung von Fligeln mit kleinem Seitenverhaltnis der ReibungseinfluB 
in der Strémung mafgeblich beteiligt und insbesondere fiir die Ablésung und Wirbelbildung an 
den Kanten verantwortlich ist, soll im folgenden bei der theoretischen Behandlung solcher Trag- 
fliigelstrémung reibungslose Flissigkeit vorausgesetzt werden. In der Tragfliigeltheorie und auch 
bei anderen Strémungsproblemen (z. B. Sekundarstrémung in gekritmmten Rohren oder Kanilen, 
Sekundarstrémung hinter Schaufelgittern) ist es iiblich, das unter dem EinfluB der Reibung ent- 
standene Strémungsfeld experimentell genau zu erforschen und danach méglichst genau durch 


1 4, Betz, Applied Airfoil Theory; in W. F. Durand: Aerodynamic Theory, Bd. IV, Berlin 1935. 

2 Siehe FuBnoten 1 und 2 von Seite 432. 

3 W. Bollay, Z. angew. Math. Mech. 19 (1939), S. 21; vgl. auch Journ. Aero. Sci. 4 (1937), S. 294. 

4C. E. Brown u. W. H. Michael, Journ. Aero. Sci. 21 (1954), S. 690; vgl. NACA TN 3430 (1955). 

5 P, T. Fink, Z. Flugwiss. 4 (1956), S. 247. 

6 D. Kiichemann, RAE-Report Aero 2540 (1955). eo 

7 R. Legendre, Ecoulement au voisinage de la pointe avant d’une aile a forte fleche aux incidences moyennes. 
La Recherche Aéronautique; Bulletin Bimestrial de 1’Office National d’Etudes et de Récherches Aéronauti- 
ques (1952). Af 

8 K. W. Maneler u. J. H. B. Smith: Calculation of the flow past slender delta wings with leading edge sepa- 
ration. Actes de 1X. Congres Intern. de Mechanique Appliquee, Bd. IT, S. 137 (1956); vgl. auch RAE TN 2593 
sowie K. W. Mangler u. J. H. B. Smith; Proceed. Roy. Soc. A 251 (1959), S. 200. 

9 M. Roy, Z. angew. Math. Phys. 9b (1958), S. 554 (Ackeret-Festschrift); vgl. auch M. Roy; Recherche 

é i : 57), 5. 3. P : : 
ceie D. taeda Ry Simpson u. W. I. Rainbird, The flow over delta wings at low speeds with leading 
edge separation. The College of Aeronautics, Cranfield, Report No. 114 (1958). 
30* 


Ingenieur-Archiv 


436 K. Gersten: Nichtlineare Tragflachentheorie 


eine reibungslose Modellstrémung wiederzugeben. Die Kuttasche Abflu8bedingung bei Unter- 
schallhinterkanten ist das bekannteste Beispiel fiir einen Reibungseffekt, der trotzdem in der 
reibungslosen Betrachtungsweise beriicksichtigt wird. 

In der linearen Tragfliigeltheorie ist es iiblich, einen gegebenen Tragfliigel zu ersetzen durch 
eine groSe Anzahl von Elementarfliigeln der infinitesimalen Spannweite dy (Abb. 5b). Dieser 
Elementarfliigel hat eine tiber der Tiefe veranderliche Zirkulationsverteilung und nach hinten ab- 
gehende freie Wirbel, die jedoch in der Fliigelebene liegen (Abb. 5a). Dieses bekannte klassische 
Wirbelmodell wird jetzt dahingehend abgedndert, daf fiir die Elementarfliigel das Bollaysche 
Wirbelmodell zugrunde gelegt wird. Die freien Wirbel liegen jetzt also nicht mehr in der Fliigel- 
ebene, sondern unter dem Winkel «/2 zu dieser (Abb. 5c). Bei dem Gesamtfliigel geht dann von 


Elementarfligel Elementartlige! 
nach Modell von nach Modell von 
E. Truckenbrodt W. Bollay 


sk 


_ Abb. 5. Wirbelmodelle fiir lineare und nichtlineare Tragflachentheorie, — 
a) Lineare Theorie; b) Aufteilung des Fliigels in Elementarfliigel; c) Nichtlineare Theorie. 


jedem Punkt der Fliigelflache, in der eine Anderung der Wirbelstarke vorliegt, ein freier Wirbel 
unter dem Winkel «/2 nach hinten ab. Da dieser Winkel mit abnehmendem Anstellwinkel kleiner 
wird, geht die ganze Wirbelanordnung im Grenzfall « —> 0 in das tibliche Wirbelmodell der linearen 
Theorie iiber. Die Hauptaufgabe besteht jetzt darin, fiir das so definierte Wirbelmodell die indu- 
zierten Geschwindigkeiten an der Fliigelflache zu ermitteln. Damit lassen sich dann die kine- 
matischen Strémungsbedingungen aufstellen, welche die Bestimmungsgleichungen fiir die Wirbel- 
verteilung liefern. Die kinematische Strémungsbedingung verlangt, daB die Richtung der resul- 
tierenden Geschwindigkeit, die sich aus der Grundstrémung und der induzierten Geschwindigkeit 
an der betreffenden Stelle zusammensetzt, tangential zur Fliigelflache verlauft. Da die Horizontal: 
komponente der induzierten Geschwindigkeit klein ist gegen die Anstrémungsgeschwindigkeit, ist 
der Winkel der resultierenden Geschwindigkeit gleich dem induzierten Abwindwinkel a«,, = — e V. 
wenn w die Vertikalkomponente der induzierten Geschwindigkeit darstellt. Die atice 
Strémungsbedingung besagt danach, dafi der induzierte Abwindwinkel gleich dem Ortlichen 
Anstellwinkel des Tragfliigels sein muB®, d. h. die kinematische Strémungsbedingung fordert: 


Oy = Karel - (3) 


Belegt man die Tragfliigelflache mit einer kontinuierlichen Wirbelbelegu 
| ng k(x, {eal Bax 
ergibt sich nach der linearen Theorie fiir den induzierten desea er or anrtiaeai ates 


Oy = L,(k) ? (4) 
1K. Gersten, Jahrbuch 1957 der Wiss. Ges. f. Luftfahrt, S. 172. 
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und damit lautet die kinematische Strémungsbedingung 
sr = L,(k)  (lineare Theorie) . (5) 
Dabei ist L,(k) ein linearer Operator mit folgenden Eigenschaften: 
Ey(ck)=cLy(k), — Ly(ky + hy) = Ly (hy) + Ly (hp) « (6) 


Im allgemeinen Fall der tragenden Flache ist dieser Operator cin kompliziertes Doppelintegral, 
das sich tiber die Fliigelflache erstreckt. Geht man nun zu dem beschriebenen nichtlinearen Wirbel- 
modell iiber, so erweitert sich die Gleichung (4) um ein additives Zusatzglied, welches nun auch 
vom Anstellwinkel abhangt. Entwickelt man dieses Zusatzglied nach Potenzen von qa, so ergibt 
sich fiir den induzierten Anstellwinkel in der Fligelebene anstelle von (4) jetzt 

a 


Oy = Ly(k) + 2 0 La(k) fee. (7) 


|a| 
Dabei hat L, ebenfalls die Eigenschaften (6). Die kinematische Strémungsbedingung (3) fihrt 
damit zu der folgenden Bestimmungsgleichung fiir die Zirkulationsverteilung k: 


Ore = L,(k) + 6 L,(k) +--+. (nichtlineare Theorie) . (8) 


|a| 


Die Lésung dieser Gleichung erfolgt durch den Ansatz 


k=ka+hk, al a2, (9) 


wobei k, den linearen und k, den nichtlinearen Anteil der Zirkulationsverteilung darstellt. Man 
erhalt dann unter Benutzung der Beziehungen (6) die beiden Gleichungen 


Lika (10) 
L,(k,) = — 2 1). (11) 


Die Gleichung (10) fiir k, ist identisch mit dem Problem (5) der linearen Theorie. Den nichtlinearen 
Anteil k, erhalt man aus der ganz analog gebauten Gleichung (11), die sich von (10) nur durch 
die veranderte rechte Seite unterscheidet. Hat man aus (10) die Funktion k, ermittelt, so laBt sich 
in sehr einfacher Weise die rechte Seite von (11) bilden. Es stellt sich namlich heraus, daB der 
Operator L,(k) viel einfachere Gestalt hat als L,(k). 

Sind die Funktionen k, und k, ermittelt, so erhalt man daraus fur die Kraftbeiwerte c, und cy 
die folgenden Beziehungen: 


2 

ca= en =p | | Heoy) ded, | (12) 
F 
2 
eu=— arp | | # Hey) de dy. (13) 
; F 
wobei F die Fliigelflache, V die Anstrémgeschwindigkeit, ferner 
Lee 
La 4 [rors 


die Bezugstiefe mit der drtlichen Fliigeltiefe I(y) und s die Halbspannweite bedeuten. 

Die Berechnung der induzierten Geschwindigkeiten, d.h. die Bestimmung der Ausdriicke 
L,(k) und L,(k), hangt nun davon ab, von welchem Berechnungsverfahren der linearen Theorie 
man ausgeht. Fiir das Tragflachenverfahren von E. Truckenbrodt werden im folgenden die Be- 


rechnungsformeln im einzelnen aufgestellt. 


5. Das Berechnungsverfahren. a) Lineare Tragflachentheorie von E. Truckenbrodt. Das 
Verfahren zur Berechnung der Auftriebsverteilung von E. Truckenbrodt kann auf Tragfligel mit 
beliebigem Grundri8 angewendet werden. Die Filiigelflache wird dabei mit der kontinuierlichen 
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Wirbeldichte k(x’, y’) belegt. Das bedeutet, daB von jedem Punkt der Fliigelflache ein elemen- 
iarer Hufeisenwirbel der Starke k(x’, y’) dx’ abgeht (vgl. Abb. 6a). Alle Hufeisenwirbel, die bei der 
Linie x’ = const. beginnen, bilden eine ebene Wirbelschicht mit der Zirkulationsverteilung k(x’, y') dx’, 
die nach (76) des Abschn. 7 im Punkt P(x, y, 0) die folgende Geschwindigkeit in z-Richtung 
induziert: 


| dw(x, y¥) = 
V 


+s 


er 9(t +2 )] 
y—y’ 
x dx' dy’. (14) 


d 
dy’ 


Dabei ist 
R=/@—«P+(y—yy- (15) 
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Abb. 6. Bezeichnungen am Pfeilfliigel 
(Tragflaichentheorie von E. Truckenbrodt), — 
a) GrundriB; b) réumliche Darstellung (von dem abgehenden 
Wirbelband ist nur eine infinitesimale Wirbelschicht 
eingezeichnet). 


Den vom gesamten Tragfliigel induzierten Abwindwinkel «,, = — w/V erhalt man daraus, indem 
man iiber die Wirbelschichten, die jeweils bei x’ beginnen, in Tiefenrichtung, d. h. in x’-Richtung, 
integriert. Nach Vertauschen der Integrationsreihenfolge ergibt sich dann 


+s “h A me (1 re ae 
Bas Ae 1 dy’ sl R | rps 
a,(k) = —% = +aar¢ | ss dx! dy’. (16) 


Ss A) 


Dabei sind x,(y') die Abszisse der Vorderkante und x,(y') = x,(y’) + [(y') die der Hinterkante. 
Fir die Wirbelverteilung k(x, y) wird jetzt der folgende Ansatz gewahlt: 


k(y, 3) = sf a po) ctg + + 4 u(?) (cre + —2sin #)) (17) 
mit 
1 rack te (orl 
yl — eos g) = =F) (18) 
und 
y’ 


cos } = y= (19) 


& . 
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Die Verteilung der Wirbeldichte in Tiefenrichtung fiir einen festen Fliigelschnitt y' = konst. 
entspricht den beiden ersten Birnbaumschen N ormalverteilungen. Die Starke der jeweiligen 
Tiefenverteilung andert sich langs der Spannweite. Sie wird durch die Zirkulationsverteilung 
y(0) baw. die Momentenverteilung (9) bestimmt. Statt der kontinuierlichen Verteilung k(x, y) 
werden nur noch die beiden Funktionen y(#) und L(A) gesucht. 


Setzt man (17) in (16) ein, so erhalt man bei Verwendung der dimensionslosen Koordinaten 
pert WR ra * ee 
sla bf rae =—, Ls Sea (20) 


den induzierten Anstellwinkel zu 


M(F59) =— a> = Lyn) = 5 roy (21) 


4 
—1 


aes; : 
a 1 oe 


Dabei wurden die beiden EinfluBfunktionen 


Xh 
: 2 y x—x'\,, 
(x, Yss i ote $ (1 +45 jae (22) 
“h 
: 8 See 
V(X, Li ten | (cts $—2 sin) (1 + “Re ) a (23) 
mit ; 
aes x —*(¥.') ye | ome i 2A 
ii Wes Ky’) oe 
eingefiihrt. 


Die kinematische Strémungsbedingung ergibt sich aus (21) zu 


‘cue we (iy +i¥) 
a(n) = Ly, m) = taf a (25) 
Da zwei Funktionen y(7) und u(y) gesucht werden, mu die Bedingung (25) auch auf zwei ver- 
schiedenen, besonders geeigneten Anstellinien auf dem Fliigel erfiillt werden. Bei E. Truckenbrodt 
sind dafiir die Hinterkante (durch einen Stern gekennzeichnet) und die 1/4-Linie (durch zwei Sterne 
gekennzeichnet) gewahlt worden. Man erhalt daher das Gleichungssystem 


d 
+1 ° a rf} 
l f G*y +j* pu) 


60(H) = L*(y, y= 2a eae) dy’  (Hinterkante) , (26) 
ei 
ape F 
1 dy! (i** y + j** u) Bi 
a(n) = LW = ag a (iL). (27) 


al 
Setzt man nach dem Vorbild von H. Multhopp fiir die unter den Integralen stehenden Funktionen 
i* y, i** y, j* u,j** u Sinusreihen an, so lassen sich die Integrale durch geeignete Summenformeln 
ausdriicken. Werden die Ansatze y = g, a, 4 =m, gemacht, dann erhalt man nach Trucken- 
brodt aus (26) und (27) 
i. 3/2 be : s 
a, y= = 2 ae b, (;;) |e Wz a Gyn Un Bin — 12 b, (i) 


3/2 m : 
Ue pa, Gyn Jon Mn » (28) 


m 2 m A 
ee 1.2180 a () Bee i inthe + 3.3000 +60, (7) Jmu— ¥ ay, j2t my (29) 
‘ a v Fl v n=1 
(py Sat, 2) eo cm) 
Dabei sind die a,,, 4,,, 6, und c, universelle, yon Truckenbrodt tabulierte Konstante. Die g;, 
und m,,, bzw. g,, und m,, sind die Funktionswerte an den Stellen 


haw. eeeeey < (30) 


Nn = COs tas fig" : 
Hs m+1 m+1 
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Die Zahlen I, sind die Fliigeltiefen an diesen Stellen. Weiterhin bezeichnet m die Anzahl der ge- 
wahlten Elementarfliigel langs der Spannweite. Der Strich am Summenzeichen deutet darauf hin, 
daB bei der Summation das y-te Glied ausgelassen werden soll. 

Die Gleichungen (28) und (29) stellen ein gekoppeltes lineares Gleichungssystem fiir die 2 m 


Unbekannten g,, und m,, (v= 1, 2,..., m) dar. Sind diese Werte gefunden, so erhalt man daraus 
die aerodynamischen Beiwerte, wie folgt: 

d td s 

(Zz) =A f s(n) dy, (31) 

ao=0 —l1 
SLi Spl 
d I : 
den) A [ gylq) dy + A [ mn) FP a (32) 
% ao=0 “ a 
a —1 
wobei ; +s 
L== | I? dy (33) 


—s 
ist. Auch fiir diese Integrale lassen sich einfache Summenformeln angeben, wenn man die folgenden 
Quadraturformeln verwendet: 


[ fen an = D'S) sin 2 (34) 


+1 

1 ¢. dyed me 

Oi ld; = pe Onda) (35) 
Bed 


b) Berechnung der Auftriebsverteilung fiir den nichtlinearen Anstellwinkel- 
bereich. Wie schon in Abschn. 4 ausgefiihrt wurde, besteht der Unterschied zwischen der nicht- 
linearen und der linearen Theorie nur in einer Abanderung des Wirbelmodelles. Gegeniiber dem 
Modell der linearen Theorie, bei dem alle abgehenden freien Wirbel in einer Ebene, der Fliigelebene, 
liegen, bilden beim Modell der nichtlinearen Theorie die abgehenden freien Wirbel den Winkel 
«/2 mit der Fliigelebene (vgl. Abb. 6b). Es seien nun alle elementaren Hufeisenwirbel betrachtet, 
die ihren Ursprung bei der Koordinate x’ haben. Wie in Abb. 6b angedeutet, bilden sie zusammen 
eine ebene Wirbelschicht mit der Wirbelverteilung k(x’, y’) dx’ langs der Spannweitenkoordinate y’. 
Der Punkt P(x, y) der Fliigelflache hat von dieser Schicht den Abstand 2 = — (x — x’) «0/2. In 
Abschn. 7, Gleichung (80), wird gezeigt, da8 im Punkt P(x, y) von der bei x’ beginnenden Wirbel- 
schicht die folgende Abwindgeschwindigkeit induziert wird: 


d x—x! 
svob fcasilag SE 
itr RE A R toon 
2 ee nm ¢ ay dy' dx 
a Leta x— x! ie 
= ocale saat) O—*) akon) oy pe 
Fiihrt man jetzt die Integration iiber den gesamten Fliigel von der Fliigelnase bis zur Hinterkante 
durch, dann erhalt man fiir den gesamten induzierten Abwindwinkel in P (x, y, 7 = — («/2) x) 
nach Vertauschen von Integration und Differentiation 
«, d — x!) 
Ts (kh |k(x’,y') ( ues * 
1 d 7 | oY 
w(x, y) = iced a ¢ | IA aad R dx! dy’ 
lL @ a ? d 
x—x 


es man wieder denselben Ansatz (17) fiir k(x’, y') und fiihrt man die dimensionslosen Ko- 
ordinaten (20) ein, so erhalt man fiir die kinematische Strémungsbedingung: 


a(&,4) = L(y, w) + Ta ° he HD (38) 
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Dabei sind L(y, u) und J(y, 4) lineare Operatoren mit den Kigenschaften (6) 


. L(y, uw) ist identisch 
mit dem Operator (25) der linearen Theorie. Fiir J(y, bs) gilt ae 


ert 
Is) = —“ ty) Sula) + uly) Jer) (39) 
mit 
A 
dG) = = i (1 + cos ¢’) (cos y’ — cos ¢) dy’ 
if (40) 
J,(n) = oa | (1 + cos y' — 2 sin? :) (cos y' — cos ¢) dq’ , 
6 
wobei 
— AE za(l cos g)s 
Sin! | (41) 
ae “ (1 — cos 9’) 
gilt. Zweimalige Differentiation von J(y, 1) liefert 
2 2 dx, dl dx, 
= filo) + Ao B) ge + Sal) + Hw S > 
+ Faby 0) (Ge) + $a (Fe) (42) 


Dabei geben die Glieder mit dl/dy den EinfluB der Zuspitzung des Fliigels und die Glieder mit 
dx,/dy den EKinfluB8 der Pfeilung des Fliigels wieder. Die Operatoren f,(y, u) bis f,(y, 4) ergeben 
sich dabei durch elementare Integrationen und Differentiationen. Bezeichnen wieder ein Stern 


die Werte fiir die Hinterkante und zwei Sterne diejenigen fiir die 1/4-Linie, so gelten die folgenden 
Formeln: 


fiir die Hinterkante 


3 Idy 1 dm 
fie.» = 4 b dy? b dn?’ 
1d d 
“ d 43 
fom =—2, < 
fi(y, 4) =9, 
Fiy.u) =, | 
fély, uw) =9, 


fiir die 1/4-Linie 


fi*(y, LB) = 373 1 dy a l i 


ert ete 
ron =—(E—-H) (22) 
fry =—(44 Be 8B &, ao 
fi*(y. mb) = es 
few =F, 
from =r, 
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Da die Funktionen y und y im allgemeinen durch ihre F unktionswerte y,, und z,, an den m diskre- 
ten Stellen 7,, = cos [nz/(m + 1)] (n=1,...,m) gegeben sind, lassen sich die erforderlichen, 
ersten und zweiten Ableitungen mittels folgender Summenformeln bestimmen: 


m 
a ee eS (45) 
dy n=1 

m 
ot 46 
dy? a Dy Aunt ( ) 


Die universellen Summenkoeffizienten A,, und A,,, sind fir m = 7 und m = 15 in den Tabellen 1 


bis 4 tabuliert worden. 


Zur Bestimmung der gesuchten Zirkulationsverteilung y und Momentenverteilung u wird nun 
die kinematische Strémungsbedingung (38) nach dem Vorbild der linearen Theorie auf der 1/4-Linie 
(zwei Sterne) und der Hinterkante (ein Stern) erfiillt. Es ergeben sich dann aus (38) die folgenden 
beiden Bestimmungsgleichungen fiir y(y) und p(y): 


a*(n) = L*(y, uw) + 


a**(n) = L**(y,p) 4 


— 0 
|) 


a | @PI*(y, 4) 


dy? 


|&| 


a | PJ**(y, 4) 
a 2 


dy 


(Hinterkante) , 


(1/4-Linie) . 


Tabelle 1. Summenkoeffizienten A,n nach Gl. (45) fiir m = 7. 


1 0,9239 0,4142 —0,7654 2,4966 | 3,1544 
2 0,7071 S19 2,2992 0,7072 —8,5238 
3 0,3827 | 2,4142 0,2242 —4,0274 4,4608 
4 0 0 —2,8286 2,0000 —2,8284 
5 —0,3827 —2,4142 1,3066 —1,1990 1,8478 
6 —0,7071 1,0 —0,7024 0,7072 —1,1330 
Zi —0,9239 —0,4142 0,3170 | —0,3318 | 0,5412 


Tabelle 2. Summenkoeffizienten A,n nach Gl. (45) fir m = 15. 


(47) 


(48) 


0,1951 0,3827 0,5556 0,7071 0,8315 0,9239 0,9808 
1 0,9808 0,1990 | —0,2530 0,3530 | —0,5518 1,0082 | —2,3516 8,9574 | 12,8848 
2 0,9239 | —0,4142 0,5354 | —0,7654 1,2496 | —2,4966 71,4540 3,1542 | —34,4610 
3 0,8315 0,6682 | —0,8902 1,3400 | —2,4218 6,3180 1,3470 | —15,7096| 19,0720 
4 0,7071 | —i,0 1,4080 | —2,3592 5,6120 0,7072 | —10,2346)  8,5238 | —13,2442 
5 0,5556 1,4966 | —2,3518 5,2058 0,4018 | —7,7596 5,4246 | —5,8990 | 10,0232 
6 0,3827 | —2,4142 5,0212 0,2242 | —6,4272 4.0274 | —3,7054 4,4608 | — 7,9178 
W 0,1951 5,0274 0,1014 | —5,6588 3,2722 | —2,7090 2,7740 | —3,5166 6,3986 
8 0 0 —5,2262 2,8286 | —2,1648 2,0000 | —2,1648 2,8284 | — 5,2262 
9 —0,1951 | —5,0274 2,9630 | —1,8374 1,5714 | —1,5374 1,7198 | —2,2904 4,2754 
10 —0,3827 254142 | —1,6304 1,3066 | —1,1842 1,1990 | —1,3696 1,8478 | — 3,4732 
Il —0,5556 | —1,4966 1,1294 | —0,9592 | 0,9000 | —0,9312 1,0790 | —1,4686 2,7740 
12 —0,7071 1,0 —0,7992 0,7024 | —0,6736 0,7072 | —0,8272 1,1330 | — 2,1474 
13 —0,8315 | —0,6682 0,5518 | —0,4952 0,4818 | —0,5106 0,6014 | —0,8270 1,5714 
14 —0,9239 0,4142 | —0,3488 0,3170 | —0,3110 0,3318 | —0,3924 0,5412 | — 1,0298 
15 —0,9808 | —0,1990 0,1692 | —0,1548 0,1528 | —0,1634 0,1738 | —0,2676 0,5098 
Tabelle 3. Summenkoeffizienten A, nach Gl. (46) fiir m = 7. 
n ™n 7] | 0 0,3827 0,7071 0,9239 
1 0,9239 0,8968 — 3,1718 15,9588 —103,6032 
2 0,7071 — 2,8284 15,6018 —25,0004 24,8700 
3 0,3827 12,6174 | —24,4024 8.5256 11,6590 
4 0 —21,0000 13,5142 —2,8284 = TNE TOS 
5 —0,3827 12,6174 — 2,8284 | 11,1114 8,8288 
6 —0,7071 — 2,8284 0,9740 —15,0002 — 5,7578 
7 —0,9239 0,8968 0,3432 10,5444 2,8286 
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Tabelle 4. Summenkoeffizienten A,,, nach Gl. (46) fiir m = 15. 


2 Ee eee 
ot BE ” 0 0,1951 0,3827 0,5556 0.7071 ~—si| 0,8315 0,9239 0.9808 
ce ee i 
1|  0,9808 0,4056| — 0,6958] —1,8054. 4,5940 8,7930 | — 21,8362| 371,3938 —1569,2914 
2] 0,9239 — 0,8968) — 1,5780| — 4,7588 — 3,0038) — 26,5648; 181,4036 —540,6458| 323.2198 
31 0,8315 1,6072/— 2,9778, _9,0586 — 11,8702] 10,5352 —268,1322| 240.9128 236.0254 
4 | 0,7071 |— 2,8284) —_5,7854) —17,1472, _ 78,5826 | —167,9976| —137,0206 | — 24.8700 — 244.4994 
5 10,5556 5,3876 —13,5242| 61,2484 —129,9428| 91,4382 — 24,7094 5,1870| 211.1428 
6 | 0,3827 )—12,6174| 54,5522' 95,6066! 79,9826 | — 19,1318 6,5300| 11,6640 |— 177,5580 
7] 0,1951| 51,5388] —88,3206) 54,4708 — 23,1596 6,7508|— 1,2444| — 12,5380 148.6024 
8 0 —85,0000| 52,5168) —13,5134 6,0540|— 2,8284|— 0,6256| 11,7216 '— 124,0242 
9 | —0,1951 | 51,5388] 12,6172) 8,8596. 55,0688 1,2338 1,2830 | — 10,3552) — 102,9086 
10 | —0,3827 | —12,6174]  5,3130 — 6,6018 — 9,2216/— 0,5046|— 1,4346| 8,8288 — 84,4126 
11 | —0,5556 5,3876| — 2.7798 29180 6,7362 0,1578 1,3510|— 7,2782| 67,8734 
12 | —0,7071 | — 2,8284, 11,6094 —_0,5722] 0,5258 0 = 1.1528 5,7578'— 52,8010 
13 | —0,8315 | - 1,6072|— 0,9632 — 1,8930|— 7,9404|— 0,0586 0,8966|— 4,2750 38,7770 
14 | —0,9239 | — 0,8968}  0,5526! 11,2434 + ~—«:10,9652 0,0624|— 0,6100 2,8286 |— 25,4718 
15 | —0,9808 0,4056| — 0,2534! — 0,3092;— 7,7088|— 0,0378 0,3088|— 1,4084 12,6196 


Dabei sind die Operatoren L*(y, ~4) und L**(y, uw) identisch mit (26) und (27) der linearen Theorie 
Die Werte d?J*/dy? und d?J**/dy? ergeben sich nach (42) aus den in (43) und (44) angegebenen 
Funktionen. 

Zur Lésung dieses Gleichungssystems wird der folgende Ansatz gemacht: 


v(n) = B1(n) & + Bln) Fa 
(49) 
p(n) = myn) + ma(n) To. 
Geht man damit in die Gleichungen (47) und (48) ein und benutzt die Beziehung 
Lc; 81 + ¢ 82» Cy My + Cy Mg) = c, L(gy, m,) + Cy Lge, ms), (49 a) 


so erhdlt man unter Vernachlassigung der Glieder, die von der dritten Potenz des Anstellwinkels 
abhangen, 


meio L7*(g,,m,) — = (50) 
und 
* d2J* A . a**(n) as s*(g A m. ) 
L*(go, me) ae toot Si a CL ey my) ; i= *(2., My) — ee cae es a (51) 


Da die Operatoren L* und L** nach (26) und (27) Integrale iiber die gesuchten Funktionen 
sind, stellen die Gleichungen (50) zwei gekoppelte Integralgleichungen fiir die unbekannten Funk- 
tionen g, und m, dar. Bestimmt man mit diesen Lésungen die rechten Seiten von Gl. (51), so 
ergeben sich entsprechend zwei gekoppelte Integralgleichungen fiir die Funktionen g, und My. 
Die Gleichungen (50) sind identisch mit den Gleichungen (26) und (27) der ebenen Theorie. Die 
Gleichungen (51) sind ganz genauso aufgebaut, lediglich die rechten Seiten sind gedndert. Man 
erhalt aus den ersten beiden Gleichungen (50) die Funktionen g,(y) und m,(y). Mit diesen errechnet 
man die rechten Seiten der Gleichungen (51) und kann dann im zweiten Rechenschritt die Funk- 
i und m,(y) bestimmen. 

Tee ds aN rae zu lésen, begniigt man sich in der Praxis damit, die gesuchten 
Funktionen nur an einer endlichen Zahl diskreter Stellen m zu ermitteln, d.h. die ee 
Strémungsbedingung wird statt langs der ganzen Spannweite nur an m diskreten Punkten erfullt. 
Man erhalt dann statt der beiden gekoppelten Integralgleichungen ein lineares Re cere aes 
fiir die 2m Unbekannten g,,, m,,(v = 1, 2,...m) baw. gy, me, (vy = 1,2,... m). Wie schon ei 
der Darstellung der linearen Theorie erwahnt, erhalt man dann statt (50) das lineare ers 
system (28) und (29). Ganz analog ist das Gleichungssystem fiir die go, my,(v=—1,2,...m 
aufgebaut. — 
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c) Berechnung der aerodynamischen Beiwerte. Aus den ermittelten Funktionen 
21> Zo, M,, m, erhalt man fiir die aerodynamischen Beiwerte 


d 5 __ [dey ate 
C4>= ( ae + Cy a a” 5 Cy = (at) _ == C, aX (52) 


[a 


die folgenden Formeln: 


dc 
= A | &1(y) dy , 
=a 
41 ar re 
(a) =—A | ly Pan + 4 | men) 
a=0 A 
re ie (53) 


41 
Ga [eo dy, 
sap 


+1 +1 
I 
Co —A f gal) Pd + A [mln Pdr 
==] —l 


Auch hier lassen sich die Integrale wieder durch geeignete Quadraturformeln ausdriicken [vgl. (34) 
und (35)]. Fiir den induzierten Widerstand gilt bei scharfer Vorderkante nach (2) 


Cy; = C40. (54) 


6. Beispielrechnungen und Vergleich von Theorie und Experiment. Nach dem geschilderten 
Verfahren wurden fiir einige Fliigelserien die Auftriebsverteilungen sowie die aerodynamischen 
Beiwerte berechnet und, soweit méglich, mit MeBergebnissen verglichen. Fiir den Pfeilfliigel mit 


fico 


acy 
“\dorfy-9 


ai 


0 7 a 3 4 ro 0 7 2 a 4 5 
Abb. 7. Lineare Beiwerte fiir Auftrieb und Kippmoment in Abhiangigkeit vom Seitenverhiltnis bei verschiedenen 
Fligelformen. — Rechteckfliigel; — — — — Pfeilfliigel, p = 45° 5 7) ee ee Deltafliigel, 2 = 0,125. 


dem Seitenverhaltnis A = 1, dem Zuspitzungsverhaltnis 4 = 1 und dem Pfeilwinkel y = 45° ist 
als Beispiel in Abschnitt 8 die Zahlenrechnung ausfiihrlich wiedergegeben. 

Fiir Rechteckfliigel, fiir Pfeilfliigel konstanter Tiefe mit dem Pfeilwinkel y = 45° und fir 
Deltafliigel mit 2 = 0,125 Zuspitzung wurden bei verschiedenen Seitenverhaltnissen Rechnungen 
nach dem dargestellten Verfahren durchgefiihrt. Abb.7 zeigt die linearen Anteile (de4/da).—0 
und (dey/da),—o in Abhangigkeit vom Seitenverhaltnis. Fiir kleineres Seitenverhaltnis gehen alle 
de4/dx-Kurven in den Verlauf der Theorie schlanker Kérper de [da = (7/2) A titber. Die linearen 
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a wachsen mit zunehmendem Seitenverhaltnis. Anders die nichtlinearen Anteile, die fiir die- 
= ben Fligelserien in Abb. 8 als Funktion des Seitenverhaltnisses dargestellt sind. Die beiden 
oeffizienten der nichtlinearen Glieder C, und C, haben beide, absolut genommen, ein Maximum 


C. 
a 


0 7 é 3 4 5 0 ee lee 2 3 4 a 


Abb. 8. Nichtlineare Beiwerte fiir Auftrieb und Kippmoment in Abhangigkeit vom Seitenverhiltnis bei verschiedenen 
Fliigelformen nach Bild 7. 


<4) a dcy a 
c4 = | —— a + C, —a?, 35) (ee SRY Geet Mee 
( pi i eli “M (3 ie ioe [a] 
14 
G8 == 
° * 
06 sae 
nichtlineare. 
Theorie 


2° 2 b OP 


: 
a4 


cs 4 6 2° 1° 2° 2 do 02 OF G6 08 @ 472 
Abb. 9. Aerodynamische Beiwerte fiir einen Rechteckfligel vom Seitenverhaltnis A = 0,5 mit scharfer Vorderkante. Ver- 
gleich von Theorie und Messung. — a), b), ec) Aerodynamische Beiwerte ¢ ,, cy)» Cy in Abhingigkeit vom Anstellwinkel a; 


d) Auftriebsverteilung fiir = Ce = So lineare Theorie nach E. Truckenbrodt; —————- nichtlineare Theorie; 
@ Messung von NV. Scholz1, Re = 10°; © Messung von H. Winter (s. FuBnote 3 von S. 432), Re = 6.105. 


1 N, Scholz, Ing.-Arch. 18 (1950), S, 84 und Forschung Ing.-Wes. 16 (1949/50), S. 85. 


bei A = 0 und nehmen dann mit wachsendem Seitenverhaltnis sehr rasch nach Null hin ab. Damit 
liegt erstmalig eine Aussage dariiber vor, in welcher Weise die Nichtlinearitat mit wachsendem 
Seitenverhaltnis abklingt. Alle bisherigen Untersuchungen hatten dariiber keinen AufschluB 


geben kénnen. 
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nichtlineare 
Theorie \ o 
+ 


lineare 
Theorie 


c °° 4° 8° ~~? is? dC CkSCGSC«iBSC‘<~‘aS*<C«*SC | 
Abb. 10. Aerodynamische Beiwerte fiir einen Rechteckfliigel vom Seitenverhaltnis A = 1,0 mit scharfer Vorderkante. Ver- 
gleich von Theorie und Messung. — a), b), c) Aerodynamische Beiwerte ¢ , C445 ¢py in Abhangigkeit vom Anstellwinkel a; 


d) Auftriebsverteilung fiira = 19:40 Sere pee toe lineare Theorie nach E. Truckenbrodt; al nichtlineare Theo- 
rie; @ Messung von N. Scholz (s. FuBnote 1 von S. 445), Re = 8.105; OQ Messung von H. Winter (s. FuBnote 3 | 
ven §, 432), Re = 4,2.10°; A Messung von O. Flachsbart', Re = 6-10°. | 


Cy 


oe nichtlineare Theorie. 
42 f 


I Vireare 
Theorie 


lineare Theorie 


0 


Abb. 12. Theoretische Auftriebsverteilung fiir einen Pfeilfligel 

vom Seitenverhiltnis A = 1(y = 45°, 4 = 1,0) mit scharfer Vor- 

derkante nach der linearen und nichtlinearen Tragflachentheorie 

(« = 25°), — ———-— lineare Theorie nach E. Truckenbrodt; 
nichtlineare Theorie. 


Abb. 11. Beiwerte fir Auftrieb c "4 und Kippmoment c,, in Ab- 


hangigkeit vom Anstellwinkel « fiir einen Pfeilfiiigel vom Seiten- 
verhiltnis A = 1 (py = 45°, 4= 1,0) mit scharfer Vorderkante. 
Vergleich von Theorie und Messung. — —-—  —— lineare H 
Theorie nach E. Truckenbrodt?; ————_ nichtlineare Theorie; ah 

© eigene Messung, Re = 7.105 


* O. Flachsbart, Messungen an ebenen und gewélbten Platten; Ergebn. der AVA Gottingen, IV. Liefg. (1932). S. 96. 
2 E. Truckenbrodt, Z. Flugwiss. 2 (1954), S. 184, 
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In den Abb. 9 und 10 sind fiir die Rechteckfliigel mit den Seitenverhiltnissen A = 0,5 und 1,0 
Vergleiche mit Messungen an vorn zugescharften symmetrischen Platten durchgefiihrt wordag 
Dargestellt sind jeweils Auftrieb, Kippmoment und induzierter Widerstand in Abhangigkeit a 
Anstellwinkel a. Zum Vergleich wurden jeweils auch die Kurven der linearen Theorie mit ein- 
getragen. In allen Fallen ist die Ubereinstimmung zwischen der vorliegenden Theorie und den 
MeBergebnissen sehr gut. Fir alle drei Komponenten wird die Abhangigkeit vom Anstellwinkel 
bis zu hohen Werten desselben richtig wiedergegeben. 

In den beiden Abb. 9 und 10 sind dariiber hinaus fiir einen jeweils festen Anstellwinkel (~ = 14,7° 
bzw. 19,4°) die Verteilung der értlichen Auftriebsbeiwerte c,(y) ttber der Spannweite mit aur 
sprechenden Ergebnissen aus Druckverteilungsmessungen verglichen. Auch hierbei zeigt die 
nichtlineare Theorie wesentlich bessere Ubereinstimmung mit den MeBergebnissen als die lineare 


=| 
lineare 
Theorie 


0° ge 6° fe° 16° 20° 24° 28° 
Abb. 13. Aerodynamische Beiwerte Cy ey © fir Auftrieb, Kippmoment und Widerstand in Abhangigkeit vom Anstellwin- 


kel « fiir einen ‘Deltafliigel vom Seitenverhaltnis A = 0,78 (A = 0,125) mit scharfer Vorderkante. — Vergleich von Theorie 
und Messung.— —-—-— — lineare Theorie nach E. Truckenbrodt; —————— nichtlineare Theorie; © eigene Messung, 
Re = 910°. 


Theorie. Vor allem ist die theoretische Kurve ebenso wie die Messung iiber einen groBen Teil der 
Spannweite praktisch konstant und fallt erst nahe am Fliigelende plétzlich auf Null ab. Das aber 
bedeutet, da®B iber einen groBen Teil des Fliigels keine freien Wirbel vom Fliigel abgehen. Das 
dem Rechenverfahren zugrunde liegende Wirbelmodell liefert also in Ubereinstimmung mit der 
Beobachtung das Ergebnis, daB die freien Wirbel in der Hauptsache von den Seitenkanten abgehen, 
wo die Auftriebsverteilung sehr rasch auf Null absinkt und daher starke Anderungen aufweist. 

Einen entsprechenden Vergleich der theoretischen Werte fiir Auftrieb und Kippmoment mit 
den Experimenten fiir einen Pfeilfliigel (A = 1 und 2 = 1) zeigt Abb. 11. Auch hier besteht sehr 
gute Ubereinstimmung. Da hier keine konische Strémung vorliegt, waren solche Fliigelformen 
nichtlinearen Theorien bisher iiberhaupt nicht zuganglich. 

In Abb. 12 ist fiir denselben Fliigel die Verteilung der drtlichen Auftriebsbeiwerte c, auf- 
_ getragen. Man erkennt einmal den sehr grofen Beitrag des nichtlinearen Anteiles. Zum anderen 
ergibt die Hinzunahme der nichtlinearen Glieder wieder praktisch konstante Zirkulationsverteilung 
iiber dem gréBten Teil der Spannweite mit starken Anderungen in Randnihe. Interessant ist der 
kleine Einbruch der Zirkulationsverteilung in der Fligelmitte, der ein Paar abgehender freier 
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Wirbel zur Folge hat, deren Drehsinn, wie in Abb. 12 skizziert, dem der tblichen Randwirbel 
entgegengesetzt ist. Tatsichlich konnte experimentell dieses Wirbelpaar bei den beiden Fliigeln. 
mittels einer Wirbelsonde nachgewiesen werden. Auch der Vergleich der theoretischen Ergebnisse 
fiir einen Deltafliigel (A = 0,78, A = 0,125) mit Messungen in Abb. 13 zeigt gute Ubereinstimmung. 

Eine grundsatzliche und wichtige Frage ist, inwieweit Verfahren zur Berechnung der Auftriebs- 
verteilung an Tragfliigeln den Einflu8 der Hinterkante auf die Strémung beriicksichtigen. Die 
nichtlinearen Verfahren z. B. von Brown-Michael! und Mangler?, setzen bei der Behandlung 
von Deltafliigeln konische Strémung voraus, d.h. auf jedem Strahl durch die Fligelspitze sind 
z. B. die Druckwerte Ap = (p,, — p,) als konstant angenommen. Der Kinflu8 der Hinterkante, 
insbesondere das Absinken der Druckdifferenz Ap auf Null, kann dabei naturgemaB nicht mit- 
erfaBt werden. Das hier beschriebene Berechnungsverfahren ist so angelegt, daB auf jeden Fall 
die Druckdifferenz an der Hinterkante in Ubereinstimmung mit dem experimentellen Befund 
auf Null absinkt. 


7. Berechnung der induzierten Geschwindigkeiten in der Nahe einer ebenen Wirbelschicht. 
Gegeben sei nach Abb. 6 ein elementarer Hufeisenwirbel der Breite dy, der in der x-y-Ebene liegt 
und die Starke J" besitzt. Befindet sich der gebundene Teil des Hufeisenwirbels in x’, y', so wird ~ 
nach dem Biot-Savartschen Gesetz im Punkt P(x, y, z) eine Geschwindigkeit induziert, die folgende 
Komponenten besitzt 3: 


a te 
ti ware to (55) 
i P 2(y—y')z a ee (y—y')?z x—x! ; ; 
pv iatig ea ( 1 R )+ GvV—yyte2 R | ay ° (56) 
PG eee eee I 
Me salgarrter (| te )— Gye | x 
mit 
R=V@—x?+ (yy P+2. (58) 


Reiht man nun viele solcher Elementar-Hufeisenwirbel aneinander, so daB ein Wirbelband ent- 
steht, dann ergeben sich durch Integration in y-Richtung von —s bis + s die Komponenten 
der induzierten Geschwindigkeit zu 


n(x, ¥,%) = a | Py) dy’, (59) 


+s 


a eas % iar = ¥) te (y—y’)z x—x'),, 
v(x, Y, 2) re LO ee (1 R epaerpneer R dy , (60) 


+s 
Bee / C7) = x— x! zt x—x! 1 
w(x, ¥; 2) a. an | Cp eae ( ci R pees R | ay re (61) 


5 


Solange man die induzierte Geschwindigkeit in Punkten auBerhalb der Wirbelschicht berechnen 
will, bereitet die Integration iiber die Fligelspannweite keine Schwierigkeiten, Fiir Punkte in der 
Wirbelschicht (z = 0) verschwindet die Geschwindigkeitskomponente u, bei den Komponenten v 
und w treten jedoch an der Stelle y’ = y Unendlichkeitsstellen im Integranden auf, die eine ge- 
sonderte Betrachtung erfordern. Zur Durchfiithrung der Integration wird das Integrationsintervall 
in drei Teilintervalle zerlegt, so daB man fiir das Integral schreiben kann 


+s y—e yte s 


ford = fev dy + fund t [-nay, (62) 


—s —s y—e yte 
1 Siehe FuBnote 4 von Seite 435. 


2 Siehe Fu8note 8 von Seite 435. 
* E. Truckenbrodt, Z. angew. Math. Mech. 33 (1953), S. 165. 
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oder mit der Abkiirzung 
r ye +s 
T dy! = [ody + f ody (63) 
—s —s Yrs 
auch 
+s yte +s 
frvdt= fea tt.g, (64) 


aes y—e 


Dabei ist ¢ zunachst eine kleine, aber endliche positive GréBe. 
Bei der Bestimmung der Geschwindigkeitskomponente v(x, y,%) in der Wirbelschicht kann 
jetzt fir den zweiten Anteil in (64) des Integrales der Grenzwert z—> 0 sofort ausgefiihrt werden, 


da hierbei keine Unendlichkeitsstelle auftritt. Da zim Zahler des Integranden steht, verschwindet 
der zweite Anteil, und man erhalt 


yte 
pos Ls ' 2(y—y')z ee! Oy )2: %—2' ae 

oe walm { TO) a (1 fete cena RS ty" 2) 

ys 


Sowohl 2? als auch (y — y')? sind klein gegeniiber (x — x’)?, so da® R durch |x — x'| ersetzt werden 
kann. Aus (65) ergibt sich dann ; 


yte ye 
hee ers %— 2! 2(y —y')z dy’ a—x' ( (y—y')zdy’ 
59> 0) = —— a Neamt, 
v(x afi ) ee (i Se ah | I(y) ly —y')? + Fal bl os |x — x’|8 (y—y’)? + -| . (66) 


4 ame y—e 
Wahrend das zweite Integral fiir beliebiges z und ¢ verschwindet, erhalt man aus dem ersten Inte- 
gral durch partielle Integration 


y+e 
ee ae x—x’'\,. ar z dy’ 
v(x, ¥, 0) 3g re +e =—x) im | dy’ (y—y’? ae 2° (67) 


y¥—e 
Ersetzt man dJ’/dy' in dem Integrationsintervall durch einen Mittelwert, dann lat sich das Inte- 
gral sofort auflésen, und man erhalt im Grenzfall z > 0 
1 x—x'\ 2dr 
Hey.) —= (+e) aa 
Da in dem Ergebnis ¢ nicht mehr auftritt, kann ¢ beliebig klein gewahlt werden, so daB statt des 
Mittelwertes dI’/dy der Wert dJ'/dy an der Stelle y’ = y gesetzt werden kann. 
In derselben Weise erhalt man durch Aufteilen des Integrales nach (64) die Komponente w 
der induzierten Geschwindigkeit. Beriicksichtigt man die Beziehungen 


yte 2 dy’ 
tim | 5 ted a27) (69) 


(68) 


Oy) ek 


z—0 
y—s 
y+e 2 dy’ 
* spl oa eh 70 
Cp ee 
ed 


so erhalt man aus (61) 


y+e 
ers I(y’ x—x'\,, a I(y’) x—x\ 7, 
Dae ae i ae al! SR jay vag f v—y'? za(t oR )ay 
: nes (71) 
x— 
R 
y—exsy' <y+e, so laBt sich das Integral sofort losen, und man erhalt im Grenzfall z— 0 
+s 
x— x’ Ty’ x— x’ . 
w(e 740) = — gale Ty (i + =R=)— tat + Ro). (72) 


* \ wieder einen Mittelwert iiber dem Intervall 


Wahlt man im ersten Integral fiir [(y') (1 + 


ia a 


_Dabei ist jetzt 


R=Ve—2“¥+0—y) (73) 
Der Querstrich iiber dem ersten Term bedeutet die Mittelung im Intervall y—e Sy’ Sy +6. 
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Zur Bestimmung der z-Komponente w aufer- 


halb der Wirbelschicht, jedoch in der Nahe der- 
selben, geht man von einer Taylor-Reihenent- 


wicklung aus: 


(77) 


(<.°y, OF} Sot tee 


ow 
Oz 


w(x, ¥, z) = w(x, y, 0) + 


(78) 


ow 
5 ar eek 


ov 
oy 


a 


ou 
Ox 


Da die Strémung quellenfrei ist, gilt 
div B= 


0 und (68) 


Daraus folgt wegen u(x, y, 0) 
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Man erhalt daher fiir die vertikal induzierte Geschwindigkeit in Wirbelschichtnahe 


+s d ; x— x! 
w(x, ¥, 2) = —z— y—y’ 


|x — x’ } dy? |z| bis 


1 1 = ae! 5 z 
ily x(t ! x— x le (80) 


—s 


Besonders bemerkenswert ist, da8 das Zusatzglied unabhangig von der Koordinate x ist. 


8. Zahlenrechnung fiir einen Pfeilfliigel. Fiir den in Abb. 11 dargestellten Pfeilfliigel mit dem 
Seitenverhaltnis A = 1 soll die Berechnung der Auftriebsverteilung und der aerodynamischen 


Beiwerte im einzelnen durchgefiihrt werden. Der Fliigel ist durch folgende geometrische Daten 
gekennzeichnet: 


A= 1; p = 45°, as = 05, 

“alla he dx A 81 
x, = (sign y) y, yaad ny (81) 
a«*(n) =a** =a =konst.  (unverwundener Fliigel) . 


Ermittelt man die EinfluBfunktionen i,,, und j,,, an den Stellen Nn = Yn/s = cos [na/(m + 1)] 
(m = 15) auf der Hinterkante und der 1/4-Linie und benutzt man die von Truckenbrodt tabulierten 
Koeffizienten b,, c,, a,, und a,, fiir m= 15, so ergibt sich das erste Gleichungssystem (50) fiir 
die Werte g,,, und m,,, der Funktionen g, und m, an den Stellen 7, = y,/s in der folgenden Form: 


n=8 n=8 
Dy Gs Sin ae 2S Mis mM, = 1 (x ame 1, a fa 6) 16) - (82) 
ad! ni : 


Dabei bilden die Koeffizienten G,,,, und M,,,, eine 16-reihige quadratische Matrix, die in Tabelle 5 
wiedergegeben ist. 

Dieses Gleichungssystem wird am einfachsten nach dem von E. Truckenbrodt angegebenen 
Iterationsverfahren gelést. Die Tabelle6 enthalt neben den Lésungswerten selbst auch die 
ersten und zweiten Ableitungen, die nach den Summenformeln (45) und (46) ermittelt worden | 


Tabelle 6. Funktionswerte der Zirkulations- bzw. Momentenverteilung und deren Ableitungen fiir den Pfeilfliigel 
mit dem Seitenverhdltnis A = 1. 


din PZ1n dmyjn, Pm py, 

re if 81n Bune ra Es i in 82n Mon 

1 0,981 0,185 0,033 | —4,726 —128,28 —0,684 —21,05 3,826 | —1,056 
2 0,924 0,361 0,061 | —2,262 — 16,90 —0,382 — 4,50 2,897 | —0,827 
3 0,832 0,522 0,077 | —1,383 — 5,55 —0,022 — 1,43 2,012 | —0,583 
4 0,707 0,661 0,078 | —0,903 — 2,92 0 — 0,77 1,525 | —0,448 
-) 0,556 0,769 0,063 | —0,553 — 1,40 +0,183 — 0,25 1,137 | —0,334 
6 0,383 0,843 0,041 | —0,315 — 1,29 + 0,060 + 0,11 0,985 | —0,301 
a 0,195 0,879 0,010 | —0,077 — 1,00 +0,323 — 1,26 0,774 | —0,232 
8 0 0,883 —0,039 0 0 0 + 3,71 0,193 | —0,163 


sind. Mit diesen Werten lassen sich dann die EinfluBfunktionen f* und f}* nach (43) und (44) 
bestimmen und daraus die rechten Seiten des zweiten Gleichungssystems (51), welches lautet 


n=8 


n=8 
me; Gr ban Da les M2, — R,, “= 1,2...16). (83) 
n=1 n=1 


Dieses Gleichungssystem besitzt dasselbe Koeffizientenschema wie das erste System (Tabelle 5). 
Die rechten Seiten ergeben sich folgendermafen: Fiir x = 1 bis x = 8 gilt R, = Ry mit x =n 
und 


ts dg, n dg, n dm, n ate 84 
Rea — 0.75 pee ade (ede 8): (84) 


Fir x = 9 bis x = 16 gilt R, = R3* mit x = n + 8 und 
amin 


Tn + 1,655 Si pe, 2. 5 18) 


dy 
| (85) 


dg n dg nm 
Rist = — 0,103 A + 0,609 Fe — 0.216 Bin — 0,333 
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Die numerische Rechnung liefert dann: 


Die Lésungen des zweiten Gleichungssystems sind in Tabelle 6 zusammen mit den Lésungen 
des ersten Gleichungssystems tabuliert. Fiir den Anstellwinkel « = 25° sind daraus die Auftriebs- 
verteilungen nach der linearen und der nichtlinearen Theorie in Abb. 12 aufgetragen.. Die aero- 
dynamischen Beiwerte ergeben sich durch Integration nach (53). Berechnet man die Integrale 
mit Hilfe der Quadratuiformel (34), so erhalt man 


9. Zusammenfassung. Es wird ein Verfahren zur Berechnung der Auftriebsverteilung an 
Tragfliigeln in inkompressibler Strémung entwickelt, das iber den Rahmen einer linearen Theorie 
hinaus auch eine nichtlineare Abhangigkeit des Auftriebes, des induzierten Widerstandes und des 
Kippmomentes vom Anstellwinkel liefert. 

Durch eine einfache Erweiterung des iiblichen Wirbelmodelles der linearen Tragfliigeltheorie 
14Bt sich die nichtlineare Abhangigkeit der aerodynamischen Beiwerte vom Anstellwinkel erfassen. 
Fir die nichtlinearen Anteile, welche besonders bei Tragfliigeln mit kleinem Seitenverhaltnis 
gegeniiber den linearen Anteilen nicht mehr vernachlassigbar sind, liefert die Theorie die Ab- 
hangigkeit vom Seitenverhaltnis. Die Berechnung erfolgt mit Hilfe eines Gleichungssystems, 
welches aus dem entsprechenden System der linearen Theorie durch eine geringfiigige Abanderung 
hervorgeht. Die Theorie ist auf Fliigel mit beliebiger Grundri®form anwendbar und liefert nicht 
nur die Kraftbeiwerte, sondern auch die Auftriebsverteilung iiber der Fliigelflache. Fir kleine 
_ Anstellwinkel geht sie in die lineare Theorie itiber. — Vergleiche der theoretischen Ergebnisse mit 
Messungen zeigen sehr gute Ubereinstimmung. 


(Eingegangen am 2. November 1960.) 
Anschrift des Verfassers: Dozent Dr.-Ing. K. Gersten, Braunschweig, Bienroder Weg 3. 
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